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The book 1s a monograph written as a result of research by the 
author. The diffraction of plane electromagnetic waves by ideally 
conducting bodies, the surface of which have discontinuities, is in- 
vestizated in the book. The linear dimensions of the bodies are 
assumed to be large in comparison with the wavelength. The method 
developed in the book takes into account the perturbation of the field 
in the vicinity of the surface discontinuity and allows one to sub- 
stantially refine the approximations of geometric and physical optics. 
Expressions are found for the fringing field in the distant zone. 

A numerical calculation is performed of the scattering characteristics, 
and a comparison is made with the results of rigorous theory and with 
experiments. 


The book is intended fcr physicists and radio engineers who are 
interested in diffraction phenomena, and also for students of advanced 
courses and aspirants who are specializing in antennas and the 
propagation of radio waves. 
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FOREWORD 


First of all, one should explain the term “physical theory cf 
aiffraction". In order to do this, let us discuss brlefly the histo- 
rical development of diffraction theory. 


If one investigates, for example, the incidence of a plane elec- 
tromagnetic wave on a body which conducts well, all the dimensions of 
which are large in comparison with the wavelength, then the simplest 
solution of this problem may be obtained by means of geometric optics. 
It its known that in a number of cases one must add to geometric optics 
the laws of physical optics which are connected with the names of 
Huygens, Fresnel, Kirchhoff and Cotler. Physical optics uses, together 
with the field equations, She assumption that in the vicinity of a 
reflecting body geometric optics is valid. 


At the start of the Twentieth Century, a new division of mathe- 
matical physics appeared — the mathematical theory of diffraction. 
Using it, rigorous solutions to the problem of diffraction by a wedge, 
sphere, and infinite cylinder were obtained. Subsequently, other 
rigorous solutions were added; however, the total number of solutions 
was relatively small. For sufficiently short waves (in comparison 
with the dimensions of the body or other characteristic distances) 
these solutions, as a rule, are ineffective. Here the direct 
numerical methods also are unsuitable. 


Hence, an interest arose in approximation (asymptotic) methods 
which would allow one to investigate the diffraction of sufficiently 
short waves by various bodies, and would lead to more precise and 
reliable quantltative results than does geometric or physical optics. 
Obviously, these methods must in some way be ccnsidered the most 
important results extracted from the mathematical theory of diffraction. 
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tn the "reometric theory cf diffraction" proposed by Keller, the 
results obtained In the mathematical theory of diffraction of short 
waves were exactly the ones which were used and gencralized. Here, 
the concept of diffraction rays advanced to the forefront. This 
concept was expressed rather as a physical hypothesis and was nct 
suitable for representing the fleld in all of space: it was not 
usable where the formation of the diffraction field takes place (at 
the caustic, at the boundary of light and shadow, etc.). Here it is 
impocsible te talk about rays, and one must use a wave interpretation. 


What hac been said above makes it clear why a large number of 
works appeared in which the diffraction of short waves was investi- 
rated by other methods. Among those applied to reflecting bodies with 
abrupt surface discontinuities or with sharp edges (strip, disk, 
finite cylinder or cone, etc.) one should first of all mention the 
works of P. Ya. Ufimtsev. These works began to appear in print in 
1957, and it is on the basis cf them that this book was written. 

\ 

P. Ya. Ufimtsev studied the scattering characteristics by such 
bodies by taking into account, besides the currents being excited on 
the surface of the body according to the laws of geometric optics 
(the "uniform part of the current" according to his terminology), the 
additional currents arising in the vicinity of the edges or borders 
which have the character of edge waves and rapidly attenuate with in- 
creasing distance from the edge or border (the "nonuniform part of the 
rand the vadiation fleld created by the additional 
currents by comparing the edge or border with the edge of an infinite 


evurrent"), Cne may 


wedge or the border of a half-plane. In certain cases, one is 
obliged to consider the diffraction interaction of the various edges 
— that is, the tact that the wave created by one edge and propagated 
past another edre is diffracted by it (secondary diffraction). 


Such an approach to the diffraction of short waves has creat. 
physical visualizability and allows one to obtain rather simple 
approximation expressions for the field scattered by various metal 


bodies. This approach may be called the physical theory of diffraction. 
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This name fs applled to many works on the diffraction of short waves 
in which the mathematical difficulties are bypassed by means of physi- 
cal conslderations. 


It is clear that the physical theory of diffraction is a step 
forward in comparison with physical optics, which in general neglects 
the additional (edge) currents. The results obtained in this book 
show that with a given wavelength the physical theory of diffraction 
gives a better precision than physical optics, and with a given pre- 
cision the physical theory of diffraction allows one to advance into 
the lonrer wave region and, in particular, to obtain a number of 
results which are of interest for radar where the ratios of the dimen- 
sions of the bodies to the wavelength do not reach such large 
values as in optics. 


In addition, the physical theory of diffraction encompasses a 
number of interesting phenomena which are entirely foreign to physical 
optics. Thus, in a number of cases the additional currents give, not 
a small correction to the radiation field, but the main contribution 
to this fleld (see especially Chapters IV and V). If a plane wave 
is diffracted by a thin straight wire (a passive vibrator), then the 
additional current falls off very slowly as one goes further from the 
end of the wire. Therefore, the solution is obtained by summing the 
entire array of diffraction waves (secondary, tertiary, etc.) which 
successively arise as a consequence of the reflection of the currents 
from the ends of the wires. It has a resonance character. Thus, the 
problem of the scattering of the plane wave by a finite length wire 
which is a diffraction problem of a slightly unusual type is solved 
in Chapter VII. The resulting solution is applicable under the condi- 
tion that the diameter of the wire is small in comparison with the 
wavelength and length of the wire, and the ratio of the lenrth of the 
wire to the wavelength is arbitrary. 


The final equations which are derived in this book and are used 
for calculations are not asymptotic in the strict sense of the word. 
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“hepofore, it is natural to pose the question: in what way will the 
subsequent asymptotle equations differ from them when at last one 
obtains them in v.e mathematical theory of diffraction? One can say 
beforehand fia: the main term of the asymptotic expansicn will not, 

in the gveneral case, agree with the solution obtained on the basis of 
physical considerations: other (at a rule more complicated) slowly 
varying functions which determine the decay of the fields and currents 
ax one poee further from the edzes and-borders, and also the diffrac- 
tion interaction cf the edges and the shadowing of the edge waves will 
fivure in the main tema. However, the refinement of the slowly vary- 
‘tne functions in the expression for the diffraction field is not able 
to seriously influence the quantitative relationships. This 1s seen 
from a comparison of the results obtained in this book with calcula- 
tions based on rigorous theory and other approximation equations, and 
also with the results of measurements. 


The relationships obtained in this book also should help the 
development of asymptotic methods in the mathematical theory of dif- 


raction, since they suggest the character of the approximaticns and 
the structure of the desired solution. 


L. A. Vaynshteyn 
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INTRODUCTION 


In recent years, there has been a neticeable increase of interest 
in the diffraction of electromanetie waves by metal bodies cf complex 
shape. Such diffraction problems with a rigorous methematical formu- 
lation reduce to an interpretation of the wave equation or Maxwell 
equations with consideration of the boundary conditions on the body's 
surface. However, one cannot succeed in finding solutions in the 
case of actual bedles of a complicated configuration. This may be 
done only for bodies of the simplest feometric shape — such as an 
infinitely long cylinder, a sphere, a disk, ete. It turns out that 
the resulting solutions permit one to effectively calculate the dif- 
fraction field only under the condition that the wavelenrtn 1s larrer 
than, or comparable to, the finite dimensions of the body. In the 

"quasi-optical"case, when the wavelenzth is a great deal less than the 
dimensions of the body, the rigorous soiutions usually lose their 
practical value, and it is necessary to add to them laboricus and 
complicated asymptotic studies. Here, the numerical methods for the 
soiution of boundary value problems also become ineffective. There- 
fore, in the theory of diffraction the approximation methods which 
allow one to study the diffraction of sufficiently short waves by 
various bodfes acquire sreat importance. 


The field scattered by a given body may be calculated apprexi- 
mately by means of geometric optics laws (the reflection equations, 
see,for example [1-3]), from the principles of Huygens-Fresnel an! 
from the equations of Kirchhoff and Cotler [3-6]. 


The most common method of calculation in t’.e quasi-optic case 
is the principle of Huygens-Fresnel in the formulation of Kirehhoff 
and Cotler — the so-called physical optics approach. The essence 
of this method may be summarized as follows. 
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Let a plane electromagnetic wave fall on some ideally conducting 
body which Ls found in free space. In the physical optics approach, 
the surface current density whlch fs induced by this wave on the 
{rradtated part of the body's surface is (in the absolute system of 


units) equal to 
j=, [nl]. (A) 


where e«@ is the speed of light in a vacuum, n is the external normal 
to the body's surface, Hp is the mapnetic field of the incident wave. 
On the darkened side of the body the surface current is assumed to be 
equal to zero (jy? = 0). Equation (A) means that on each element. of 
the body's irradiated surface the same current is excited as on an 
ideally conducting surface of infinite dimensions tanrent to this 
element. The scattered field created by the current (A) is then 
found by means of Maxwell's equations. 


It is obvious that in reality the current induced on the body's 
surface will differ (as a consequence of the curve of the surface) 
from the current 5°. The precise expression for the surface current 
density has the form 


Ja HS we 
where jt is the surface density of the additional current which 
results from the curve of the surface. By the curve of the surface, 
we mean any of its deviations from an infinite plane (a smooth curve, 
a sharp bend, a bulge, a hole, etc.). If the ba ‘s canvex and 
smooth and its dimensions and radii of curvature are lnrre in compari- 
son with the wavelength, then the additional current is concentrated 
mainly in the vicinity of the boundary between the illuminated and 
shadowed parts of the bocy's surface. But if the body has an edge, 
bend, or point, then the additional current also artses near them. 
The additional current density Is comparable to the density 4° as a 
rule, only at distances of the order of a wavelength from the corres 
opouding edge, bend, or point. Thus, if the body's dimensicns 
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sleniflcantly exceed the wavelength, the additional currents occupy 
a comparatively small part of its surface. 


Since the current excited by the plane wave on an ideally cuon- 
ducting, surface 1c distributed uniformly over it (the absolute magni- 
tude of its surface density is constant) then the vector 3° may be 
called the “uniform" part of the surface current. The additional 
current jt which 1s caused by the curve of the body's surface we will 
henceforth call the "nonuniform" part of the current. In the physical 
optics approach, only the uniform part of the current is considered. 
Therefore, it is no wonder that in a number of cases it gives unsatis- 
factory results. For a more precise calculation, it is necessary to 
also take into account the nonuniform part of the current. 


In this book, the results of the author relating to the approxi- 
mation solution of diffraction problems are discussed and systematized. 
Essentially, these results were briefly discussed in a number of 
papers [7-14]. Roughly at the same time, the works of other authors 
devoted to similar problems appeared. We will discuss them in more 
detail (in §25) after the reader becomes accustomed to the concepts 
being used in diffraction problems of this type. For the present, let 
us only note that in these works, as a rule, other methods are used. 


In the book, problems of the diffraction of plane electromapmnetic 
waves by complex metal bodies, the surfaces of which have discontinul- 
ties (edges), are investigated. The dimensions cf the bodies are 
assumed to be large in comparison with the wavelength, and their 
surface is assumed to be ideally conducting. 


Obviously, if the edges are sufficlently far from one another, 
then the current flowing on a small element of the body's surfacc in 
the vicinity of its discontinuity may be approximately considered to 
be the same as on a corresponding infinite dihedral angle (a wedre). 
In fact, in Chapter I it its shown (see also [5] §20) that the nonuni- 
form part of the current on a wedge has the character of an cdre 
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wave which rapidly decreases with the distance from the edge. 
Therefore, one may consider that the nonuniform part of the current 

is concentwvated mainly in the vicinity of the discontinuity. By means 
of this physically obvious assumption, the field scattered by a strip 
(Chapter I), by a disk (Chapter II), by a finite length cylinder 
(Chapter III) and by certain other budies of rotation (Chapter IV) 

1s calculated. 


For a more precise calculation, however, it is necessary to keep 
in mind that the actual current distribution in the vicinity of the 
body's edges differs from the current distribution near the edge of 
the wedge. Actually, the edge wave corresponding to the nonuniform 
part of the current, propagated along the body's surface, reaches the 
adjacent edge and undergoes diffraction by it, exciting secondary edge 
waves. The latter in turn produce new edge waves, etc. If all the 
dimensions of the body are large in comparison with the wavelength, 
then as a rule it is sufficient to consider only the secondary dif- 
fraction. This phenomenon is siudied in Chapter V using the example 
of a strip and disk. 


In the case of a narrow cylindrical conductor of finite length, 
the edge waves of the current decrease very slowly with the distance 
from each end. Therefore, here it is impossible to limit oneself to 
a consideration only of secondary diffraction, and it is necessary to 
investigate the multiple diffraction of edge waves. Chapter VII is 
devoted to this problem. 


The uniform and nonuniform parts of the current are more than 
auxiliary concepts which are useful in solving a@iffraction problems. 
In Chapter VI it is shown that one is able experimentally to separate 
from the total fringing field that part of it which is created by 
the nonuniform part of the current. There, it is also shown that the 
depolarization phenomenon of the reflected signal ts caused only by 
the nonuniform part of the current. 
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Let us note the following feature of the method discussed in the 
book. A physical representation of the nonuniform part of the current 
is widely used in the book, but nowhere are its explicit mathematical 
expressions cited. This part of the current is generally not expressed 
in terms of well-known functions. Obviously a direct integration of 
the currents when calculating the fringing field is able to lead only 
to very complicated and immense equations. Therefore, we will find 
the fringing field created by the nonuniform part of tre current on 
the basis of indirect considerations without direct integration of 
it (see especially Chapters I = IV). 


The method by which the diffraction problems are solved in this 
bock may be briefly summarized as follows. We will seek an approxi- 
mate solution of the diffraction problem for a certain body by first 
having studied diffraction by its separate geometric elements. For 
example, for a finite cylinder such elements are: the lateral surface 
as part of an infinite cylindrical surface, each base as part of a 
plane, each section of the base rim as the edge of a wedge (the curva- 
ture of the rim in the first approximation may be neglected}. Having 
studied the diffraction by the separate elements of the body, we will 
obtain a representation of the nonuniform part of the current and or 
the field which is radiated by it. Then secondary, tertiary, etc. 
diffraction is stuuled — that is, the diffraction interaction of the 
various elements of the body is taken into account. 


This method appeals to physical considerations, not only when 
formulating the problem but also in its solution process, and in this 
way differs from the methods of the mathematical theory of diffraction. 


Therefore, such a method may be referred to as the physical theory 
of diffraction. 


A whole series of other diffraction studies which appeared In 
the last five to ten years also are able to relate to the physical 
theory of diffraction. The first work which contained the idea of 
the physical theory of diffraction is evidently the paper of 
Schwarzschild [15] which was published at the beginning of thls 
century and was devoted to diffraction by a slit. 
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One should note that approximate solutions of diffraction 
problems would be impossible without the use of the results obtained 
in the mathematical theory of diffraction. In particular, the 
rigorous solution to the problem of diffraction by a wedge which is 
attributed to Sommerfeld [16] is widely used in this book. In 
Chapter I this solution is obtained by another method. The works of 
Fok [17, 18] served as the starting point for numerous studies on 
aiffraction by smooth convex bodies. The rigorous solution of the 
problem of diffraction at the open end of a wave guide [19] revealed 
the mechanism for the formation of primary diffraction waves, and 
their shadowing by the opposite end of the wave guide. The rigorous 
theory as applied to a strip and disk allows us to examine the 
precision of the approximation theory (see Chapter V). 
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CHAPTER I 


DIFFRACTION BY A WEDGE 


As was already said in the Introduction, the ffleld scattered by 
a body may be investi;ated in the form of the sum of the fields being 
radiated by the uniform and nonuniform parts of the surface current. 
The uniform part of the current 1s completely determined by the geo- 
metry of the body and the magnetic field of the incident wave. The 
nonuniform part generally is unknown. However, one may arproximately 
assume that in the vicinity of the discontinuity of a convex surface 
iv will be the same as on a corresponding wedge. Therefore, it is 
necessary for us to begin by studying the diffraction of a plane elec- 
tromagnetic wave by a wedgfe. This chapter will be devoted to this 
problem. First we will investipate the rirorous solution of this 
problem (§ 1 and 2). Then we will find its solution in the physical 
optics approach (§ 3). The difference of these solutions determines 
the field created by the nonuniform part of the current (§ 4). 


§ 1. The Rtgorous Solution 


The rigorous solution to the problem of diffractton of a plane 
wave by a wed;ye was first obtained by Sommerfeld by the methed of 
branching wave functions [16]. Later, the diffraction cf cylindrical 
and spherical waves by a wedge also was studied. A rather extenstve 
bfibllopraphy on these problems may be found, for example, in the 
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paper of Oberhettinger [20]. Since the problem of diffraction by a 
wedge lies at the base of our studies, we considered it advisable not 
only to present the results of its rigorous solution, but also to give 
them a new more graphic derivation. The idea for this derivation 
follows directly from the work of Sommerfeld. Sommerfeld found the 
solution to the problem in the form of a contour integral, and then 
he transformed it to a series. However, one may proceed in the oppo- 
site direction: first find the solution in the form of a series and 
then give its integral representation. Such a path seems to us more 
graphic, and is discussed in this section. The necessity for a 
detailed derivation is caused by the fact that the results of Sommer- 
feld [16] are not represented in a sufficiently clear form, which 
hinders their use. 


Let us assume there is in free space (a vacuum) an ideally con- 
ducting wedge and a cylindrical wave source Q parallel to its edge 
(Figure 1). Let us introduce the cylindrical coordinate system r, 9, 
z in such a way that the z axis coincides with the wedge edge, and 
the angle $ is measured from the irradiated surface. The external 
wedge angle will be designated by the letter a, so that O0<p%a. The 
coordinates of the source Q we will designate by ros y° 


Let us investigate two particular cases for the excitation of an 
electromagnetic field. In the first case, it is excited by a "fila- 
ment of electric current" 


j-== — iop,3(r — ty, 9 — 9%) (1.01) 
in the second case, it is excited by a "filament of magnetic current" 


7 2 
(= —iwm,3 (¢ — Fr, & — %)- (1.0?) 
The quantities Py and m, here designate, respectively, the electric 
and magnetic moments of the filament per unit length along the z axis, 
w is the cyclic frequency (och me), 8(r— ry 9-9) == 3(r--r 31 (@— gall 

is a two-dimensional delta function which satisfies the condition 
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Figure 1. The excitation of a wedge- 
shaped region by a linear source. 


- source; P - the observation point; 
L=- the integration contour in Equation (1.10). 


ffie—r. 9 — %) rdrde=1 
with integration over the neighborhood of the point Tg» 9° 
Here and henceforth, we will use the absolute system of units 


(the Gauss system), and we will assume the dependence on time is in 
the form e~™ 


In the first case, the "electric" vector potential Ay satisfies 
the equation (see, for example, [4]) 


= —!*_ 1.0 
AA’ + kates — "2? (1.03) 


and the boundary condition 


Als with p=0 and ¢--4 (1.04) 
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In the second case, the "marnetic" vector potential a satisfles the 
equation 


tr 5 
AA" + h*A™ xa -> r (1.05) 
and the boundary condition 
oat ; 
“jg 70 with g==0 and g=a. (1.06) 


It is natural to seek the solution of the nonhomogeneous 
Equattons (1.03) and (1.05) in the form 


ad, (&r) HO (hry) sin’y, 9ysin 49 WIth l< for 
'? '* 


v a,J, (Ary) H'" (kr) sin v.99 sin ve withe > foi (1.07) 


( 5 6,J, (kr) HY (Rrq) COS ¥,Fq COS YP WIth FS Fos 
Ac= | a0 : 
s oo 
ENT B,J, (Aro) HH" (kr) cos v9) COS v9 WIth 7 >to (1.08) 
| +20 °'¢ e 
= $ ile é 
The products - a 
J, (kr) sin v9 J, (kr) cos vp 
and 
HY (Ar) sin ¥,9 HY ke) COS ¥,9 (1.09) 


are the partial solutions of Equations (1.03) and (1.05) without the 
right-hand member which sattsfy the boundary conditions (1.04) and 
(1.06). The remalning factors entering into Equations (1.07) and (1.08) 
ensure the observance of the reciprocity principle and the continuity 


of the fleld on the arc r 2 To: The Bessel function a (kr) enters 


these equations when r < To because it remains finite when r » 0, 


a 


and the Hankel function HY (ke) 1s taken when r > rg in order that the 


solution satisfies the radiation condition. 
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The coefflelents a. and by may be determined by means of Green's 
theorem 


2 Og 
ea =: [ auus, dS =: rdrdy (1.10) 
ul . 

for the contour L In the plane z = const which 1s shown in Figure 1. 
Bere, the external normal to tie contour L is designated by the letter 
n. Applying Equation (1.10) to the functions AY and AD and performing 


the limiting transitions r) * ry and r, + ry in it, we obtain 


( Pomel” ON 


eae 


re ele git mn, fa(e— —~ Fa) “ef? 


dekh 
: ) rye =i oe ps f¥e—sdrede, 








Since here the integration limits are arbitrary, 1t follows from the 
equality of the intearals that the integrands are equal: 





ong _? Ag =A8 | re ee 
Sarg Sa wr Kaa sq— —%), ( 1 ° 11 ) 
0A vas 4ckm, 
1 atts (1.12) 


Now let us substitute Expressions (1.07) into Equality (1.11) and 
multiply both members of the latter by sinvyg . Then tnterratine the 
resulting equality over ¢ in the limits from 0 to a, we find 

a=" kp. (1.13) 
In a similar way, let us determine the coefficients 


by == eyo hime (1.14) 


where 


w=. A eo (1.15) 


PT D\eHCH3 3-25-71 


wn 


Consequently, tne electric current filament excites, in the space 
outside the wedge, the fleld 


t @ 
. 423 : . . 
i~— tps y) nH” (kroy J, (kr) sin veg sin vo 
aml 


with r<f 


r- ri 
tS aD Shr) HM (hr) sia ggsin vg 
axl 
with r>% 
E,=E,=0, H=-1 rot E, (1.16) 


and the magnetic current filament excites, outside the wedge, the 
field 


= ikAT= 
4a? - 
i- #n, y “Hy (krg)J, (41) 60S »4F9 C08 4 
an0 


with r<lo 


= 


) 
i as k’m, y td, (kr) Hn (kr) COS ¥, Fy COS ¥4P 
ame 


with r>le 4 
1 
-H,-=H,=0, Em —-— roth. (1.17) 


Now using the asymptotic equation for the Hankel function when 
krg> CO [21], we have 


ifae— aoe 


HY (hr) =Vure ( 2 7 HO ergot » 1.13) 





Then Expressions (1.16) and (1.17) in the region r < r. take the form 


0 


12 1 Bp H" (hr) X 


hed —t s. * 
x y e 7's , (hr) sin v.99 sin v9, 


(equation continued on next page) 
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Hyatt ktm! (he) X 


oo 
te. ) 
x y s,e * 4, (kr) cos », 9) COS v9 


8:30 . 


or 

E, = ine’ p,H (kr,) [a(r, 9@— Pdl—4lr, P+). 

H,==ink'm,H) (hr) {u (r, = %) + u(r, e+). (1 . 19) 
where 

— ad -1% % 
u(r, =F Yue . 4, (&r) cos vp (1.20) 
s=f . i 

(p=9=®) > 


Let us note, furthermore, that in free space the field of the 
electric filament with a moment P, is determined by the relationship 


t 
E, = ink" p,H" (kr), (1.21) 


and the field of the magnetic filament with the moment m, is determined 
by the relationship 


H, = izk’m,H,'(kr). (1.22) 


Therefore, the expressions in front of the square brackets in Equations 
(1.19) may be regarded as the primary field of the filament — the 
cylindrical wave arriving at the wedge edge. Now removing the fila- 
ment of current to infinity (ro +o), let us proceed to the incident 
plane waves 


£,= E,g.e tt e890), 1 E,==E.=0 (1.23) 


and 
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Hy=Hgy-e ~sd. Hy == Hy 0. 7 (1.24 ) 


The field arisinzs with the diffraction of these waves by the wedge 
will obviously have the component 


E,==Ey,[4lr, 9--9,) —4 (0, 9-F PI (1.25) 


and 


H,=Halu(r, 9--%)+-4lr, 9+ %)]- (1.26) 
Let us find the integral representation for the function u (r, 
v). For this purpose, let us use the equation (see [16], p. 866) 


J, toma fe {« Con P+, (-F)ue 
1.27 


where the limits I - IIi mean that the integration contour goes from 
region I to region III (Figure 2). The cross-hatched sections in 

the plane of the complex variable 6 (8') shown in Firure 2 are regions 
in which Imcosjs>0 (Imcos#<0), Therefore, in the sections of the con- 
tour extending to infinity the intexrand strives to sero, ensuring 

the convergence of the intesral. Substituting Expression (1.27) into 
Equation (1.20), we obtain 


a(r, Y= 
mat 


moe ike cos 3 o iv, P~aty) \ iv, (3-8) dp. 
\s (az? . a y : J 


s=l 


After summing the Infinite seometric progressions and replacins the 





variable 8 by B' = 68 ~ m7, the function u (r, W) acquires the forn 
%) == 
we m9) 
=). f goiter cos Sasa eae a’, 
ae) 1h +e) -itw—w 
l—e l—e 
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Fivure 2. The invegration contours Figure 3. The integration contour 
in the complex plane 8. in Equation (1.28). 


As a result, we obtain the well-known integral of Sommerfeld 


a(r ») e— iar cos 8 
= | pee iba - (1.28) 


The integration contour C is shown in Figure 3 and consists of 
two infinite branches. Since the integrand expression has poles at 
the points Bxa=2em—% (m=0, +1, 2,...), then for the values of $ corre- 
sponding to the space outside the wedge (09<¢<a) the function u (r, p) 
may be represented (with t<as<2n, 0<9,<2 ) in the following way: 


u(r, Y=olr, eT" with—ecpce, 
u(r, ¥)=o(r, 9) With ec 2a—z, 
u(r, 9) =9(r, ) e fhe cos (2a—9) 

ve with 2a—s8<y<2a, (1.29) 


where 


ew lbrces |] 
o(r, non | Gamer APRON 
1 e 
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or 


at n 
‘ 608 = — cos (+9 


aera the cost 
ee en 


The integration contours D and Dy are shown, respectively, in Figures 
3 and 4. * 


With an arbitrary incidence of a plane wave on a wedge, one of 
two cases may occur: (1) the plane wave "illuminates" only one face of 
the wedge (0<9,<a—#) , and (2) the plane wave "illuminates" both 
faces of the wedge (2—*#<9,<2) . Let us write out in more detail the 
functions u (r, ¥) corresponding to these cases. In the case %#<2—* 
(Figure 5), we have 


u(t, 9 — F_) = } ’ 
=o/(r, eee ela con (%~ Zeb 
peas )=: binmeecesy, 
’ >) == 
=0(r, 9+ 9,)-fen i emierin | 


u(r, 9— P= 
=O (r, P— Py) + eit con; —¢,! \, ithe-- PoE Pes 
u(r, PA P)=Ilr, PT Po) 
u(r, 9 — Po) == 9(7, P— Po) \ 


< ° 
u(r, P49) == 0(r, P+ Fo) with t+%<79<e 


, (1.31) 


and in the case *4—-"*<%<e (Figure 6) we have 


u(r, 9-- P=. } 
=o(r, 9 — Fa) pe” ier corte te | 
withd tH, 
u(r, 9P+%) = thocp<cn--% 
s=0(r, a-+r)te Meret | 
a(t, 9--%) = 
( ? Pa) = C08 (¢ — Oe) | with 


s=9(r, 9--Pdte 


Re Pir PLLA—-— 9, 
u(r, PEP) =O(r, P-E Mo) J ; : 


OS itt 
=0(r, — — Pa) e' 7 60d (¢ —Ca) w t , 
oe poe |. (1.32) 
=u(F, 9-E a) fer : 
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The direction $¢ =m - 5 corresponds to the ray reflected in a specular 
fashion from the first face (we will consider as the first face that 
face from which the angles are calculated), and the direction = 

=s 2a—=— 9 corresponds to the ray reflected specularly from the second 
face (Figure 5 and 6). The functions e~"™** describe plane waves of 
unit amplitude: e~"s*"-~ describes the incident wave, e~ raiser 
describes the wave reflected from the first face, and gq ~treos(e-e-e) — 
the wave reflected from the second face. 





Figure 5. Diffraction of a plane 
wave by a wedge. The plane 


wave irradiates only one face of 
Figure 4. The integration contour the wedge. $, is the angle of 
in Equation (1.30). 


incidence. The line 6 = 7 - % 
1s the boundary of the reflected 


plane wave, and the line 
ee a 95 is the boundary of 


the shadow. 


Figure 6. Diffraction by a wedre. 
The plane wave irradiates both 
faces. The line 6 = 2a - 7 = oo 


is the boundary of the plane wave 
(ca from the second face 
62a). 
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§ 2. Asymptotic Expressions 


The integral 





t at ¢ gitr cos tte 
o(r, $) = 5, sla EL ee 
. | ene ters (2.01) 


in Equations (1.31) and (1.32) generally is not expressed in terms of 
well-known functions. However, when kr >> 1 it may be calculated 
approximately by the method of steepest descents [21]. In integral 
(2.01), changing for this purpose to a new integration variable 


sa fte™ sin =, s*=i (1 —cosQ), 


we obtain 


sia ‘ w+) entre 
0(7.9)= Fee ds ; 
SexStys con e0t (con $ con EE) con (2.02) 


where 
n=. (2.03) 


It is not difficult to see that the point s 2 0 is a saddle point: 

as one goes further from it along the imaginary axis (Re s = 0) in 
the plane of the complex variable s, the function e-"" most rapidly 
increases, and as one goes along the real axis (Im s = 0) it decreases 
most rapidly. Therefore, when kr >> 1 the main contribution to the 
integral (2.02) is given by the integrand in the section of the 
contour in the vicinity of the saddle point (s = 0). 


The method of steepest descents is carried out by expanding the 


est 


integrand (except for the factor ) into a Taylor series in powers 


of s. This series is then integrated term by term. If the integrand 


FTD~HC=23-259-71 12 


expansion converges only on part of the integration contour, the re- 
sultant series obtained after the integration will be semiconvergent 
(asymptotic). Limiting ourselves to the first term in it, we obtain: 








cop =~ cos — V2etr ; (2.04) 


The remaining terms of the asymptotic series have a value on the order 
1 


of ant and higher. 

Expression (2.04) is valid with the condition (cos = — cost.) / ir >1 
and describes that part of the diffraction field which has the charac- 
ter of cylindrical waves diverging from the wedge edge. With the 
incidence of the plane wave (1.23) on a wedge, the electric vector of 
which 1s parallel to the wedge edge, the cylindrical wave is determined 
in accordance with (1.25) and (2.04) by the equation 


£,=—H,=E,-[0(r.9—%)—9(r, Pt) = 








T a 
= P,,f ——--—, 
AS : ¥ 2akr (2.05) 
where 
j “*( 1 i 
® Von con FE con con EE (2.06) 


When the wedge is excited by the plane wave (1.24), in which the mag- 
netic vector 1s parallel to the wedge edge, the cylindrical wave has 
the form 


H, = E, =H lo(r.9— Wt etl = 


(a +) 
= Het ae 
‘ (2.07) 
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where 








sio — ; ‘ 
&=— a Le (2.08) 


In the vicinity of the shadow boundary (¢ unt 5) and near the 
directions of the mirror-reflected rays (p=*—,, 9 2a—*— ?e) 
Expressions (2.04) - (2.08) are not valid, since the poles 


12 
sin*>", s,=Vie “sin (ne —*3* 


= fie 


of the integrand in (2.02) are close to s = 0 and, consequently, its 
expansion in a Taylor series loses meaning. Physically, this result 
means that in the indicated region the diffraction wave does not re- 
duce to plane and cylindrical waves, but has a more complicated char- 
acter. An asymptotic representation of the function v (r, ¥) in this 
region was obtained in 1938 by Pauli [22]; here we will present the 
derivation of the first term of the asymptotic series obtained in [22]. 


Let us multiply and divide the integrand expression in Equation 
(2.02) by the quantity 


cos + cost =i (s* — is’) (si=2cos*-3-) (2.09) 


and let us expand into a Taylor series in powers of s the function 


cos $+ coa§ 
* con cone 


g° 
(cos eos A joo 
which no longer has a pole at the saddle point (s = 0) when wv = ¢ + 
+ 95 = 7. Limiting ourselves in this series to the first term, we 
obtain 





s*—ist 


o (rig Ve ttt aes a 


(2.10) 
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The integral here may be represented in the form 


ir enfe® ea ° ra - (ots) c 
-— ~~ ds=e ol ds J e dt. 
s*~is2 2 


Changing the order of integration here, we find 








es inet © et ato 
£ ds=e ° fe ae Ome a 
st—ig 2 vt 
Vu _—thrs? i) 
Tye) a (2.11) 
ae ed 
and finally 
s 4 eden 
ln cos “5 
2 a 2 arcade © 
er ie 
; _ 
x J e“dq. 
Vite |cos (2.12) 


The next term of the asymptotic expansion for the function v (r, yp) 
has a value, whose order of magnitude depends on the observation 
direction: in the vicinity of the border of the plane waves (7=#*+%) 
its order of magnitude is Ji » but far from it the order of magnitude 
is 1/kr in comparison with the term written in (2.12). 


It is convenient to represent Expression (2.12) in the following 
form: 


p= 
2 sin —- cos > bs van 4 
on g=—— re? thr cond, va x 
a a 
V2br cos — 
D4 j e'“dg. 
«cos (2.13) 
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Here the absolute value of the lower limit of the Fresnel integral 
always equals infinity, and its sign is determined by the sign of 
cos ~/2. Therefore, when passing through the boundary of the plane 
waves (p = 4 + 5 = 1) the lower limit changes sign and the Fresnel 
integral undergoes a finite discontinuity, ensuring at this boundary 
the continuity of the function u (r, ¥) and consequently of the 
diffraction field. Actually, by means of the well-known equation 


@ —— a 
fetaq=> Ry 4 


F (2.14) 
it is not difficult to show that 
ofr, e+ 0)=, o(r,a— = — 2 (2.15) 
and consequently 
a (r,n-20)=+e™, (2.16) 
In view of the asymptotic relationships 
felag=4r., { edg = —S (with p> 1) 
: . (2.17) 


Pauli's Equation (2.13) 1s transformed with ¥ 2k cos 3] >1 to the 
Expression (2.04). As was already indicated above, it determines the 
cylindrical waves diverging from the wedge edge. 


By means of Equation (2.13), one may also calculate the field in 
the vicinity of the direction 6 = 2a - n - $6 —that is, near the 
boundary of the plane wave reflected from the face 4 = a; for this 
purpose, it is sufficient to replace $ by a - ¢ and $5 by a - o>: 


It is also interesting to note that in the case of a half-plane 
(n = 2) Equation (2.13) gives the expression 
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~i > Vite co $ 
o(r, y) =e ome = x e'“dg, 
| wen $ (2.18) 


which completely agrees with the rigorous solution. Actually, with 
a = 2m, when the wedge is transformed to a half-plane, integrai (1.30) 
equals 


i gltr cost 
v7, p= — ye see 
cane a (2.19) 


and it may te reduced to a Fresnel integral. For this purpose, let 
us divide the contour Do into two parts by the point > = 0. Summing 
the integrals over these parts of the contour, we find that 


+ i0 


d=——a a ocr 
a a 





oa 


— 


Now changing to a new integration variable s=/Ze * sin 5 and taking 


into account Equation (2.09), we obtain 


paar ents? 
eee aa Sema (2.20) 


The integral here was already calculated by us. Turning to Equation 
(2.11), we arrive at Expression (2.18) which — together with rela- 
tionships (1.25), (1.26) and (1.31) — give us the rigorous solution 
to the problem of the diffraction of plare waves by in ideally 
conducting half-plane. 
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§ 3. The Physical Optics Approach 


In the physical optics approach, the fringing field is sought 
as the electro-magnetic field created by the uniform part of the sur~ 
face current 


y=,~[nH,)- (3.01) 
Let us recall that here n designates the external normal to the body's 
surface, and Hy designates the magnetic vector of the incident wave. 
First let us investigate the case 0:°.4,<a--7, when the incident plane 
wave irradiates only one face of the wedge (Figure 5). 


From Equation (3.01) it follows that the density of the uniform 
part of the current being excited on the irradiated face by plane 
waves (1.23) and (1.24) has the following components, respectively 


fo = Bys-sin pee", p=fi=0 (3.02) 


and 


—thky . 20 of 
V C03 Fo 0 


i= a Hee i, =i, =0. (3.03) 
For the purpose of calculating the field radiated by this current, we 
will use the following integral representation of the Hankel function 
(see [16], p. 866 ) 


' Riz a 
Hy” (p)= = j phen ge es ( et ch dt 
~t+los a 


(0<2: 8S x). (3.04) 


Assuming here p = kd and changing to a new integration variable 
Z =dsht, we obtain 
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H” (dye (OEE ae 
u i* : Vd pet (3.05) 


It is easy to show by means of Equations (3.02), (3.03) and (3.05) 


that the vector potential 


[-") @ ee ae 
1 7 ei* Ve ARUP 
A(x,y,0)-2 © ( jocna: ie Ae ee ee acde 
‘ 5 Vu Fae +H (3.06) 


has the components 
Av=-p Eorsinte I, As=A,=0, (3.07) 
if the wedge is excited by plane wave (1.23), and 


A= Hel, A, =4,=0, (3.08) 


if the wedge is excited by plane wave (1.24). Here, I, designates 
the integral 


@ 
ieee (UG Vy +(x — 8 at. (3.09) 


Let us transform it by using the relationship 





d>0), (3.10) : 


It is not difficult to establish the correctness of this relattonship 
by verifying that it changes into Expression (3.04) with the substi- 
tution w= k sin t, v 2 k cos t and kpd'fziz=p. As a result 


o 

t ef (o lul--wa) 
I, { p(k cos gy, — a) = (3.11) 
—0 
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where the integration contour passes above the pole w = k cos oo: 


Let us note that integral (3.11) is a function of |y|, and let 


us change to polar coordinates according to tne equations 


X= 70S 9, 
lyj=rsing with ¢<, 
lwi=—rsinewith e>*. 


Furthermore, by carrying out the substitution 
w==—k cost {o==ksin6), 


we obtain 


1 faite cos Ro) 
ia remit withecs, | 
4 giltr cos (te) 
t= inp! COS 99+ cose 
y . 


dG withpe>a, 


The integration contour F is shown in Figures 7a and 7b. 


(3.12) 


(3.13) 


(3.14) 


In Figure 7a 


the cross-hatched areas indicate the sections in the plane of the com- 
plex variable £ in which Imcosi—9>0 ; in Figure 7b, the cross-hatched 
areas indicate the sections where Incos(t-+.9)>0 . Now let us deform the 
contour F into the contour G (G,) for the values ?<*(7>r), and let 


us change to a new integration variable 


C=t—» withe<z, 
C=3§—(22—g)with p>. 


As a result, we obtain the following expressions: 


Owithgy>r—», 


1 efrcos lye \ 
=i | cos yPeoseby tT gt itr cos ie +e0) 
'e hain yy 
with PEP, 


if ¢@ < 7 and ‘ 


FPTD-HC-23-259-71 20 


(3.16) 


x 


4 
Yi; 
GH 

2 





Figure 7. The integration contour in 
Equations (3.14), 


Owith pete, 


po: l elas cos Cae 
ve inh | cosy, + cos G&—¥) So 2 en itr cong — ey 


A sing, 
with? >=z+49,, 
(3.17) 
1f o>. The integration contour Dy is shown in Figure 4, 
By means of Equation (3.06) and the equality 
E,=tkA, (3.18) 


let us find the field which is radiated by the uniform part of the 
current excited on the face ¢ = 0 by the plane wave (1.23) 


0 (9, 9.) — enlir Fonte tee) with0<pce_p, 
£, | % ( %) with®—~a%<e<e, 
Es | OF (p %) withe<pcet a, 


O (% Me) — OO With Re Coa, (3.19) 
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where 


e! de cos cae 


+ _i. 
mae massing | ieee (3.20) 


It is not dirficult to see that with a—7z<9,<« , when both faces of 
the wedge are illuminated, the field excited by the second face 


(¢ = a) will be described by the same equations if one replaces ¢ by 
a-=- $, and %5 by a - 5 in them. 


Adding the field being radiated by the uniform part of the 
current with the incident plane wave (1.23), we obtain the diffraction 
field in the physical optics approach. It equals 


of (?, Pe ) + ei con(¢—te)__ ent 08 (¢+¢—) 


with0«p<e—y,, 
O (72 Me) perv cmn— v4 


om with §—%<9<#, 

"Torte dpe 

with e<e- ety, 

oF (9. Ge) with §+%<?9<% (3.21) 


if one face of the wedge (0<%&<e—*) is illuminated, and 


OF (F %)+o, (2—%, 2—H%)+ 
+ the cos (¢ —9e)__ e ihr 003 (9 > %e) 
with 0<9<e—p, 
OF (7, %)-+ 0) (@— 9, 2—%)+ 
fevered with r—gceca—e, 
E, | Of (ts %) tot (2— 9, 2%) + 
pe ieee with a—acece, 
OF ('P, Pe) + 0; (2— 9, 2— Fe) + 
pe Preto with ac gc 2a—2— Py 


OF (» Fe) +9; (a— 9, 2%) 
a ei" cos(¢ — Fe) enitr cos (2a— 9 ~ FQ 


cy 


with 2a—"2—@%<pas, (3.22) 
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1f both faces of the wedye (2--4r<9,< 2) ave illuminated. 


Now let us calculate the field arising during diffraction by a 
wedge of a plane wave (1.24). The field seattered Ly the first face 
(9 = 0) 1s determined by the relationship 


One may write the component Hy in the form 


Hi=—~FHad I, (3.24) 


* ay 


or 


He=— eal, with P<, 


He= heh vith e>e (3.25) 


The quantity I, introduced here is the integral 
— de 
kco3 9, —-@ (3.26) 


1 
L=-s 


along the infinite contour which passes above the pole w 2 k cos o- 
This integral, precisely the same as integral I); is transformed to 
an integral along the contour Do: As a result, we obtain 


OF (9, Pe) pen re tO PCE — PH 


on - Oy (% %) with *—%w<er<e, 
Hee vy (9, Pe) with r< P< rte, 
e (Ge) — eth wh ge FcR, (3.27) 


where 


‘ i (cin CS pei areont 
Oy (> Po) = = a cosyetcosié ty 
. : 
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In the case when both faces of the wedge are illuminated, the field 
scattered by the second face also is determined by Equations (3.27) 
and (3.28) in which one need only replace $y by a - oo, and ¢ by 


a- >. 


Then adding the field radiated by the nonuniform part of the 
current with the incident wave (1.24), we find the diffraction field 
in the physical optics approach. This equals 


OF (9, G4) fe em OIF ~ PL. = UT onlt #90 

with 0ceer— Pe 
oy (9 %) en itr cos ¢-%) 

ag with *#—%®<9<® 
OF (9. Pare ny one 90 

with *Ce<ete, 


O (% Pe) with t+%<Pas 
(3.29) 


if one face of the wedge 1s illuminated, and 
03 (9, %) +0; (a—9, a—%)+ 
+ en far cong ey en ltrens ig +10 
with Vcgaqr—%, 
oF (9. 9)-+-0; (a— 9, a~9,)-+ 
fermen ith ree ca—s, 
Ha __ | 22 (% %) +07 (2 —9, a—y) + 
fer ee— eo ith aac pce, 
OF (9. @) + of (a—9, a— 9) + 
fen fr e089 = ved withe< p< 2a— n—,, 


0; (% %) +0; (a— 9, a—%) + 
+ et cos (9 7 ete cos Qu— 9 — 9) 


with 2a—1—9,< 9a, (3.30) 


if both of its faces are illuminated. 
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The interrals Vie v5 generally are not expressed in terms of 
well-known functions. However, by usin the method of steepest 
descents, it Is not difficult to obtaln thelr asymptotle expansion 
when kr >> 1. Far from the directions ¢ = 7 + 4 and @ = 2a - 7 - %Q» 
the first term of the asymptotic expansion gives us the cylindrical 
wave diverging from the wedro edie. In the case of wedge excitation 


by a plane wave (1.23), these cylindrical waves are determined by 
the equation 


hy 
E,=—H,= =E,,-fo.* =e oD) 
E,=H,=0, (3.31) 


and with the excitation of the wedge by a plane wave (1.24) they are 
determined by the equation 


(w+ >) 

=F =H,,.g.$  __ 
H,=E,=Ha & ae 

H, = E,=0, (3.32) 


The functions f° and 2° have the form 


f*=— sin 9 
COBy + C03 ye ° 
eae ne 
= COS p + COS ye” (3.33) 


if one face of the wedre cp ca—ez) , and 


esas tN soa I 
cus ¥ cosy, | cos(a— r+ cus (ted ” < 
a sing sin(a —— ; 
s° ~ cos p +cosy,  cos(a— ,)+ cos(a—ya) ° (3.34) 


if both of its faces (a-1<qo<%) are {lluminated. The index "0" for 
the functions r° and a means that the cylindrical waves (3.31) and 
(3.32) are radinted by the uniform part of the surface current cy”), 
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§ 4. The Field Radiated by the Nonuniform 


Part of the Current 


In § 1 and 3 we represented the rigorous and approximate expres- 
sions for the diffraction field by integrals along the same contour 
in the complex variable plane. By subtracting the approximate expres- 
sion from the rigorous expression, we find the field created by the 
nonuniform part of the current. It is determined by integrals of the 
type 


{ 6. 9: Par Cele “ak, (4.01) 


te. 
which, with the replacement of the vartable ¢ by s=//2e * sin + » are 
transformed to the form 
@ 


e'™ j 9 (2. % Ge se” ds (4.02) 


—o 


and may be approximately calculated by the method of steepest 
descents. 


For this purpose, let us expand the function q(s) into a Taylor 
series 


7(2, % Pe S)=Got Us tas*-+... (4.03) 


Let us note that expansion (4.03) does not have meaning only in the 
particular case 


=r? with 9,=0; z, 
p=2a—x— 9, withe=—a—*, 


(4.04) 


when the observation direction (¢) coincides with the direction of 
propagation of the incident wave glancing along one of the wedge faces. 
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Substituting series (4.03) Into Equation (4.02) and then perform- 
ing, a term by term integration, we find the asymptotie expansion for 
the field radiated by the nonuniform part of the current. We limit 
ourselves to the first term of the asymptotic expansion, omitting 
terms of the order (&r)~"* and higher. As a result, the required field 
from the nonuniform part of the current will equal 


o ‘ (» +4) 


pO a et 
E.=—H, = Eul Visa [ 
E,=H,=0 (4.05) 


with wedge excitation by plane wave (1.23), and 


_t +7) 


H,= 2, =Haug' rae 
H,=E,=0 (4,06) 


with wedge excitation by plane wave (1.24). 


By calculating, with the help of Equations (4.05) and (4.06), 
the nonuniform part of the current, it is not difficult to see that 
it is concentrated mainly in the vicinity of the wedge edge. But the 
field created in the region kr >> 1 by this part of the current has 
the form of cylindrical waves, the angular functions of which are 
determined by the relationships 1 


on eee) 1 9. ge 
where in accordance with § 1 and 3 we have 
a 
isso DEN ee SemetS 
g ” (eae ae wat eats) 


= (r8) 





(4.08) 


(1) 
Footnote appears on page 42. 
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and 


f° — sin 
~~ C28 g + COS Ge 


si 2 
g=- TEI (4.09) 


if one face of the wedge is illuminated (that is, when %<m<e—= ), 
and 


eo, ___ 210 sin(a —¢,) 
f GS e608 ty + COS (0 — 4) + C03 (a—y_) ” 
e=— aln? xf sin (a — 9) 
COS P+ COSG,  cos(a— y+ cos(a—y,) ” (4.10) 


if both faces of the wedge are illuminated (that is, when a—m<qo<x ). 
Let us recall that the functions f and g describe the cylindrical 
waves radiated by the total current — that is, the sum of the uniform 
and nonuniform parts, and the functions f° and 2° refer to the 
cylindrical waves radiated only by the uniform part of the current 
(3°). 

Let us note certain properties of the functions ri and gl. The 
function rh = l(a, >, $4) is continuous, whereas the function gt = g 
(a, 9, 5) undergoes a finite discontinuity when 5 z=ae-= mn. The 
reason for this discontinuity is that the uniform part of the current 
differs from zero on the face along which plane wave (1.24) is propa- 
gated (with $) 74-1). In the case of radar, when the direction to 
the observation point coincides with the direction to the source 


1 


(9 = 0) both functions et and gi are continuous. There is no dis- 
continuity of the function gi with $ = $5 = a@- 1, because the current 
element does not radiate in the longitudinal direction. 


On the boundary of the plane waves (that is, when 9=x+q@0 and 
g=2a—n—Go ) the functions f, f° and g, e? become infinite, whereas 
the functions rl and gt remain finite. In accordance with Equations 
(4.07) = (4.10), they take the following values 
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t ® 
s a ilas F 1 
= —___ ___ + -cig9, = ,-ctg—, 
g' cos = — cos ze nes (4.11) 


if o2tn- %o>» and $9 <a- 7, 


+ sio% fae t m } 
P=—— aa te te tacts 
cos — cos ——F 
= sin (e—9,) 
cos(a— 4) eos (e— aed” 
+ 0+ : ‘ 
= +7 te. — 3, cts f+ 


& y—TVe 
Coss — cos a 


sin(a— 9) 
ary as | ESTA) Le (4.12) 


if ¢ =" =- ¢) anda-t < $) < 7, and 


Pl _ ag tl oe 
ees a a 
” = (4.13) 


if g¢=m7 + $4, and $, <a- 7. The value ¢@ = 1 + $) with a- 1 < $, 
< 7 corresponds to the angle inside the wedge, and therefore is not 
of interest. In the direction of the mirror-reflected ray ¢ = 2a - 7 
- $9» the functions ct and gi are determined (with a = 7 < 5 <7 ) 
by the following, equation: 


a a 98 % 
f= = + 
Le ot oT 
' 
+4 ctg(a—e) +508 2, (4.14) 
; fun = + sing 
eS ean 
cos — cos — Se ‘ 


eg . 
+7clg(2—9 —petg = . 


FTD-HC=23-259-71 29 


The functions ri and gi have a finite value everywhere, except 


for the particular values ¢ and %% enumerated in Equation (4.04). The 
graphs of the functions ri and gi (Figures 8 - 13) drawn in polar 
coordinates give a visual representation of the effect of the nonuni- 
form part of the current which is concentrated near the wedge edge. 

In particular, they show that this effect may be usbstantial for the 
fringing field not only in the shadow region (1+%<9Sa,), but also in 
the region of light (0<9<2+q) . The continuous lines in the figures 
correspond to the functions rh crt <0). The dashed and dash-dot lines 
correspond to the functions gt — the dash lines refer to the case 

gh < 0, and the dash-dot lines refer to the case gi > 0. 


cg 


Ge 





-00° dst 


— 
Figure 8. The diagram of the pag fr 


field from the nonuniform part 
of the current excited by a 
plane wave on a he1f~plane. The 
function fl (or g+) corresponds 
to the case when the electric 
(or magnetic) vector of the 
incident wave is parallel to the 
wedge edge. 


Figure 9. The same as Figure 8 
for the case > = 9° 


Let us turn our attention to the next important aspect. As is 
seen from § 1 and 3, the nonuniform part of the current on the wedge 
is described by a contour integral which is generally not expressed 
in terms of well-known functions. But in order to calculate the field 
scattered by some convex, ideally conducting surface with discontinui- 
ties. (edges), the indicated expression still must be integrated over 
the given surface. Obviously, such a path is able to lead only to 
very cumbersome equations. Therefore, henceforth, when calculating 
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iti 


Figure 13. The functions rl and 


gt for a wedge (¢ = o> a = 210°). 
Figure 12. The same as Figure 10 
for the case ¢ = $9: 


the field scattered by composite bodies, we will not integrate the 
explicit expressions for the nonuniform part of the current, but we 
will endeavour to express these integrals directly in terms of the 
functions rt and gi which have been found. 


§ 5. The Oblique Incidence of a Plane Wave on a Wedge 


Above, the diffraction was studied of a plane wave incident on 
a wedge perpendicular to its edge. Now let us investigate the case 
when the plane wave 


Em Eye! ise e + corp + secs (5.01) 


falls on the wedge at an oblique angle 1 (0<1<) to the wedge edge 
(Figure 14). 


From the geometry of the problem, it follows that the diffraction 


field must have that same dependence on the z coordinate as the field 
of the incident wave, that is 
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E=E(x, yet, ° 
ae a (5.02) 


CEES Yy 
ETD 


Using Maxwell's equations 





rot =z —ikE, rotE=ikH, = (5.03) 


Pigure 14 Diffraction by a one is able to obtain the follow- 


wedge with oblique incidence of ing expressions for the radial and 
a plane wave. y is the angle : 
between the normal to the incident azimuthal components of the field: 


wave front and the z axis. 


! 10H, o£, 
be — Panty (+ “ag TOs we 


u 1 02, OH, 
He =a (5 Gp SST 3, )» 


_ 1 OM, cosy OF, 
E,= waar (“or et) 


ane | OE, , cos 70H, 
y= waa (Cae tS): 


(5.04) 
The functions E, and Hy in turn satisfy the wave equations 
bE, +h) E,=0,°4H,+ 8} H,=0, (5.05) 
where 
best ge and A= ksint. rene 


In § 4 we found the fields (4.05) and (4.06) which satisfy the 
equations 


AE, AE, =0, 4H,-+4'H,=0 (5.07) 


and which are created by the nonuniform part of the current excited 
on the wedge by the plane wave 
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E = E,e it seen te + vein gd | (5.08) 


Representing Expression (5.01) in the form 
E am Eel 0 1 shits cone + sin ed (5.09) 


and comparing Equations (5.05) and (5.07), we easily find the field 
created by the nonuniform part of the current with the irradiation of 
the wedge by plane wave (5.01). For this purpose, it is sufficient to 
replace in Equations (4.05) and (4.06) k by k,, and E,, and Hog by 
Ee "* and He" 2 As a result, we obtain 


a (% + +) baa 


V inky t 
i a+) 


H,=E,= Hyg" (9. 90) a aaa | (5.10) 


E,=—H_=E,l* (9, ») 





The angle $5 introduced here is determined by the condition 





elttrcmes ¥COdS) __ g— h(x con gy + y sing.) ' (5.11) 
hence 
__. cos B 
C= iee (5.12) 


The remaining components of the field created by the nonuniform 
part of the current witn the oblique incidence of a plane wave are 
found from relationships (5.04), and when kr >> 1 they equal 


&,=—ctgyf&, H,=—ctgyH, 
E.=—-H, H.=——'_B£,, (5.13) 


¢” siny ? sing 


The equiphase surfaces for these waves have the form 


rsiny -+-2 cos y=const (5.14) 
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and are conical surfaces, the meneratrices of which form the angle 

a/2 + y with the positive direction of the z axis. Thus, with oblique 
irradiation of the wedge by a plane wave, the field created by the 
nonuniform part of the current is a set of conical waves diverging 
from the wedge edge. The normals to the phase surfaces of these waves 
form an angle y with the positive direction cf the z axis and are 
shown in Figure 15. These waves may be represented in a more graphic 
form if one introduces the components (see Figure 15): 


E, = £,cosy — E,siny, 


H,=H,cosy—H,siny. (5.15) 


Then the final expressions for the fringing field in the far zone will 
have the form 


£,=4,=— Ey 


I 
way Ae (5.16) 


eb 
aay 


H,=—E,=— 


Now we are able to proceed to the application of the results 
which have been obtained for the solution of specific diffraction 
problems. The simplest of them is the problem of diffraction by an 
infinitely long strip which has a rigorous solution [23] in the form 
of Mathieu function series. However, in the quasi-optical region 
when the width of the strip is large in comparison with the wavelength, 
these series have a poor convergence and are not suitable for numeri- 
cal calculations. Therefore, the requirement arises for approximation 
equations which are useful in the quasi-optical region. The derivation 
of such equations for a field scattered by a striv will be given in 
the following section. 


§ 6. Diffraction by a Strip 


Let us investigate diffraction by an infinitely thin, ideally 
conducting strip which has a width of 2a and an unlimited length. 
The orientation of the strip in space is shown in Figure 16. 
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Figure 16. Diffraction of a 
plane wave by an infinitely long 
Figure 15. The cone of diffracted strip. The section of the axis 
rays. y (-a < y < a) shows the trans- 
verse cross section of the strip 
with the plane z = 0; @ is the 
angle of incidence. 


Let a plane, electromagnetic wave strike the strip perpendicular 
to the edges. Let the direction of propagation of this wave form an 
angle a(lal< 3) with the plane y = 0. The field of this wave is 


represented in the form 


E= E,e!* cosa} Ge H=Hyel® eo e+ytin ” (6 01) 


The uniform part of the current excited by the plane wave on the 
strip has the components 


fp =0, 
i = 5H os eltv ela * 
thy sin * 


I= go Eocosae (6.02) 


Substituting these values into the equation for the vector potential 
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A= =f r6 yay [3 elt Vere Oar ee Fa 
€ 3 ney, 7 +e 
4 Paes +6— + (6.03) 


and taking into account relationships (3.05) and (1.18), we obtain 
the following expressions in the region r »> ka’: 


i(«*7) 
un [ta (s!n o-— sin g)} © 
"slaa—sing Y2:r * 


i deo 


_ 2 sin (ta (sina — sin ¢)} e 
A, == E,, cosa ar ee Yi-tr * (6.04) 


A, =n, Bobet 








The components of the fringing field in the cylindrical coordinate 
system equal 


E,=—H,=ikA, Hi=E,=ikA,, (6.05) 
where 
A, =A, cos — A, sing. (6.06) 


Substituting Expressions (6.04) here, let us determine the field 
radiated by the uniform part of the current 








ieee) 
Ey —H, = 2 yy coss Mn ne SEE , 
(ara) 
E y= Ha 2H, cose SOE Gin — tinal VEE : e200) 


This field may be represented in the form of cylindrical waves diverg- 
ing from the strip edges 
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Eyes Hy Eye [fechyelteoine steed p 
ort) 
v2. ya 
EM “a Hya-fge(ly eters veh 
ster 


View * J (p.08) 


+ /*(2ye sha (etna. aesi 8 : 


+¢° Qe -(da otro --sin me 


Here the first term: correspond to the waves Crom edge lo (yaa), and 
the sevond terms correspond te the waves from edge eo (Cy = ea). The 
Pune t Lons r ana a are determined Uno the right: habf-plane (iv1<3) 

by the equirt Pons 


PP (Uyer -- f2Q2)-< Mee Gan! 
(pes p92? i ial A 
g(t) g°(2)- *pne—sng’ (e,09) 


Now fetus ff nd the ffeld radtated by the nonunhform part oof the 
current. Assuming the strip ta suff fefent ly wide Cka se 1D), ome fs 
able to approxtmately conagfder Chat the current near Eis upper ede 
Is the same ascon the fdently conduct tng half-plane --esysa, and 
near the lower edpre ft fe the same as oon the liad t-plane sas wea, 
Thevefore, En accordance with § Ay the tletd from the nonunt form port 
Of the eurrent: Plowfoe on the atetp may be reprenemted foo othe form of 


the sam oof the edere eyiindrteadl waves. 


Eye My es Fu lft (iyi ee 
+s" (aye ~ ita (sine: inane ame 

ot My~ Mee let) ea noe 
tha (ni e-= sing) ees ‘ 

-patQye }° om ror 


where the unmet Poms ! and ie are determined Tmo the rteht detltep dane 


-& 
(Ivis 2) by the eqiunt Poms 


SOL eet | Gee om oar 


PO)= FI —79(1), re! 





e'(Y=g()—2"(I). a2) =£(2)— 2°) coe 
in connection with which 
cvs” 2 tis sin zt 
10)*= ~ “noe Sing 
cus ot? +in —" 
P(2)s* = sig ging 
coe ae --- sha —t 
a) ~~ as aing 
— cos 2 4 ain 5! 
6) -——igeaing 
(6.12) 


The functions f° and ge? are described by the relationships (6.09). 


\ 
As a result, the fringing field (the sum of the fields radiated 
by the uniform and nonuniform parts of the current) will equal 


E,=—H, = Eve [P(e eine 


+ f (2) evita (sin esin a ea) 
binkr * 


E, me bla Es Helge (elt ors ae 


rg (2 
+e (eye Fae ies 


Consequently, the resulting field is expressed only In terms of 
the functions f and -¢ which determine the cylindrical wave in the 
rigorous solution (see § 2). The field is the superposition of two 
such waves which diverge from the edges 1 (y = a) and 2 (y = <n). 


Substituting into Equations (6.132) the explicit Expressions: (6.17.°) 
for the functions f and f, we obtain 


FTD=$HC=2 3-259-71 ja 


E,=— H, =, ‘ia [ta Gis sin ¢)] fs 


af! 


in [ta (sta 1 —~ sing)} | @ 
Sacre es Carer ara 
stat | 
ay cos [ka (sin « — sio yi) 
E, == HH, = Hy» ‘a {ha aie 9) + 
2 


Viner 


+i sin (4a (sin 2 — sin 9)} | 


iatee 


(6.14) 


These equations are valid when r >> ka? and lel<z. Moreover, it is 
assumed that ka >> 1, since only under this condition is one able to 
consider the nonuniform part of the current in the vicinity of the 
strip's edge to be approximately the same as on the corresponding 
half-plane. In the case of normal incidence of a plane wave (a = 0), 
Equations (6.14) change into expressions corresponding to the first 
approximation of Schwarzschild [15]. 


From relationships (6.03) and (6.05), it follows that the elec- 
tric field is an even function, and the magnetic field an odd function, 
of the x coordinate measured perpendicular to the plane x = 0 (in 
which the current flows) 


E.(x)=E(—2), Hs (x)=— Hs (—2). (6.15) 


Therefore, on the basis of Equations (6.14) and (6.15) one is 
able to write the expressions for the fringing field in the region 
x < 0 (where ¥ <lpl<e) 


E,==— A, == Fy ie [ta(sino — sing _ | 


sis ae 





cost V 2nkr 
E, = Hye Hy, [onlin sto 
ela 2 2 


(6.16) 
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4 psintiawins slaw ye (v3) 
cos pt } 


(6.16) 


Here one must select the upper sign in front of the braces when 
@ > 0, and one must select the lower sign when @ < QO. 


The resulting Equations (6.14) and (6.16), 1n contrast to 
Equations (6.07), satisfy the reciprocity principle. It 1s not diffi- 
cult to establish this by verifying that Equation (6.14) is not 
changed with the simultaneous replacement of a by $ and of ¢ by a, 
and Equation (6.16) is not changed with the replacement of a by 1m + 6 
and of @ by a- 7 (if execs) and with the replacement of a by 


nm - 9 and of ¢ by n- a (if 3<9<e). 


However, the indicated equations lead to a discontinuity of the 
magnetic vector tangential component H, on the plane x = 0. This is 
connected with the fact that, by considering the nonuniform part of 
the current in the vicinity of the strip's edge to be the same as on 
the corresponding half-plane, we actually assume the presence of 
currents on the entire plane containing the strip. In order to refine 
the resulting expressions, it is necessary to solve the problem of 
secondary diffraction — that is, diffraction of the wave travelling 
from one edge of the strip to its other edge. In other words, it is 
necessary to take into account the diffraction interaction of the 
strip's edges. As we see, it is also necessary to take into account 
the secondary diffraction in the case a= when the Hy component 
of the fringing field must equal zero. 


In Chapter V, we will return to the problem of diffraction by a 
strip, and together with the investigation of the secondary diffrac- 
tion, we will present the results of the numerical calculation based 
on Equations (6.07), (6.14) and (6.16). 
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FOOTNOTE 


Footnote (1) on page 27 The designations used here differ 
Slightly from those used in the 
papers [7 - 11]. The functions f 
and f! there were designated by 
fy and f, respectively. 
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CHAPTER II 


DIFFRACTION BY A DISK 


The problem of diffraction by a disk has a rigorous solution 
[24-26]; however, it is not suitable for numerical calculatiens in 
the quasi-optical region when the dimensions of the disk are large 
in comparison with the wavelength. The physical optics approach 
used in such cases sometimes gives erroneous results. In particular, 
the fringing fleld calculated in this approach does not satisfy the 
reciprocity principle. 


In this Chapter a refinement of the physical optics approach is 
carried out. First the diffraction of a plane electromagnetic wave 
by a disk with normal incidence (§ 7-9) is investigated, and then 
(§ 10-12) diffraction by a disk with oblique incidence of a plane 
electromagnetic wave is investigated. 


Normal Irradiation 


§ /. The Physical Optics Approach 


Let an ideally conducting, infinitely thin disk of radius a 
(Figure 17) be irradiated by plane wave 
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E,=—H,=—H,,e", 
(7.01) 


£,=H,=0. 


The uniform part of the current 
excited on the disk by wave (7.01) 
is determined by Equation (3.01) 
and has the components 





p= —ae Mes 2 =2=0. (7.02) 


Figure 17. Diffraction by a Let us find the field created by 
disk of a plane wave propagated 
along the z axis. this current. 


Since the diffraction fieldin the far zone (R >> ka“) is of 
interest to us, the vector potential 


e ae 
A(x, y. z= (pap ( 349, 2 a 
ve rf j ahaa (7.03) 
may be simplified by using the relationship 
r= // K+ 9" — 2pR cos. = R —pcosQ, (7.04) 
where 2 is the angle between p and R, and 
cos OQ = sin § cos ('» — 9). (7.05) 
As a result, we obtain the simpler equation 
_ @ 22 
1 ft? 7 ik, con @ 
A(x y. =a cel Ne yy even? dy, (7.06) 
0 
Continuing by using the equations 
H=rot A, rotH = —ikE, (7.07) 


it is easy to show that in the spherical coordinate system the 
fringing field components with R >> ka? equal 
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E,=H,=ikA,, 
E,=— H,=ikA,, 
E, =H,=0, 


(7.08) 


where 


A, = A, cos 9 — A, sln 9, 
: } (7.09) 


A, = (A, cos ¢ +- A, sin p) cos ®—A, sing, 


Substituting here the values 


Ay Hex zit, A(hasingy SS, (7.10) 


A, >= A,=0, 


which result from Equations (7.02) and (7.06), let us find the field 
radiated by the uniform part of the current in the form 


E,=H4,= ~ iaHys 5 yr cos 0/,(ka sin dS, | 
(7.11) 





aie eee 


rs, (tasin’) 2. | 


The function Jy (ka sin #) is a first order Bessel function. By 
using its asymptotic expression 


ho sini 


J(ka sin) =|/ 2 cos (ta sind — *), (7.12) 


which 1s applicable when ka sin §>> 1, one 1s able to rewrite Fqua- 
tions (7.11) in the following form: 


; a sin 
Ey =— iM V apy OOO TEX 


> cae Coe. 


Emcee a ee at (7.13) 
Sew iby = nee | 


$ Ink sind a x 
i ae 
—l [ea vl 8 —-— i taaind— 2) 
fr ce | 
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The resulting equations show that in the region Ra’, kasind >! 
the fringing field may be investigated as the sum of the spherical 
waves from two "luminous" points on the rim of the disk, the polar 
angles of which respectively equal y »= 6 and » = 7 + 4. It is not 
difficult to see that these waves satisfy the Fermi principle. 
Actually, of all the points on the disk's surface, the point pe # a, 
¥ = @ is the closest to the observation point (R, 8 9) , and the point 
eo @a, » = 1 + ¢ is the furthest from it. 


However, Equations (7.13) describe the radiation not only from 
the two "luminous" points, but they determine the field radiated by 
the entire "luminous" region which is adjacent to the line connecting 
the points p 7 a, » = ¢ andp =a, punt 6. 


Let us show that the luminous region actually makes the main 
contribution to the fringing field. For this purpose, let us cal- 
culate the field radiated by the currents which flow inside the 
sector encompassing the line » = @ (Figure 18). Let us take the 
angular dimensions of the sector in such a way that its arc, which 
equals 2aey, would occupy the first Fresnel zone. When this is done, 
the angle % will satisfy the equation 


a(! cos @) sind ae, (7.14) 


In the case being investiga- 
ted by us, when the condition ka 
sin §>! is fulfilled, we have 
from Equation (7.14) 





co eo 1 tet (7.15) 
hence 
Figure 18. Calculation of the 
field radiated by the "luminous" = (7.16) 
region of the disk. a= V poy. : 
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The vector potential of the currents flowing In the Ind!eated 
sevtcr ts determined by the equation 


i} oft 7 4 
Ay --% He Qe dy funn O cos +e dp, 
oe 6 


A, = A,=O 


(7.17) 


Taking into account the condition ka sin ®@>1t, one may show that the 
field created by the currents of this sector will equal 


Eo Hi~He spaay 089 X 
sia gel?” —thasin 9 HE. ry 
ane Rr? +0(Fa5)- 
r) 18 
Eo —My~ Huy 8 x (7.18) 


cose ef#® = ~itasin Ode I 
xis +0( saa): 








The amplitude of the expressions which have been found is approximate- 
ly ¥2 times larger than the amplitude of the first terms in Fquation 
(7.13). Moreover, expressions (7.18) and the corresponding terms 

in Equation (7.13) differ slightly in their phases: the first have 


the factor e® , and the latter — the factor—e¢*, The result obtained 
is similar to the well-known thesis in optics that the effect of a 
wave 1s equal to the effect cf half of the first Fresnel zone (see, 
for example [27], p. 132). 


In the vicinity of the directions 0-0 and 0=2, when the 
azimuthal components lose their meaning, for the purpsse of studyvine 
the fringing field it is more convenient to use the Cartesian 
components . 

E, = (EF, cos§-+-£, sind) cos ¢ — EY sing, 
E, = (E,cos§ + E, sind) sine 4- E, cos9. (7.19) 


Turning to Equations (7.11), we find that when §=—OQand d=» 


E,=0, E,= iH. “FS. (7.20) 
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Conseyuently, in the physical optics approach the field scattered 
in the directions 0=0 and %-=" preserves the polarizaticn of the 
incident wave. 


§ 8. The Field From the Uniform 
Part of the Current 


Let us proceed to calculation of the field created by the non- 
uniform part of the current with normal irradiation of the disk. 
Since the latter is concentrated mainly in the vicinity of the disk's 
edge (p = a), the vector potential corresponding to it will equal, 
tn accordance with Equation (7.06), 


aon F 
Ane. 5 arf Fp, qyerenbem—nag, (8.01) 


The inner integral is calculated with ka sin®>I1 based on the 
stationary phase method (see, for example [21], p. 256), and Equation 
(8.01) is transformed to the form 


PD ee eee 
Aix 2=2+V aa ee x 


i ; e (8.02) 
xf 3 (?, ,) aiid) (p. »,) elt? ea] : 


which allows one to interpret the fringing field as the field from 
a luminous line on the disk. This line is a diameter, the polar 
angle y of the points on which equals 


=P and $,=2-+9. (8.03) 


Assuming the diameter of the disk 1s sufficiently large in com- 
parison with the wavelength (ka >> 1), one may approximately assume 
that the nonuniform part of the current near the disk's edre will be 
the same as on the corresponding half-plane (Figure 19). On the basis 
of § 4, the field from the nonuniform part of the current flowing 
on the half-plane ~°%y¥,<a@ may be representedinthe form 


FTD-HC-23-259-71 48 


i nase ’ 
E, (=A, (= Es, PY ae oe, 
H,, (y= ik, (1) cos ® = 


ifaaee 


(8.04) 


=H,,- ‘ e ° 


g 0) ee 


and similarly the field from the current flowing on the half-plane 
—a<Sy,<O, may be represented in the form 


‘ ahee 
Eien. (2) = Ey, 1) ae =e". 
: arse (8.05) 
H, 2) — thA,, eee Mon g ens Ear oS giaasind 
Here 
ae (sete) 2 
A(t) yf Fre Lrretmtas 
(8.06) 


A(2) = + Ve (ere) hy (9) elt 1806 an, 


1 


and the functions rt and g” are determined for the right half-space 


(o<0<+) by the equations 


com + sla 
rM=fN+54,. as 
Wg (8.07) 
: cos -3>.—‘sla 
ra=1\-c, Ff ia Aas 
eos-g-— slag 
a (1) = g(1) + oeh 10—-— a 
, (8.08) 
coos +} tle 
e'(2) =e), g2)—= — | 


From relationships (8.04) ~ (6.06), 1t follows that 
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f ii, (1) em tay sin® d= 


8 
-o 
=—- é,,;1' (2) e/*2 and 


and 





t ii, (1) eit y, ain *da= 


=a BEM (emene , | 
7 "" 
Figure 19. Diffraction by a disk. | 


The half-plane L lies in the a“ # “gi (Ie ~ite mn 
plane of the disk. The edge of Bak cosd °° % (8.10) 
the half-plane is tangent to the © reer : 
circumference of the disk at the fi,a@e om aa 
point ¥Y; 2 a, xX, 2 0 (a is the ~c 

le He. g*(2 ¢ itonine 
radius of the disk). Tikcost /e,,°8' (2) 


In accordance with the assumption of equal currents on the disk and 
on the half-plane, one may consider the following equalities to be 
valid: 
fro. $,) ete ein® 5 oo i J (») cinta, | 
— 
fro. get"? dpm j F (el 8 dn, (8.11) 
—* 


Therefore, the fleld from the nonuniform part of the current flowing 
on the disk will equal 





la£, ‘ ta cin 9 
a a ries [r are a 
~i [rasine~™ 
peer Dyae, 
lal 1 {2a tn 02 (8.12) 
ya, pats [at ae | ) 


sae elt 


where in view of (7.01) 


E,, = —H,,cos9, H, =— Messing. (8.13) 
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For the direction§=0, we have according to Equation (8.01) 


rs} ae 
smb [are #dp, (8.14) 


but in accordance with equalities (8.09) - (8.11) 


f1i,(0» #) p= 35 Hea-cosy, 
e 


A (8.15) 
fit. $) dos 75 Hy sing. | 
‘ 
Consequently, 
te 
A,= <. he Jeoseas ff (p. Ydp+- 
+f vy fr. 1,0 9) a]. 
ft? (8.16) 
Ame -fensafe 0. Nar 


+ fom vv, (. na] = 


that is, in the direction of the main fringe (0=0) the field from the 
nonuniform part of the current equals zero. 


By using the Bessel functions Jy and Jo for the field from the 
nonuniform part of the current, one may write the equations 


B= — Hy =F {IN ()—P (DA (Ra sind) +- 
FHP +P (J, (easin 9) 
Ey=H, = 5 (1a! (2) — 2" (1)] 4, (hasin 8) + 


“file! (2) +9" (IJ, (ea sin SS , 


(8517) 





which with ka sin ®>1 change to Expressions (8.12), which were already 
found. In the direction0=0, these equations give a field which 
equals, in accordance with (8.16), zero, and with intermediate value.: 
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they are dnterpetated. Stnee the frame tt bem from Bauat bers (s.r 
bo Ramat fons C8. bo) tis met comteteds untque, ft the ateutar foteae ct 
DOs F Banat bene Ga. Ee) may eive a eettaln error,  Uhi: errr Le 
eto very ctpnteteant, stnee dn this taterval the rledd Crem the ant 


form port oof the current fa barye. 


WO, Phe Toul Fre td Bete feat tered by a 


Mek with Nerreul Trrid fat ton 


Tarntie: tee Kqnat tens Ga.ary, Cede and GALETY, Pet us tepresent 
the Cheda Prem Che oenmunf Corm paret ef the enrrent fi the ted bewlne 
Porm 
Boe My: stattes S2t8 J, (ha sin) "e oa - 

alle (2) -- LOD) (aa sin 8) . 


BAU) PEO (ha sin O)} cos 9 ® e ’ 
an 


Ey Hy. iatlys “in Fcus ! (ka sin) 9 a ee 
— tates tty (2) — a (NY) J (ha sin 8) + 


ila (2) Fe 4 ke sin &)} sing ee . 


Here the (hpet terms, ihe ts oreadt ty apparerit . pep renent three (hehe 
Crom theount form pout oat the cirrent tahenm with the oerpentte loti 
Aso a result, the fetal “Medd seattered by the dtek cutie Vee ase 
cam et the flelds tadbated ty Che oute! Sorin ind monn fee gaps ah bite 
eurrent . with be exprencsedh only Ino teprms of the Cunet Gere i ond 
woateh addetermine loo the ebeeton: sobuttotm the ey birdied waves eres 


the TabPeplane a edpe 


BE, Hy cose He) 
~- fh sin 8) i (. |- 
LE) 4, (4a si Oy) © Ro 
Feoeth, os ae sing (KC) > yt} 


— p(y (Air sins 0) \ vile a) | 
wees (Ra sin O° a 


Subst ftut fap here the explfeit expresstons Cor the funet fone 


Poand oy we arelve at the fini) express fons for the Crinrtny tleld 


jee f Jy (ha ala 
E,= -Hy=- aa | ea oat) 
oan 
~j: J, (ta o. Lae: 
con ‘r 
glee J, (ta st (9.03) 
=Hi,= : ml ¢ oly sof) S 
ey 
aes J,(kas aah) ‘|e ‘ne. 
cos =" ; 


These equations are valld in the right half-space (0<0< +) - In the 
left halfespace (4 <?<*) » the fringing rleld ts easily found by 
assuming that tts electric field ts an even function, and Its mayrnet te 
fleld an odd funetfon of the so coordinate: 


E,(2)=E,(— 2. 


H(z) =—H, (—2). (Q,0n) 
Consequently, In the region 2 < 0 (that fs, when She) 
lays | Iy(ta stad) 
E,=—Hy=-— "4" Iy(a vt — 
cot ~s° 
if ka fo 0 
ra ( oe reo 
sla H 
ial, | J,(ta sta (aan) 
Ey H, = lattes, | Jufte sine 4. 
con 
tt 
+i Ais sia vsing 
sian : 


Assumtne that tn Equations (9.03) and (9,05)@==0 and Qa, respectively, 


we obtain 


ika® Vlad 
Ey == —"p Hea "ge E,=9, (Ose) 


whieh Is equivalent to the physteal opt les appreach [seer Faquat bon 
(7.20) ). 


dO | a a La HS 


Expressions (9.03) and (9.05) agree with the result obtained by 
Braunbek [29] for the scalar fringing field in the far zone. It ts 
also interesting to compare these expressions with the precise numeri- 
cal results obtained by Belkina [34] by the separation of variables 
method in the spheroidal coordinate system. It turns out that even 
with ka = 5 a satisfactory agreement is observed between our approxi- 
mation method and the rigorous theory. In Figures 20 and 21, graphs 
of the functions V"(0) and Vv" (8), are presented which allow one to 
calculate the fringing field on the basis of the equations 

Ey = — Hy = Buy (8) cos 9, ] 


a . glt 
Ey=H, =F BV 0) sin 9. f 


The continuous curve corresponds to the rigorous theory [34]. The 
dash-dot curve corresponds to the field from the uniform part of the 
current, and the dashed curve corresponds to the field calculated 
according to Equation (9.03) and (9.05). 


(9.07) 


Oblique Irradiation 


§ 10. The Physical Optics Approach 
Let us investigate the general case when the plane wave 


E=E, eft (wate ye cos a (10.01) 


falls on the disk at an arbitrary angle to its xis. Let us take 
the spherical coordinate system in such a way that the normal to the 
incident wave front, n, would lie in the half-plane => and form 
an angle y (<1<-+) with the z axis (Figure 22). Adhering to the 
investigation procedure used in the previous sections, let us first 
calculate the fringing field in the physical optics approach. 


The uniform part of the current excited on the disk by wave 
(10.01) is determined by Equation (3.01) and has the components 


—¢ thysing 40 ¢ 
fg Mae, Pah 


so Meet "T, j=0. (1000) 
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Won nang 
“ CETTE PP TT) tP iri eq 
CCC 
SANS ERR Ae 
NT TINE LET ETL Wi LA 
ACCC EEE 
UNM VERSE 
TRL Tp TT Td 
BLISS ee 
TTIW TITEL TTT Wt oy 
PAN TTT a a ea 

PCC ee 
PL ALERT TE TT Tat Rt 


















Figure 20. The function 4 for a disk 
with a normal incidence of the wave. The 
varicus curves correspond to different 
approximations. 


The field radiated by it 1s found, as was done in § 7, by Interrating 
(with the condition R >> ka*). In the case of E-polarization of the 
incident wave (E.lyoz), this field equals 
=H = ; J, (tap) e* 
fe =H =iaF,, COS 1 COR CONG a 


See T (10.03) 
Pah a= + I: (ta Viet pey aft ® : 
E,= y= (aE y_-cos y sing ‘ 
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tae CA fa-$ 
Yt tt rT Reg’ Wet? ["]- 
NI} 7 | | (ATT TE RT | TT fe 
NSS Ree ARR Re 
NT PPT TEE ea ATT TT OV 
BIBS ae eaNIeee 















SUG) RRR 
2ST Se 
Coc 
CCE Co 
CACC CCC 
ROC /AR SR 
POR 
PCA a 

















BhCEE 


SCORERS 
SeenkeeN/-s8t Aneee 













Figure 21. The function ¥"() for a disk 
with normal incidence of a plane wave. 

The various curves correspond to different 
approximations. 


and in the case of H-polarization (H.1yoz) 


ikR 
Ey = H,=—iaH,,cos 8 sing 1G2V BED =, | 





Yee 

= -_ 10.04) 

= a Fytka) Vi Ea it® 
E, Hy = —ialt,, cos p Vitae 


The quantities a and uw in Equations (10.03) and (10.04) are deter- 
mined in the following way: 


posin sing —siny, (19.95) 


4==sin 9 cos 9, | 
VY Pp >0. 
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Assuming 92-5 and ae—1(3<0<*),1n the resulting expressions, 
let us find the field scattered by the disk in the direction toward 
its source. With E-polarization of the incident wave, it equals 


£,=-—H,=- fafee tet 4, (2hasin 8)e 


£,=H,=0. (10.06) 
and with H-polarization 
E,=H, = leet" ka sinh, } (10.07) 


E,=H,=0. 


Using the asymptotic expressions for the Bessel functions, one 
is able to show that when R >> ka” andke//J*’Pp'>1 the fringing field 
is radiated from a luminous region on the disk. In the case when 
Yi +p'+0, the luminous region is increased and in the limit (when 
A =» = 0) the entire surface of the disk starts to "shine". 


§ 11. The Field Radiated by the 
Nonuniform Part of the Current 


Let us calculate the field in the nonuniform part of the current. 
J*(p. G=J(p, gpelte sintee?, (11.01) 
Its corresponding vector potential 


=.2 —ihpV Ftp pteon(t—) ,, 
=e ferfre we ibpV 7 *toteon(; Muy (11.02) 


_by means as the stationary phase method is transformed with 


« a 
Sor = nee ae La. ae 
- YY wieee Vat par [Jee we de 
' 


‘aes : (11.93) 
- ifs. 4) nea 4h 
» 


Here 
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w=, Y= 2-bs (11.04) 


are the stationary phase points 
and the quantity 6 is determined 
by the equalities 


. i 





Figure 22. The oblique inci- From Equation (11.03) it 


dence of a plane wave on a disk. 2 
n is the normal to the incident follows that with R >> ka" and 
wave front. ka fps 1 tne main contribution 


to the fringing field is given 
by the luminous region adjacent to the line y = Vy> ¥ * Vo: Thus, 
the stationary phase points vy» V5 physically correspond to the 
luminous line on the disk surface. 


In order to calculate the vector potential (11.03), it is 
necessary for us to first express the nonuniform part of the current 
on the half-plane in terms of its field in the far zone. For this 
purpose, let us introduce the auxiliary coordinate systems Xp» Vy 
and X55 Yo (see Figure 23), and let us take the following designations: 


ais B, (ay, Bo) are the angles between the normal to the incident 
wave front and the coordinate axes Xy> Vy (x5, Yo)3 


ated * -99) is the angle between the z axis and the projection 
of the indicated normal on the plane xX, = QO; 


(6) = -$5) is the angle between the z axis and the direction 
from the coordinate ortgin to the point P(y,> z) which lies in the 


plane x, = 0 and is the projection of the observation point P(x, v, 2); 


mr is the distance from the origin to the point P(Yy> a) 


The quantities introduced here are determined by the equations: 
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c0S a, 5 —-sin cos },, cos sf, ==sin ysin4,, 
ae 3,=e—~§,, 
sing) = 5 ahs » cos, Sel 


% 
we. , Mew Oy — oe) 
sin %, = Vieni tsint(—9) ” (11.06) 
ieee 


=A, i=—F , 


n=RV iano sin 9) - 


Furthermore, let us write the expressions for the field from the 
nonuniform part of the current excited by wave (10.01) on an ideally 
conducting half-plane —~co<y,<a. In accordance with § 5, they have 
the form 


oltre) OW, 
EB, =aelntmne I @. my bere aE oon ant (11. 07) 


ls ay +) iha(sta @? stag) 


H, mn eltemayy or (Ps %)° V aby, 7 








where 


ky ==hsina,, 11.08 
k, (sia 9, — sing! =k YP Pp, } : 


ate tt Pept ud 
2 wser 





(%» 9 — 
Pw % sin 9} — sin y, sing? —sing, * 


‘ 

h-9F nt | 
2 

| 

} 





—sin —3— +08 
2 
a' (#1. of = ———__,_______— -. a 


sin p? — sing, siny?— sing,’ 


(11.09) 


. Je 
(—4<9.<3). 
On the other hand, this field may be expressed in terms of 


the vector potential 


it FP Pe 
Ving vtG—y * 





a @ 
Aaa [senate ray | eens (11.10) 
-—2 


-~e 
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By means of equation 








~ pero 
1 ip YREG 
Hy) @D)= 75 (5 

—- 


q==Vp—F, Img>0, D>O, 


which follows from Equation (3.10), if one substitutes z= t, w= C, 
ad = -ip, k = -1D, in it we find that 


a 
Im the, coe ‘st Llp ~ u\st pitnecs 
A= jel a OL PY — yer dy. (11.12) 


Taking the fact into account that 
fizy.ay the nonuniform part of the current 
is concentrated mainly in the 
vicinity of the half-plane edge 
and using the asymptotic repre- 
sentation of tne Hankel function, 
we obtain 





a ik(z, cos a-r, sin ett 


© @ 
’ 2 T thy, 
A=— hare e fre ‘dy, (12235 
-@ 


Figure 23. Diffraction by a disk where ‘, = cos §, -— sin a, sing, = 
with oblique incidence of a plane 

wave. The half-plane L lies in =sin {sin}, — sin 9 cos(p, — P) 

the plane of the disk. Its edge Seige ee 

is tangent to the circumference V ees (11.14) 
of the disk at the point x, = 0, Penn's sia’ (yy = 9) 


y, =a (a is the radius of the 
disk). In the case when v= é [see 
Equations (11.04) and (11.05)] 
the function 4 takes the value 
0, = —fPpp (11.15) 


Starting from expression (11.13), it is not difficult to shew 
that the fringing field in the far zone is described by the following 
equations: 
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E, ==— thsin a, cosa, sin9,A,, + ik sin"a,A,, (11.16) 
H,= — ik sina, cos pA, , 
where 
A,, =, (1) A, =O, (h), 
ne seve (11.17) 
gat ze 1 COs B+r, maa 
é 
1, (Qt J J,,(0 er Pay 
— (11.18) 


f é (%,) = j 4, (3) eames 


Then by equating expressions (11.07) and (11.16), we find the desired 
connection between the nonuniform part of the current on the half- 
plane ~®©<y,“2@ and its field in the far zone 


Ly) — nia ae tage (Ou 9h) oe FR 


jee Eq. f* » % = 


(11.19) 
— cosa, tg eH B" (Pu 9° ener, 
One may show in a completely similar way that the nonuniform 
part of the current excited by wave (10.01) on the half-plane 
—oo<y,Ga creates, in the far zone, the fringing field 


E,=—ik oe 2, cosa, sin ?,A, + fksin"a,A, , (11.20) 
H,, =='—iksin 2, cos 9, A, , 
where 

A,,= 8,1 (ads A,, = 0,/,, (Y2), { 
1 Te itr cin, =e, coe pti > } (11.21) 

= V ine my 

, 
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Py ee 
L(t) = f J,,Capet? Pun, | 
--@ 








“ VT (11.22) 
LWad= J cet dy, 
-” 
On the other hand, in accordance with § 5 this field equals 
~ikz, cos 5) Py nee 7) shaolin s)—oln Gy) 
E, =e fy ” Bord! (1%) Sra e ‘ 
. +) (11.23) 
eT) satin ul 
ihre », e saGin ¢,, wn Ge) 
He ike, 9% HBP Fs ) : Vick, : 5 
Here 
k, (sin 9, ~sin #8 ) = — VR Pye, (11.24) 


and the functions /'(9,, 5) and g'(%,93) are determined by the equations: 








tate coe EM ; 
a 2 cus ¢, 
4 0 a es 
Pu %)° sing) — sin % +n go —sing,” 
W) 9 
a-th 7) + A 2) 
8 (Pe 9% cee ae Ed aaa Se COs fy faby25) 
ae sig — sin gy sin 9? — sing, 
‘3a a 
\ 7 <9, <3) . 
Equating the quantities (11.20) and (11.23), we find 
inaV is 
1, (= = war nea Hg! (Pa, % )s 
' ce eltaVingw : ” 6 
1, (0) = jaz sinta; — !Eng," (Pes 9. )— (11.26) 


--- cosa, te 9,/1y, 8" (Par #2 ))- 


In this way we established the relationship between the nonuniform 
part of the current on the half-plane and its field in the far zone. 
Now let us return to a calculation of the field from the nonuniform 
part of the current flowing on the disk. 
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Since the disk is assumed to be large in comparison with the 
wavelength, the nonuniform part of the current in the vicinity of its 
edge may be approximately considered the same as on a corresponding 
half-plane. Consequently, the integrals in Equation (11.03) will 
approximately equal the corresponding integrals from the current on 
the half-plane: 


i J,.(0. $4) e™ VET Py wn 1, (ys 
J+. (, yer PF Map == 7 (>), 
(11.27) 


j Jl ga) Pde =I, (od) 


fn (0, g) elt” F Papa t (4,). | 


As a result, the vector components of (11.03) may be represented in 
the following form: 


1 z J 
4,27 VS oT, G+, Oak | 
A= TV ao ee WEA, Cdl | 


Then substituting these values into the equations 


(11.28) 


E, =ikA, =ik|A, sia (9, — 9) — A,, cos(, —9)], } (11.29): 


Ey =ikA, =ik[A, cos (y, — 9) + A, sin (4,— 9)] cos 8, 


we find the field from the nonuniform part of the current flowing 
on the disk 


i= 
‘ tar 
= pera | Sake, ase, ae sin(; = 9) __ 
E, = Hy= ¥ /takaVia R | av, aineceaee 
to 
— Fata 608 tee |X 
X la! (Ps 9) penta FH gt (p,, 9) eit en (11. 30) 
See e ——— 
~&,,, Sar eee (Pir 2 Jemite FT 


te ee 
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s 
iy 

__ae_— cos d eft® 
E£,=4,= ects 


-~ cus Oi — 9) 
V2rka | Vita +ne R | H.,,[ | SiN 2, COB 9, + 
aa! sin (+, — 9) 


sin ay te 9 Jie" (71. 7) peel ee 


—ig' (5. 9, ell’ Wea] 4g - sin(4, 


- (11.31) 
“Sinta, 2 “tints, 
x<if (?,. ?, Sen iaaV ite ee os if* (Par ?, ) elite ive ep ] } 





The resulting expressions are valid when sa/i*f_">1 
may be slightly simplified to 





—_—_ 





They 
ute eine 
E =—H —— 
: * VistaYh tae RX 
X{— Hay cos tts) tat (Ps *) enteaV ie 
— ig! (am 5) eteYPar] — E,, SC — x 
IF (,, of emt PF. if (9, 9h etary . 
; i2 (11.32) 
¢ ita 
Bia eee a 
=e } Qsta Viper X 
x{- was EM Ie" (Ps. %; * ) entero 
—la'en eye | + E,, cosd 2G — 0x 
x (f (71. ?, sents VTA if (Pa 7) efte? giz ] } 
if we use the identities 
sta (¥,—9)__cos(t — 9) cos # sin (2,—9) 
ain a, c039, Tinta, COS% te9,= n'a, 
cos (4 —¢) 1 $10 (4-9) 0. == £08(%)—9) (11.33) 
sin sina, C03 v1 a aaa cos @ tp = “aint a; couD ‘ 


The operations carried out above may be briefly summarized in 
the following way. 


The field from the nonuniform part of the current 
on the disk 


20 a 
IaR ; 
Axe Se {ay[sie, gel ae 
) 0) 
is found (without direct calculation oa 


1? the current) in terms of the 
known field of an auxiliary haif-plane 
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Je{s,c18 G29, ole aes 2 
Ant “freemen 
-——@ 
by a replacement of 


$30, 9 cap vy Jseredy 


in those cases when ¢ # 1° The functions ¢ and 1 are determined 
by the equations 


@=sin ysin y — sin B cos (¢— 9), 
&, = sin y sia y — sin 8 cos (py Tasers ae | (11.34) 


Solution (11.32) was determined exactly in this way with kay’i'+y"> 1, 
when, for auxiliary half-plane whose edge touches the rim of the 
disk at the pcints »y = 6, » = 7 + &, the phase 4 was equal to 9. 


A solution to the problem using this method also is possible 


in the case 
\ 


%=Le= yy (11.35) 
when ies 0. The direction §=y7, => corresponds to the princi- 


pal maximum of the scattering diagram, and therefore is of special 
interest. Substituting the relationships 





sing? —1 





| 
(11.36) 
} 


¢ u 0 | 
= hic sin? a, (E,,, -L Cus a, tg 7, Ay,) cos ° 


which follow in this case from (11.19) into the equations 
te 


A,= = _ ((1,,c0s y-+/, sin 4) dy, 

: (11.37) 
ith * a aad 
A,== 7) (7, sing —J, cos 9) dy, 
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we Thad fhe held radbated by the nota bpori: pore oe thie current fa 
the direst fon er the gauioetpadl maxtmaem, With the @epoturtsat ben of 


the Tueldent waive (Ey; vec), Pte quarts 


eR 
Fy. H, 0, 


aed with othe tHepotaetsar ton (if, Laosypot the taeddent waive. Pt equerds 


Foo Hie Rigs oe 


“18 nm sia®y 


X (2K cosy Shetty) 
FL H, °0, } 


where 


r 2 
dd " 
“ f vi ~sint yz sin? | pe in sin? y sin? > dy cir. 
8 0 


wihetoree e@emetete ae bitppr te tategiate, Mr fle oresudbt trae anpresoheme, 


Pe Voldowe that with the mata dem ero the Paetdent wave modlaetostt Pom by 


VO the phase of the fheld Crem thea i fea pearb ooet tie eurrent: 


Psochaneed by bSsehy cms dtocwn tu tte ease om on hade@ephinme. Fr oy ors 
them tine di eeercnce feetweet fhe pebartoat Port iisoappss ara, ci 


Vimit qwhen yoo OY wip dve at fae pemathorss: aasdt qe. ded. 


. : , ' : at) 
SEE EOE pe a Bel APE Ce ee Pn ay us pore ie BRN erRCY Eo G obec 
ad Chia wed, hayes waned © Tyrant Moves qr oy ee er ghey Get thay 
{Aha SEBO lich SATE t poteh AM op de® aecpe at 8 SES SCE leah 
ae $21 tline et i . tlyes ve genes Y Paes a tho tase eee r 
Lae sai aa nm fost : \ 
pile my y t 1 
ta 4 se ESHA Tea a ; . ' . 
‘ ‘ q vege. ft i te . ret 
: tae ees . ¢ 
. . 


detall the fringing fleld in the ineldent plane (=0, ga 7) where 
the expressions (11.32) take the form 


E,=a-H= 7 re —f'(ye fipaye™™)Se 


Y 2stap 
withe=—, &>y, 


ae . ae 1 itap —thep tk 
EN aol (ye +P Mer rg (12.01) 


with t=F- o<y. 


a7 ion 
7 3 hen} CU 
=-H,= y raenaiie* ff (2) e"* p(y Pa Ra 
with 9-5: 
ae ike shay 
E,= += 1 Tada | g(lye T+ ig! (2)e i 


withe=+, O>1, 


i- 
xi jag 
_- Hor thay ’ item @ 
E,=H,= ea ar(2je dig (ey S— eacias 
with p=t, <1, a 


4 
E,=H,= Feaapile (2) ¢* — igh (el 


s 
with g=— 7: 


The funet dens rr) and ray) correspond toe the fleld of the auxiliary 
half-plane —ec@vGa , and the Cfunetfons rh(s) and rho) cerrecpeonad 
to the fleld of the half-plane m@S"Soo, Tn aecerdance with 


squat homes CLE. y and (11.29), Chey ore deterrftned ty the oxpresc tens: 





' cosy an 7 eos : f 
PW=10— jaaine + A= Sa ae 


§+1 
iin T4031 
PVs aes, gia : 


siny —sind* sing—sme* Choral 


(Matert Dem voentinuged: ono aext perce.) 


sin alee | —tos Sty 


2 2 
'Q=(O+a5—ae 12=—asmam 
oy O+y (12.03) 
® sin z + cos —3—- 
8 (2)=5 8 (2+ 5 82) — -— aa 
if =F and ®<4 ; and 
e+ —y; 
ee sin? TT ogg "51 
PdQ=f— way + 3in0" (w= ae ie 
4 ae sin — +1 cos = 
ea siay Fane: a()= "sny+siov' 
(12.04) 
re ie a a 
1 cosy 2 2 2 
f (2) = f (2)-+-—- sing + 3nd’ it )= sy siay4-sind = 
ino $7 cos * Ft 
sin -~cOs 


g'(2) = g(2)4- amar » #Q2)= sing+sne° 


if e=—-F and <+ : 


It was mentioned above that when ka >> 1 in the direction 
8=y (0° 7<55°), P= the field from the nonuniform part of the current 
is negligibly small in comparison with the field from the uniform 
part. Therefore, for the field from the nonuniform part of the 
current one may write, with the help of Bessel functions, the 
following interpolation formulas: with ee 


Ve Ay 8 (HL — 


iE MNOS 


Cite) 
Ey= Hy = "3" {9 (— a! CAO — 
— ila (1) +g QIAO» 
and with g=—> 
E,=— Ay = FQ AQ+ Cen 


F ; ela 
tHP OD +P A: OF R > 
(Equation cont fiaed on tes 


ae of . ” . . cae | 
BPP Ce 3a PG AT] OS 


£y =H, = he ttg'()— a NCH 


ff (12.06) 
bile (V+ e214, 0) 5» 
where 
&= ka (sin 6 — sin y), 
t= ka (sin 6 +- sin 7). (12.07) 


These expressions are valid in the region 0<@<—; when 9>1 andi >! 
they change to Equations (12.01) and (12.02), and in the direction 
O=y, e=+ they give a field equal to zero. 

Using specific expressions for the functions rl and ge, it is 
not difficult to establish that the total field scattered by the 
disk in view of Equations (10.03) and (10.04) may be represented in 
the following form: 


with P= 
E,=-- Hy =Se{t (12 40— 


eftt 
— HN +124, Olan. 


ten (12.08) 
£,2H4,=— 3 {le—eQV,.0— 
e 
— Hg (1) end} » 
and with p=—~——> 
E,=— Hy = SP) — 1014, + 
+ iI) FAS ON 
A 
(12.09) 
Ey=H = "Fle ()--2 M404 | 
+ileM+eQOr ge. | 
It is convenient to write these expressions as follows 
iof,. — IR 
E,=-- Hy= “3 MA, Dp 
(12.39) 


* sar 
E,:aH = teva, yt 
r . » ( {) R . | 


PT eH CH? 3-259-71 69 


where the functions £ and £ are determined in the region istsy by 
the equations: 














50, 1) |= — GL sj - at) with p= >, 
ra, 1) ain Ft cos- yt 
58, 1) JA 2-7; «with P=—-—=, (12.11) 
58, Y) sia cosy 
' 
and in the region > <h<e 
3, 1) == 16) — iO) with e=F. 
x8, ?) cosy a 
sa, oe + AG) A =with @=—s. pares 
wo, y) * eoaty! or 


= 


Here assuming y 2 0, we obtain the previous relationships (9.03) and 
(9.05). 


In the directions 8=y and ®8=s—7 (with o=s), where the 
scattering dlagram has a principal maximum, it follows from Equations 
(12.11) and (12.12) that 


Si) =S(1) = — kacosy (12.13) 
and 


Se — y= —E(x--- 1) = — kacosy. (12.14) 


In the direction toward the source (@=2—y. f=='5)s the functions 
Si) and 5(8) take the values 


£8 
oy eH H, (®). (12.15) 


Here considering f==2, we obtain 


S(r):-— —S(r) = ka, (12.16) 
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which corresponds to the physical optics approach [Equation (7.20)]. 


The functions £ and £ allow one to calculate 


open lElt, 2,005)? (12.17) 


which are the effective scattering surfaces with the E- and H-polar- 
izations of the incident wave. Let us recall that, by definition, 
the effective scattering surface is a quantity equal to 


om dent! 3h (12.18) 


where 


S= 4, RefEH’ (12.19) 


which is the energy flux density averaged over one oscillation cycle 
(the Poynting vector) in the scattered wave, and So is a similar 
quantity for the incident wave. 


In this way, we obtained the expressions for the fringing field 
which approximately take into account the nonuniform part of the 
current. In the incident plane (e=*+) , they have a form which is 
rather simple and convenient for calculations. It is also interesting 
that in this case they satisfy the reciprocity principle as distinct 
from expressions (10.03) and (10.04) which correspond to the uniform 
part of the current. It is not difficult to prove this by verifying 
that Equations (12.11) are not changed with the simultaneous replace- 
ment of y by ® and of ® by y, and Equations (12.12) are net chanred 
with the replacement of 9 by m - y and of y by e—0 [in the case 


g=-> }. 


However, Equations (12.11) and (12.12) lead to a discontinuity 


of the magnetic field tangential component Hy on the plane z = 0 in 
which the disk lies. As in the case of diffraction by a strip, the 
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reason for this is that we did not consider the interaction of the 
edges. It is also necessary to take into account this interaction 
in the case of glancing incidence of the plane wave (T=-}-), when 


the fringing field components E, and H, must be equal to zero. 


¢ 
Let us point out once again in conelusion to this section that 

expressions (12.11) and (12.12) near directions O=y, Ome #©~—-7 (with 

225) have an interpolation character, but in return they allow 

one to represent the fringing field in the incident plane x 2 0 in 

a convenient (uniform) form which frequently is of greatest importance 

(compare § 24). 
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CHAPTER III 


DIFFRACTION BY A FINITE LENGTH CYLINDER 


The distinctive feature of this probiem is that, in addition 
to the nonuniform part of the current on the cylinder's surface 
which is caused by the discontinuity, there also exists a nonuniform 
part of the current arising as a consequence of the smooth curve of 
the surface. This part of the current has the character of waves 
travelling over the cylindrical surface along geodesic lines [36] — 
that is, along spirals on the cylinder. These waves, which as they 
move strixe the edge of the cylinder, undergo diffraction and 
excite secondary surface currents. In turn, the nonuniform part of 
the current resulting from the discontinuity undergoes diffraction 
while being propagated over the cylindrical surface. It is clear 
that specific consideration of all these effects is a verv comriicated 
problem. 


However, if all the linear dimensions of the cylinder are 
sufficiently large in comparison with the wavelenrth, these effects 
may be neglected when calculating the fringing field in many cases 
which are of practical interest. In particular, they may be nerlected 
when calculating the field scattered in the direction toward the 
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source [5, 37]. In this case it 1s sufficient to consider only the 
nonuniform part of the current which is caused by the discontinuity 
of the surface, and we will do this in this Chapter. The equations 
obtained in this way are generalized to the case when the observation 
direction does not colncide with the direction of the source. 


§ 13. The Physical Optics Approach 


Let us investigate the diffraction of plane electromagnetic wave 


E= E, eittvetn T+ecoe 7 ( 1 3 .O1 ) 


on a finite, ideally conducting cylinder of radius a and length 7. Let 
us position the spherical coordinate system in such a way that its 
origin is at the center of the cylinder, and the normal n to the 
incident wave front lies in the half-plane 9= 5 and forms an angle 

" (0<1<+) with the z axis (Figure 24). 


An incident wave having an 
arbitrary linear polarizaton always 
may be represented as the sum of 
two waves with mutually perpen- 
dicular polarizations. Therefore, 
for a complete solution of the 
problem, it is sufficient to in- 


Figure 24. Diffraction of a plane vestigate two particular tases of 
wave by a finite cylinder. n is ; 
the normal to the incident wave 
front. 


incident wave polarization: 


(1) E-polarization, when the incident wave electric veeter is 
perpendicular to the plane (Eo _L yoz)and 


(2) H-polarization, when HoLyoz . 


The uniform part of the current excited on the evlindrical sur- 
face hy wave (13.01) has, with the F-polarization, the compene its 
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f=— = £,,:sinysinge™, 


c 


= £ E,.sinzcospe™, (13.02) 
; fm % Exa-cosycosye™,: 
and with the H-polarization it has the components 
L=f=0 | 
f= — gy Heasinge’™®, | meee 
where 


b = asin ysin 9 + € cosy. (13.04) 


Let us calculate the fieia created by these currents in the region 
n 
sdeleaae 


The vector potential of the fringing field is determined by 
the equations 


ae 
4 
A=t [dy { re. ya with 7=0 (13.05) 
a oak 
and 
is 
eee ae eit 
are) dy J #6 de with y>0, (13.06) 
T 
where 
r=/Pfoy—y+rae-o . (13.97) 


Since the field in the far zone (R> ka’, R >a) is of interest to us, 
these expressions may be simplified by using the relationshi!n 


r=R+a sin) sins ~€ cos 8. (13.08) 
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As a result, we obtain a simpler equation 


U 


r . 
an , 
As tS [eteuntonegy fe werner, (13.09) 
t 
“Tt 


Since the current components are described by the functions 
ie". then the problem of finding the field reduces essentially 
to a calculation of integrals of the type 

ui 


3 Q 
i) eid (cos 1 — cs * ae j Piper 4 dy = 
é —s 


2 
[ IF (cos Tow &) _ 


= jk (608 y— cos) (13.10) 
- Me cvs O : . 
- ‘5 (cos y—cve ‘| f fcpeleni dy, 
The integral 
2 } 
J fever say, p=ka(sin {+sin 9) (13.11) 
—s 


when p >> 1 is easily calculated by the Stationary nhase method. The 
stationary phase point is determined from the condition sing == 
and equals 


~p=— > * 


(13.12) 
Then assuming y=—-3-+8, we find 
° r 
Jimermrayxsae’ [e'F a= 
= -$ (13.13) 
. oe 7 -1 +5. . 
=A Frome {ena V Fre! 


AS a result, we obtain the following expressions for the 
vector potential: 
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with E-polarization 
5S eye 
2-3 dint Bnd, A, =A, =0 (13.14) 


and with H-polarization 


Ang Mad Ap Ay =O (13.15) 


where 


iF eos Woe re “i oT (cus y—cos 8) 


e 
fs hea — vos 7) x 


4 ~te tn rosin ty +t (13.16) 
xV da(elay fein dy © ° 


The fringing field in the region 9=—y is determined by the 
relationships 


E.=—H = ikA,, 
. | (13.17) 


B= H, = —ikA, sind. 
Therefore, with the E~polarization it equals 


ia 
E,=—H,= Bt E,ssiny G1 


(13.18) 
£,=H, =0, I. 


and with H-polarization 


ike 
E,=4, =— Fe Ms sind I, 
£,=H,=0. 


The resulting equations show that the field scattered by the cylin- 
drical surface 1s created mainly by a luminous band adjacent to the 
cylinder's generatrix with q=H=— > . The radiation from this 

band may be represented [see Equation (13.16)] in the form of spheri- 


cal waves diverging from its ends (points 2 and 3 in FPirure 24). 


Now let us write Expressions (13.18) and (13.19) in a form which 
is most convenient for calculating the effective scatterinr area 
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iAR gad 
y= Ha BP 0 1h 
Y.' ") 


(13.20) 
lth @ 
EH, =F He )) ® - 
Here 
J easing, YY = —Gsind 
oa re (13.21) 
and 


~ « 2 
G =2 V sears 5) x 
sia (eon — cos 0) —iha (ain ¢+sin 41> (13.22) 
x cos ¥ — cos 7 . 


The index "0" on y. and x means that the field was calculated in 
the physical optics approach (based on the uniform part of the 
current), and the index "c" shows that this fringing field is created 
by a cylindrical surface. The effective scattering area, in accord= 
ance with (12.17), 1s determined for a cylindrical surface by the 
relationships 


oot Sa aatsintr ia, | 


of y= ra"| S| nat sin? (Gl | (13.23) 
In the direction of the mirror-reflected ray (®#=y), we have 
ert a halen ba Si aint, (13.74) 


In the direction toward the source (8 =<—y) , the functions ~ 
C 
and i", equal 
tika sin O+i— 


aly = sin} sin(k/ cus #) 7 a 
5. =-\i= =f. “cost ae s (13.25) 


These expressions are valid if ka sin @>!1. [It is not diffieult 
to see, by means of equations (13.02) - (13.05), that the frimeini 
field equals zero if y = 0 and =a. Thus in the case of radav (that 
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is, in the direction toward the source) we find an expression for 

the fringing field in the region ka sin 9%J anc in the direction 
@=a1, Naturally the desire arises to write interpolation equations — 
that 1s, equations which would provide a continuous transition from 
the region ka sin8%1 to the direction g=x. Now let us note thit 
the field scattered by a cylinder is comprised of the fields scattered 
by the lateral (cylindrical) surface and the base (end) of the 
cylinder. In the pnhysical optics approach, the fleld scattered by 

the end of the cylinder 1s equivalent to the fleld scattered by 2 
disk. But the field scattered by a disk is described by Bessel 
functions. Therefore, it is also advisable to express the field 
scattered by the cylindrical surface in terms of Bessel functions. 

As s result, the desired interpolation equations for the field 
scattered by the cylindrical surface may be represented in the form 


x ee y=-- aL i) (elttecs O_ eM Ty (Cy id, (6)], 
(13.26) 


C= 2kasin 8. 


From this it follows that 209 “le = 0 in the direction @=«, and with 


the conditions ka sin ®%1 we obtain Equations (13.25). 


The field being scattered by the cylinder's end (by the disk), 
in accordance with equalities (10.06) and (10.07), is descrited in 
the physical ontics approach ty the equations 





2a | assim (13.27) 


Consequently, the field scattered by the entire surface of thie 


cylinder will be determined in the plane e=—- > by the equations: 


i (eR exe 
E,= Ay =S > Eg )) (5), | 
if eth? rye 
=H =F Me cy (9), | 


~ 
| al 
wd 
e 

ro 
™~ 
— 


where 
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py (®) sin’ 


=> (el* vos8 ei cus *) lJ, (e) be id, (6)] = 


° 2 cosa 
yi) (13.29) 


t fate ,, . 
saerAmer® (&- 2xasind). 





Tnese equations allow one to determine in the physical optics 
approach the effective scattering area of a finite cylinsler. 


§ 14, The Fleld Created by the Nonuniform 


Part of the Current 


Let us find the field from the nonuniform part of the current 
caused by the surface's discontinuity. Figuratively speaking, the 
field scattered by the cylinder is created by the "luminous" regions 
cn its end and lateral surface. Mathematically this field is 
described by the sum of spherical waves from the "luminous" points l, 
2 and 3 (see Figure 24). Obviously the field from the nonuniform 
part of the current also will have the form of spherical waves diverg- 
ing from these same points. 


In the case when the length and diameter of the cylinder are 
sufficiently large in comparison with the wavelength, one may 
approximately consider that the nonuniform part of the current near 
the discontinuity is the same as that on a corresponding wedge. The 
field radiated by this part of the current in principle may be found 
in the same way as in the case of the disk. However, such a method 
is rather complicated. We will find the desired field by a simpler 
and more graphic method, starting from a physical analysis of the 
solution obtained for the disk. 


For this purpose, let us investigate the structure of waves 
(12.01) and (12.02) which are radiated by the disk. These equations 
include the factor 

ia lan ith i 1 
Viska@iny tia) | Vou VR isa (14.01) 
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Here Vi 


is the unfoldiny, coefficient of the wave. It shows how 
the fleld is formed with increasing distance from the disk: the 
agiffracted wave which is cylindrical near the disk unfolds intn 2 


spherical wave as the distance from it increases. The coeffictlent 
(sing -{-sin 8)-"s is proportional to the width of the luminous reston on 
the disk or, in other words, to the width of the first Fresnel zone. 
Thus, tn Equaticns (12.01) and (12.12) the functions f! and g! depend 
only on the body's geometry — more precisely, on the character of 
the discontinuity. 


Therefore, it is entirely natural to assume that the similar 
waves which are being scattered by a cylinder have the same structure 
and differ only in the functions rl and gt which correspond in this 
case to a rectangular wedge. Consequently, in the direction toward 
the source, the field from a nonuniform part of the current flowing 
on the cylinder may be represcnted 


when ka sin @>!f in tne followinr 


way: 


: (:-¥) + tht con® 


=H waite Ba fraye | 


po (i? (2) eltteand +i (3) e -ihleos se ~1 (7 Phe ltt 


E,=H, = ag Ho te’ (I)e « ssid 


a 
e 
; (18.02) 


pam {g (2) elt 4g (3) e7 ‘eo Se to), : 


1 


In accordance with § 4, and rt are determined by 


the functions fr 
the equations 











to— 
ayy fe =" ( —— 
ey e ar ; 
Pa ! = coe @ 
is cos 2 = con alee! (14.03) 
8 
rey) _, onl sets fil )= 
g' (2) cos I cos cone 
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(Equation continued on next race. 


81 


_. co98 __ sio& 
~“Qgine + 2cos8 ° 


P 3) 4 ue ( ee eee 





. 
. 
, Sp ee? EE CS 


g& *(3)) @—~.—} 
(14.03) 
os ! lob 
pala) ait ' 
cos — — ene n 
where 
3 
a (14.04) 


In Chapter IV, we will show [see Equation (17.25)] that in the 
direction #=a—y=# one may neglect the fleld from the nonuniform 
part of the current flowing on the cylinder in compartson with the 
field from the uniform part, tfka>I. Therefore, for the field 
from the nonuniferm part of the current, one may write with the help 
of Bessel functions the following interpolation equations: 


aR = 
a ee . wk (8), 
(14.05) 


YS" (8) = LAT, (0) EIN, ype“? — 
—PALO—i Qe *, 
Si (8) = iat", () ANY, (er? — (14.06) 
= 8' (3) [AG —U, en, 


and the functions Mr, v2 and Mm), nt respectively equal 


M 


it SR -f"(2), Mo = (1) =F g" e (14.07) 


= =| 


or 
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a 
mL sia-y _ ' 
M (a | eat 











8 w—?2 
cos ~~ cop 
' cose? sink 
ar ao). iin’ * Seon * 
cos —- — cos—— 
wi ‘ita 
N \- 2(1= (14.08) 
N’ " coo—- = 1 
a 
a3 ZZ 3s 
rr rr Py | iad 
cos =~ — cos —— 
= u sind 
‘s 28 e con 8 
roa — con 


The resulting Equations (14.05) change when ka sin®>I into Equations 
(14.02), and in the direction ®=* they give a value equal to zero 
for the field. \ 


§ 15. The Total Fringing Field 


Summing Expressions (13.28) and (14.05), it is not difficult to 
see that the total field scattered by a cylinder will equal 
5 a 
E, = —Hy =} Eva SY). 
aR 15.01 
Ey =H, = "2 Ha SY. C ? ) 


where 


— f (3) [4, 0) — #49] ieee ° 
»y (9) = (MJ, (2) EEN, (6 elt? — 


—g (3) 14, (6) — iJ, (Spe, 
C= 2kasind . 


J; 8) = 1M, (0) + INS, (0) et? — | 
| (15.02) 
J 


and the functtons M, N and M, N are expressed only in terms of the 
functions f and g which correspond to the asymptotic solution for a 
rectangular wedge 
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=) = 


pros == / (2), who =g (l= (2). (15.03) 


or 








M85 ae ’ oe ’ 
M{™ a ( = a a—— # 6) * 
coe -—-— con cos cos 


(15.04) 





a acre gee eee 


——— | cos ¢03 ——~ =" cos 


The functions f(3) and ¢(3) in turn are determined by the equation 








132) wp, a 
g3)}~ 4 ( See eee aa): (15.05) 


e 
cos —} _ 608 — C08 


Thus only the functions f and g are included in the final expression: 
for the scattering characteristic of a plane wave by a cylinder. 


In the direction %=2, the functions 2(6) and sj take, as in 
the case of a disk, the values 


E (2) =: — 5 (2) == hae, (15.06) 


and with d=: 5 they respectively equal 


(15.07) 





4 ow 3 
Y(=)= [s me — eg cshan (+. 
cos 


+(4 ctg Z+itl) AO, 


where ¢ 2 2ka. The terms in this equation which contain the factor 
kl refer to the field from the uniform part of the current, and the 
remaining terms refer to the field from the nonuniform part cof the 
current. 


FTD-HC~-2 3-259-71 


In accordance with (12.17), the effective seattering area of 
the eylinder is determined with the E-polarization of the incident 
wave by the function 


a, == 0a"/E(8)|", (15.08) 
and with the H-polarization of the incident wave by the function. 
9, = aa" | &(8))’. (15.09) 


Let us note that Fxpressitons (15.02) for the seattering field 
may be obtained directly on the basis of an analogy with the 
Equations (12.06), omitting the calculation of the fields from the 
uniform and nonuniform parts of the current. In the same way, one 
may obtain the expressions 


_ a 
im 


Z(®, %) = (MJ, ()+INS, Je’ 


1 $F (eon 6 + con by) 


—/A)4,0—Mm@le ; 


(cos © + cos 8) 


9 eg (15.10) 
£(0, 0) =[M/,()+iN,QJer 


eM con® + ens 
— 2), 0— Weyer 


which are suitable for calculating the fringing field in the region 
2 


!=—F. F<h O,< (0,=2—1), The quantities here equal 











§= ka (sin 8 +- sin 9,) (15.12) 
ea 
2 M ein 
uo aX 
= t 
si bg = n—0—6, ~~ ® =m 
cos ~_~— coe n cos -~—"Ccos 
N Ay 2 -_ (15.12) 
oor (qEctesn 
cos =~ cos 
85 
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£ 
F(3) ea ' ae ! 
ge peo =i se athe’ 4 
a(3) " (nee * oceeencars) (15-13) 
Expressions (15.10) satisfy tne reciprocity principle — that 


ls, they do not ¢hange their values if one interchanires ® and %. 
When @=8, they change into tne previous Expressions (15.02). 


Equations (15.02) and (15.10) describe the radiation from the 
currents flowing’ only on part of the cylinder's surface: on vhe one 
end (when z 2 By and on half of the lateral surface (-" < » < 0). 
Moreover, these expressions do not take into account the nonuniform 
part of the current caused by the curvature of the cylindrical sur- 
face. Therefore, they must be refined with values of 8 and & 
which are close to - and ns. However, in the case 4 =f, — that is, in 
tne direction towards the source — these corrections may be nerlected 
if the parameters ka and ki are sufficiently large. ‘Xumerical cal- 
culations performed by us on the basis ef Equation. (15.02) show 
that this evidently may te done already when ka = 1 and k? = 10 7. 

The graphs of the functions HE, = E(b))? and oh, ==) Ei) * constructed for 
this case in Figures 25 and 26 agree with the experimental curve ‘}) 
(the dashed line): the position of the maxima and minima basically 
agree, and the number of diffraction fringes is the same. For the 
purpose of illustrating the effect of the ends, we constructed a raph 
of the effective scattering area for those same values of ka and kl 
taking into account only the uniform part of the current on the 


3 


A 


cylindrical surface (Figure 27). A comparison of Firures 25, ac and 
27 shows that tne effect of the ends begins to aprear when &® = 120°. 


Footnote (1) appears on page 89. 
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Figure 27. The effective scattering area 
of the lateral surface of a cylinder in 
the physical cptics approach [see (13.26)]. 
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1. on pare 86. The experimental curves shown in Figures 25 and 26 
and also those in Figures 31, 32, 65 and 71 were 
obtained by Ye. N. Mayzels and L. S. Chugunova. 
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CHAPTER IV 


in this Chapter we will refine the physical crties arrreach 
for certain ctheriodies of rotation, whose surface: have cireoular 
discontinuities. We wlll limit curselves te the case when a plane 
electrcemagnetic wave falls on the todfes alens thetyr cyvmmetry axts 


nz befere, we will assume that the linear dimenctions <f the 
bedies are large in cemparison with the wavelencthn, In this care *: 
currents in the vicinity of eireular disecntinulty of anv eonvex 
surfs 


e of retaticn may be aprpreximately considered to be the same 
crrespondinzg conical tedyv. Consetuentiv, it 


u 
hich £5 caused ty the circular diseccntinulty of the aur 
u 


ad Created Cy the Nonunt form 


1 
Part of the Current 


Let a plane electromagnetic wave fall on a conical bedy in th 





tent to study the field from the nonuniform part of the curren 
£ 


oO 


Pend 


1 
resitive direction cf the z axis (Firure 28), Prem the relationch! 


B= —-) (graddivA+#A), | 


‘ (16.01) 
H==rotA | 
we find the following expression. fer the fringing fleld in the wave 


Zones 


22=H,=ikA,, \ ree ee ee 

E,«—HeaikAy § "°"? (16.02) 
and 

£,=—H,=ikA,, ‘ e— 

E,=H,=ikA, \uttn . (16.03) 


The vector potential is determined by the equation 
a : 


A:= + 7] [fs (pee ™ “(a —Ssin to) do 


1€.94) 
4 m ( / 
+} Jaye 8 (a — Sin (2) dl] dy. 


Here r is the distance from the discontinuity to the observation 
point, J, 05) {s the surface current density flowing on the trradiated 
side of the body, and Jos) is the current density on the shadowed 
side. The upper sign in the exponents refers to the case 8m, and 
the lower sign refers to the case ®8=0. Since the nonuniform rart of 
the current is concentrated mainly in the vicinity of the discontin- 


uity, the vectcr potential ccrrespending to it may be represented in 
the form 


A= a (j Htpet ae + 
+f Oe? *2) dy. Ree 
) 


Obviously the nonuniform part of the current near the discontin- 
ulty of a conical surface may be considered to be approximately the 
same aS on a corresponding wedye (Figure 29). In the local cylindri- 
cal coordinate system Tis o> Zi» the field from the nonuniform part 
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Figure 29.0 The dthedral ant le 
corresponding to the diseont ine 


Fipure 28. Df ¢fraction of a uity of a confeal surtace. 


plane wave by oa conleal body. 
The plane wave ts propayated 


the ppent ff soon sue 
alone the 2 axis. of the current flowing on such 


a wedre ts determined tn the 
far vone (kry >> 1) by the followtny equations: 


oo? | 


Mia boty EV Ga | (16,06) 


where 


Ae. / te [{ I aad 


. " aie (16.07) 
+ fire rae ax I 


Here the upper ofan fn the exponents: refers te the ease ay son tiwy 
und the lower siym -— to the ease o) aw. On the other hand, In 


gs 4 it was shown that this fleld equals 


ee Clee) 
Hy MWR! Se oes 
: ork, 


where Ey. a). Ay, uy) tre the values ef the Inefdent wave amplitude at 
« oa? 
Jd | . ; 
the wedire edge, cid fo and et are anyular funetdens eharneterts fae 


tie seattertoe diagram. 


adler eee a 


Let us introduce the designation 


. : - iit cos 
1) ieet™ a Hee? ty. (16.09) 


Equating Expressions (16.06) and (16.08), we find 


J = fF, py 


= eH g,, (4) 
hy 1k2a jaa 


t= Tae (16.10) 
The components Jz, and Jo, are mutually perpendicular, and when 

®=0 and ®8=s they are parallel to the plane xOy (Figure 30). The 
different orientation of the unit vector ee, when 8=0 and @=«a Is 
connected with the fact that the angle oy 4s measured from the trrad- 
lated face of the wedge. In the original x, y, z coordinate system 
the vector J has the components 


J,=4,si09—J,cos¢, | a=0 (16.11) 
y=— J, cosy —J, sing yay 7 


and 


with ®=«. (16.12) 


J,=-J, sing+J, cosy, 
Jy=— J, cose+ J, sing } 


Substituting Expressions (16.10) here, we obtain 
Ja= Faq (P Eee, (9) sin» — gH, (4)c08 4], 


with 8=08 
Jy = Fj lI "Eon, (4) 008 9 +.2'M,, (9) si0 $] : (16.13) 


and 


Jn = ayy UE os, (SiN 4 + BHg,, (9) £084), 


itr o-- a ‘ 1 
Joa phic WE yc <M, sing) [EE OE OSU) 


woe yee ee 


Now identifying the current near the conical surface discontinulty 
with the current on the wedre, we find the components of vector 
potential (16.05) 
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Figure 30. 
orientation of the 


The relative 


a eres rey 
eS ee 7 fire 
rf 


“us, (ppsin gy ~ 
ng (Hew 2h, 
A, =— f° “efie E,, (phous bf 


Fu (9) sin 9] hp 


with @ =0 (16.15) 


0 ee © eee oo 


and 


unit vectors Cg) and 


ez) in the cases 
and d=., 


Furthermore, let the Plane wave be 


Eg! |ox. Then 


Eu, ($) = E,. sing, Ay, (¢) = — E,, cosy. 


60 


A, =: Tor wu, lysing 4 


+g! 


. @ 
Ay = age 


(¢) cos 2}, 
ia ‘a with ®=® (16,16) 
aae patate 
e frente cose. “ 

C) 
=a", (sind, 


Polarized in such a way that 


(16.17) 


Considering these relationships and substituting Expressions (16.15) 
and (16.16) into Equations (16. 02) and (16, 03), we find the fteld 


from the nonuniform part of the current which 


1s caused by the 


circular discontinuity of the conical surface 


and 


Equation (16.18) is 


Equation (16.19) for the values 


PTD-HC=2 3225927] 


FE 7 H, POA 3 -(f pers 
= ne * | with ®=0 (16.18) 
Fi =>H,=0 | ! 
Ff, =: H oa wee ( tr. yer 
b 9 ? & ’ ten Mss: ge, (16.19) 
Ey >H,=0 f. : 


applicable for the values OSS Fes OS2 , and 


0<e<t .O<SHese. In the ease of 4 


94 


td 


disk(#- >. # ), the Pleld frem the nonuntform part of the current 
J ] 
4 


equals zero on the 7 axte, cinee £o = 2° @ =1/2 when 9—0, and 
pho tw -2/2 when @aae feotpure (E11 

Using the result ine ce ot os onfi: fs the following seetions,we 
will calculate tae 2 ff toe ‘corte. area (in the direction 4=« ) 
for specif!le btiudfes. oa that they are lrradftated by 


the plane wives 
Fe, WA Fie. (16.20) 


and their linear dimenctoans are lorce fn comparison with the wave- 
‘leneth. 


§ 17. A Cone 


Let a cone (Figure 28) be trradiated by plane electromarnetic 
wave (16.20). The uniform part of the current which is excited on 
the cone's surface has the components 


f=m £,, sine! *?, 1 
f=% on (17.01) 
f= % E,. cose cosge’* 
and creates in the direction =a (with R>-ha’, R>RP) the fleld 
; ge 
2,2 —H, = — Eva tee + 


+E ve (A ete t tge)ty etm, (17.02) 
F, = H,=0. 


Here the first term describes the spherical wave divereing from the 
vertex of the cone, and the remaining terms describe the spherteal 
wave from its base. 


Tne field caused by the discontinuity of the surface at tne 
cone base is a spherical wave, and is determined in accordance wit!: 
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(16.19) by the expression 


3 at iar 
a "et ' 
E,=—H,=- $n toot a ic erie (17.03) 
] 
| 
J 


cos a cos . 


£,=H,=0, 


where 


naipets, (17.04) 
An asymptotic calculaticn of the rigorous diffraction sertes 
for a semi-infinite cone [38-40] shows that in the directionO=a 
one may neglect the effect of the nonuniform part of the current 
caused by the conical point. Therefore, summing (17.02) and (17.03), 
we obtain the following expression for the fringing field: 


E,=—-H,= ~uf tet (1 — etary 


2 s = 
a a ia (17.05) 
+ ka enar fo 
con e082 | 
8, >H,=0. 


Let us point out the following important feature of the resulting 
equation. In the problems which were investifrated in the prevtous 
chapters, the edge waves of the fringing field were expressed only 
in terms of the functions f and ¢. But now intie equation for the 
spherical wave from the cone's base, in addition to the term which 
depends on f and ¢ [the last term in the bracket of Fquation (17.05)], 
there is an additional term [term --@/te’me?"" tn Equation (17.05) ] 
which does not depend on these functions and is determined by the 
uniform part of the current. Therefore, it 1s impossible to rerresent 
the resulting spherical wave from the cone's base only in terms of 
the functions f and g which characterize the total edre wave dlarram 
from the corresponding wedge edge. This important fact was not 
considered in (41, 44], as a consequence of which their authors di!d 
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not succeed in obtaining correct results for a cone with an arbitrary 
aperture angle u(0 «< w « n/2). 


The effective scattering area in accordance with (12.18) ts 
determined by the equation 


a= e2(Z/°, (17.06) 


where the function £ is connected with the fringing field by the 
relationship 


Ey =H, =~ $B te 
oN, = g Con BR (17.07) 
and equals 
2 
== Ligtesin ki eitt + ~ "5 the, (17.08) 
cos cos 


The analogous function in the physical optics approach may be written 
in accordance with (17.02) in the form 


Beas tygtesin lem —tgaetn, (17.09) 


With the deforming of the top part of the cone into a disk 
a 
(o> F- /=9), Equations (17.08) and (17.09) are transformed, respec- 
tively, to the form 


f= —tha~ + ctg—, ) 


TP = — jha. 


(17.10) 


Furthermore, it follows from (17.08) and (17.09) that for larre 
values of the parameter ka(ka >> teow) the functicns £ and a may 


be represented in the form 


= sin — (17.17) 
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2 = tgw eres, (17.12) 


Thus even in the case of short waves ( ka >> teu, but R >> ae 
our Expression (17.08) does not change into the physical optics 
equation, but substantially differs from it because 





3 
— sia -— 
pers x A = (17.13) 
cos ne 
and 
3° == ea" tg*e. (17.14) 
With this 
1 ~ sin = i 
an a 
1=2| pope “tay ’ (17.15) 
! 


that is, for sufficiently short waves (or for sufficiently large 
dimensions of the cone) the function o is proportional to 0°, The 
coefficient of proportionality here does not depend on the cone 


dimensions, but. is determined only by its shape. 


This result is graphically illustrated by the curves riving 
the effective scattering area of a cone (w = 10°25', k = nr, 2 = 90°) 
as a function of its length (Figure 31). Whereas our equation (the 
continuous line) is in satisfactory agreement with the results of 
measurements (the small crosses) ‘2) | the physical optics approach 
(the dashed line) gives values which are smaller than the experimental 
values by 13-15 dB. For sharply pointed cones, the nonuniform part of 
the current has an especially large values. In Figure 32, a curve is 
constructed for the effective surface of a cone (ka = 2.75 7, 2 = 90°) 


Footnote (1) appears on page 113. 
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with its deformation into a disk 
(w + 90°). The discrepancy. be- 
tween our curve and the physical 
optics approach here reaches 
almost 30 dB when w = 2°. 


Expression (17.08) obtained 
by us also allows one, ln con- 
trast to the physical aptics 
approach (17.09), to evaluate 
the role of the shape of the 
shadowed part of the body and 
shows that the reflected signal 


Figure 31. The effective scatter- Will be larger, the closer this 





ing area of a finite cone as a shape 1s to a funnel-shaped form 
function of its length. The 

function 3g (the continuous line) (Q = ~ uw). Thus, for example, 
was cece ee the basis of in the case wu 2 10°, ki = 10 r 
equation (17.06) which con- 

siders the nonuniform part of (k = 7) the sipmal reflected by 
the current in the vicinity of the cone may exceed by 15 dB the 
the circular discontinuity. 

The function 0° (the Aeenca value corresponding to physical 
line) correspcnds to the physi- optics (see Fipure 33) if 2 = 


cal optics approach. 170°. 


Let us note that our Expression (17.13) 41s equivalent to the 
expression presented in the above-mentioned papers [41, 44]. Hewvever, 
the latter expression is applicable only for sharply pointed cones, 
whereas we have, in addition to (17.13), Equation (17.08) which is 
suitable for cones with any aperture angle w (M fes op 

The calculation method discussed may be peneralized in the ense 
of asymmetric irradiation of the cone. However, with asymmetric 
irradiation, generally speaking, it is necessary to take into account 
the nonuniform part of the current caused by the point of the cecne. 


In concluding this section, let us calculate the effective 
scattering area for a body which is formed by rotation arcuni the 
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Figure 32. The effective scattering area of a finite cone 
as a function of the vertex angle. 


z axis of the plane figure shown 
in Figure 34. Integrating the 
uniform part of the current, it is 
not difficult to show that the 
field scattered in the direction 
d=" by the lateral surface of the 
VORB eR Ree eee ee truncated cone (Figure 35) is 


ee 
Ht tT determined by the equation 


E, ioe ha: H, a Lys: =i tg? a, -}- 3 tg “) em + 





Figure 33. The effective scatter- oe \ cael. [elk 
ing area of a finite cone as a +( ten byte; & oth iy 
function of the shape of the ; ; 
shaded part. (17.16) 


FTD-HC=2 3-259-71 toe 





Figure 34. The generatrix of Figure 35. The generatrix of 
the surface of rotation. of a truncated conical surface. 


Summing this expressing with (17.02), where the quantities 7% and a 
must be replaced by Ly and aj,» we find the field from the uniform 

part of the current flowing on the entire illuminated side of the 

body 


£,=—H,=— ce ie ty? wsin ki, eo — tge eh 


+(e: tg?o, sin kl, ef 4. (1 _- a eM ltg o, |e \ ae ‘ 


(17.17) 


The field radiated by the nonuniform part of the current is 
determined in accordance with $ 16 by the equation 


vip / asin ai 
E,=—Hy= — 83" (oe +tge—tgeyem+ 


EN ti 
where 
maIte eM al p Ee, (17.19) 


Now summing (17.17) and (17.18), we obtain a refined expression 
for the field scattered in the direction 8=--s 


penta igteainaler or (17.70) 


(Equation continued on next pare.) 
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2 
= sin-- 
al aint, ! ‘ jig & ely 
bam ag Oh gg, tems sin Ae 8 OF 
cos ay on 
2s (17.20): 
— sin — i ° 
a 4° mu | eit +0 tt 
a 2u, Ye 
cos 5 —~ Coa 
e a 
Consequently, the effective scattering ares will equal 
pres) 
2) , ial, my Ms aikt 
a== RQ) ga, '8 sin ke * of --—= see + 
cos 5 ~~ cos 
1} t 
: t) (17.21) 
I ak, , @ nla, Qik (f+, 
+5, te? sir biel Mert 4 ns Qt 
L * cos .- — cos — 
Ly m 


In the physical optics approach, the analogous quantity equals 
of = za; Fe tytasin kle™” — tg we 4. 
t] 
1 . : i (17.22) 
+[is; tg", sin ki,e" + (1 a0 enn tg «| et | ‘ 


\ 
When the top part of the cone is deformed into a disk (e+%- i, +0) : 
Equations (17.21) and (17.22) take the form 


2 : ! . j 
= ea} | —ika,-- x cg. tas tg? sin kl, -} 


2 I) 
tn. Sh (17.23) 
+% i an eikle : 
af cus = — cos 28) 
My 


so = ra} ! — ika, + as ig? a, sin kle™ — 


+(1— je) tgs, |. (17.24) 


In these expressions assuming w, * 0, we find the effective 
‘scattering area for a finite cylinier 
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4 . t = my Ay sit 
== RQ" |--ita,— -ct --o-— ey, 
37: 8A) ira, mB aT ee (17.25) 
ur) 
3.2 aa’ (ka,)’, (17.26) 
in connection with which 
3 2 
A=yeA=l-f oS. (17.27) 


Equation (17.25) 1s more precise than Equation (15.06) which 
was derived in § 15, where the value of the field itn the direction 
0=x was taken in the physical optics approach. 


§ 18. A Paraboloid of Rotation 


Let us calculate the effective scattering area of a paraboloid 
of rotation r? = 2pz (Figure 36) which is irradiated by plane wave 
(16.20). The uniform part of the current excited on the parabcloid's 


surface has the components 


= 5 Eessinae™, 
i= (18.01) 
i = Ey + COS 2 COS szelM, ) 


Integrating this current, it is not difficult to show that in the 
direction =e it radiates the field 


E,=--Hy= — Eqs ral =e) tg 0S"). (18.02) 
E,==H,=0 j eo 
Here a ts tne radius of the base of the paraboloid; toa = clgw is 


{ts length; a is the angle between the z axis and the tangent to the 
generatrix of the paraboloid (r? = 2pz ). At the point z= 7, the 


angle a takes the value a=e/(tge= 2) ‘ 
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{Site 


Fissure 36. Ciffractioan 
of a plane wave fe 4 
paratolcid of rvtaticn. 


The ffletd from the nonunifern mart 
of the current eaused by the discontinuity 
of 


Ly Equation (17.02). 


the parstolcid's surface tu determined 


The ffeld from the 


nonuniform part cf the current which is 


eaused by the smaotn curve af the parne- 
surface eyusls zero itn the case 


boletits 
of symmetric radiation [45]. ‘Therefore, 
summins (18.02) and (17.03) we find the 


expression for the resu:ting fringing 


field. 
-sin- = \ ing 
E,=—Hy= — te tga-- ee ieee ee 
. cos — cos 


Conseaquentiy, the 
will be determined by the relationship 





2 ie « 2 

= sin -- 

a a f 
srr@ tga =~ ae ent 


cas a7 cos Tn j 


(18.03) 


effective scattering area of the parsboloid 


(16.04) 


which, when the paraboloid is deformed int. a disk (m+ 7,1 -+0, 0 


the form 


const} » is transformed te 


" ! 7 
a=: ra? |ika -|- nce x | 


Comparing Excressicn (18.94) with the equation 


a°-= rat tg*w fl —e it, 


which physical optics siv for the effective 
see that they differ sianificantly from oné another. 


of the funetion a? 


First 


the oscillating character draws 


Seatterine are 


(18.05) 


(18.06) 


a, we 


of all, 


our attenticn: 


ected enal equals zero tf a whele number of halfevaves 
' 


e dost 
A --7) \ ‘ ; , 
gf Ae 1 By eS, tes Fitted inte the leneth of the garatceletd, and 
tC tances 4 maximum value Lf a half-{ntesral number of nalfesayes 

A 

2 


{ 
(! sy {at Pye al 2itice:) 1s contained in this lenrth. 


A calculation performed by us on the basis of Fyuation (16.04) 
for parabcloids with the pararneters 2 = 90°, teu = 9.1 (x = +) shows 
(Fifeure 37) that, although the ceetllatinz character cf the effective 
Scattering area is preserved, the amplitude of the oscillations 4s 
only 2 dB, and the maximum values of the function o exceed the corres- 
ponding values in the physical optics approach by almost 13 4E. A 
Still stronger divergerse between the results of our theory and 
Physical optics is detected when the paraboloid is deformed into a 
disk (Figure 36, ka = jn , k #7, 2 = 90°, w + 90°). 





Figure 37. The effective scattering area of a finite 
paraboloid as a function of its length with a constant 
value of the angle w (tgw = 0.1). The diameter of the 
base varies. 
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Pheupe 38, The effective ceatterine area of a 
finite paraboloid as a function of *ts leneth 
with a constant radius of the Lase. 


As in the case of a cone, the shape of the shadewed part turns 
out to have a substantial influence on the retleeted stenal. For 
exumple, for a paraboloid with the paramcters ka = on, ki = 10 a, 
tow = O.1 (k 


inerease of 2(w < 2 eo nm ew) CRPicuivs 49). 


nm). the reflected sienal tnereases by HE di with an 
’ 
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Figure 39. The effective scatter- 
ing area of a finite paraboloid 
as a function of the shape of 
tne shaded part. 


In concluding this section, 
let us dwell for a moment on 
the question of calculating the 
effective scattering area for 
bodies of rotation of a complex 
shape, whose elements arc the 
lateral surfaces of truncated 
The field from the 
nonuniform part of the current 
arising in the vicinity of 


paraboloids. 


circular discontinuities may be 
determined without difficulty 
from Equation (17.03). The field 
from the uniform part of the 
current is found by quadratures. 
Thus, the field being created in 
the direction t@=s= by the uniform 


part of the current which flows on the lateral surface of the truncated 
paraboloid r? = 2pz(p = a,tgu, = agtgus; see Figure 40) ts determined 


by the equation 


at. vp tf 
Ey=—Hy=— + Eye (ayty 0 -— azty ov, 07 yee, 


£,=H,=20 
Here 


' 
= 7 (aactg a, — a, ctg ») 


a (18.07) 


(18.08) 


{is the height of the truncated paraboloid (the distance between its 


bases). 
corollary of Expression (18.02): 


Let us ncte that Equation (18.07) is a simple alrebrate 
it is the difference of the fields 


scattered, respectively, by the paraboloid of heirht h+h= 5 and 


2 
by the paraboloid of heirht b= 35. 
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§ 19. A Spherlcal Surface 


§ <a 

Incident wave (16.20) excites a 
a surface current on the surface of an 

— ,ideally conducting sphere (a radius 

—le p and a center on the z axis at the 

point z = 9). The uniform part of 

Figure 40. The generatrix of this current has the components 
the lateral surface of a - ie 
truncated paraboloid of f= —§ E,scos Ge, | 
rotation. 

§=0, (19.01) 


j° =: 5 Ey, sin Ocosye'™. | 


The currents flowing on a spherical ring cut from the sphere's 
surface by the planes z 2 Ly and z = Ly + ls (Figure 41) create, in 
the direction 9==2, the field 


eH Ey [-(3 tg, —H)e eth + | 


+(S1¢¢,—4) aad ae | (19.02) 
£,=H,=0, 
where 
f,=p(1 —sin»,), } (19.03) 
t, = p(sin , — sin o,); 
P= assur eos a . (19.0%) 


Here ay is the radius of the first cross section; a, is the radius 


of the second cross section; wi (u.) ds the angle between the = axts 


and the tangent to the meridian at the potnt s = Lj (2, = ty + fy). 


Furthermore, assuming in Equation (19.02) w, 2:2 5 (, == const) , we 
obtain in the physical optics approach an expression for the fleld 


scattered by the spherical segment (Firure 42) 
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of 





Figure 42, A spherical segment 
with a conically shaped base. 
Figure 41. A ring cut from the 
surface of its sphere by the 


planes z 2 ty and z = ts + Lo. E,=—Hy=Evs (-sct4ht 
} ak 
+(Stgn— 4) e™ |. (19-05) 


Here we used the new designations 


\ 
@=4,, 8% 


t= 1, =p(1—sine). (19.06) 


Equations (19.02) and (19.05) are simplified if kay >> 1 and 
ka, >> 1. Thus, the field from the spherical ring will equal 


a hl, , eft 
Epa — Hy = —Eu($ 18 — Fe mS, (19.07) 

and the field from the spherical segment will equal 

aray elt 
Ey =—H,=:— “22 (seca—tg oe) — (19.98) 
If here one assumes w # 0, then equation 

aF,, ef .69 
E,=— Ay FR e109) 


gives us the field scattered by a hemisphere. The value of the 
effective scattering area corresponding to it will equal, in 
accordance with (17.06), 


of = na’, (19.10) 


PTD-HC=23-259-71 109 


Now let us find the field scattered by the spherical segment 
considering the discontinuity of the surface; one mav neglect the 
perturbation of the current. as a consequence of the smooth curve of 
the surface 1f ka >> 1 [74]. The nonuniform part of the current 
which 1s caused by the discontinuity creates in the direction #=a 
the field (17.03). Summing the latter with the field (19.08), we 
find the desired field 


2 2 
—= sin— 
eed Coen co (19.11) 
C 


Consequently, the effective scattering area of a spherical segment 
will equal 


— sin — 
i: a a 2iht 
a= a" cos @ 2 Te * : 
cos-,; — cos (19.12) 
= a+2 
a=1+——. 


In the physical optics approach, a similar quantity is determined 
by field (19.08) and equals 


1 zzall* 
# == 2a*| -— —tgwe™'|. (19.13) 


With the deforming of the spherical surface into a disk 
RY . 
(o> F..4+0. Q=; const ) » Equations (19.12) and (19.13) are transformed, 
respectively, to the form 


' 1 
a= ta" | ika + — ctg = 





2 
| (19.14) 
ae | : 

It follows from Equations (19.12) and (19.13) that the effective 
scattering area of a spherical segment is an oscillating function of 
its length. The oscillation period equals 3. Numerical calculations 
performed on the basis of these equations showed (Fisure 43) that, 
with small angles of the discontinuity (2 = 15°), cone may still 
neglect the field from the nonuniform part of the current. In Figure 
44, graphs are constructed for the effective scattering area of a 
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Figure 43. The effective scatter- Figure 44. A comparison of the 





ing area of a spherical segment effective scattering area of 

as a function of its length with a spherical segment (continuous 
a constant radius of the base. line) and a finite cone (dashed 
The function o (the continuous line) which have the same bases. 


line) ts calculated on the basis 

of Equation (19.12) which con- 

siders the nonuniform part of 

the current near the discontin- spherical segment and a finite 


uity. The function o° (the 
dashed line) is calculated from cone (the dashed curve) which 


Equation (19.13), and corres- have the same diameter and tase 
ponds to the physical optics shane 
approach. Res 


The results obtained in this Chapter show that the reflected 
signal depends suostantially on the shape of the shaded part of the 
body, and increases with an increase of the concavity. However, 
since the nonuniform part of the current is concentrated mainly near 
the discontinuity, that part of the shaded surface which is several 
wavelengths away from the discontinuity evidently will not have a 
noticeable effect on the reflected signal and may be an arbitrary 
shape. 
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It is interesting that our expressions, which agree satisfactor- 
ily with experiments, even with large (in comparison with the wive- 
lengths) dimensions of the bodies, do not change into the physical 
optics equations, but differ from them substantially. At the same 
time, physical optics, contrary to the widely held opinion concerning 
its reliability in such cases, leads to a significant discrepancy with 
experiments. 


The method used in tis Chapter allows one to calculate the 
effective scattering area associated with the symmetric irradiation 
of any convex body of rotation, the surface of which has circular 
discontinuities. It may also be generalized to the case of asymmetric 
irradiation. However, when doing this it is necessary to take into 
account the nonuniform part of the current caused by the point and 
the smooth curve of the surface. 
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FOOTNOTES 


1. on page 98. See footnote on page 86. 
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CHAPTER V 


SECONDARY DIFFRACTION 


In the previous chapters, an approximate solution of diffraction 
problems was carried out which was based on the representation of the 
fringing field in the form of the sum of the fields from the uniform 
and nonuniform parts of the surface current. The first field was 
found by quadratures, and the second field by approximation; it was 
assumed that the nonuniform part of the current near the discontinuity 
(edge) of a surface is the same as on a corresponding wedge. 


However, the fields found by such a method are actually the 
fields from the currents flowing, not only on the flat and curved 
parts of the body's surface, but also to some extent on the geometric 
extension of these sections. The error in the expressions for the 
fringing field which is thus introduced is most significant with a 
glancing incident wave, when the edge zone occupied by the nonuniform 
part of the current is noticeably broadened, and also with a glancing 
radiation, when the direction to the observation point forms a small 
angle with the given section of the surface. In these cases, the 
results obtained earlier are in need of substantial corrections. We 
already talked about this briefly in § 6 and § 12. 
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Por the purpose of refining the solutions which were found 
previously, it is necessary to assume that in actuality the currents 
flow only on the body's surface, and that a wave travelling from one 
edge to the other will undergo a perturbation at the latter. The pro- 
cess of forming the fringing field when this occurs may be investigated 
in the following way. The edge wave propagated from one of the edges 
is diffracted by the other edges; the waves arising with this in turn 
are diffracted by adjacent edges, etc. In this chapter, we will 
investigate the case when the dimensions of the surface faces are so 
large in comparison with the wavelength that it is sufficient to limit 
oneself to considering the diffraction of only the primary edge waves. 
This phenomenon we shall call secondary diffraction. 


In this chapter, secondary diffraction by an infinitely long 
strip (§ 20 - § 23) and by a circular disk (§ 24) is studied. The 
solution of these problems may be obtained by means of the principle 
of duality from the solution of the diffraction problems for an 
infinite slit and a circular hole ina flat, ideally conducting screen. 
In the latter case, the physical treatment of diffraction of edge 
waves is significantly simpler; it 1s exactly for this reason, therefore, 
that almost all diffraction studies of edge waves are related to holes 
in a plane screen. However, we will not take such a path, but we 
shall investigate a strip and a disk directly. This approach has the 
advantage that it is easily generalized to the case of three- 
dimensional bodies. 


§ 20. Secondary Diffraction by a Strip. 
Formulation of the Problem. 


Let an infinitely thin, ideally conducting strip of width 2a and 
unlimited length be orientated in space as shown in Figure 45. A 
plane electromagnetic wave incident normal to the strip'’s edges is 
directed at an angle a to the plane xoz and has the following form: 


E=E,e* rms? velar H=H,e*" tose + wala) (20. 01) 
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In § 6 approximation expres- 
Sions were found for the fringing 
field in the far zone which did 
not consider the interaction of 
the edges. In the case of E- 
polarization of the ineident wave 
(E,l| 02) , these expressions May be 
represented in the form 


Erm — H, = Eqs lf (yet ere sinsy 








(aes) 
FAR) erBoee sens : 
Figure 45. The transverse cross = 
section of a strip with the £,=H,=0, 
plane xoy, x = 0, y = a and 
x = 0, y = -a are the coordi- (20.02) 
nates of the strip's edge; n is 
the normal to the incident plane 
wave front. and in the case of H-polarization 
(H, (02) 
He= EB, = Hyslg (1) el ine 8) 
‘ (era) 
— dee (sine — 4 : 
Fate meme mny ce, 
4, = &,==0. (20.03) 


Let us recall that the functions f and € included here are determined 
in the region (elegy (when lai<c>) by the following relationships: 








cos “FT stg SP . cos thy sin 
Mh= ine— ang 1Q)= - —ay , (20.04) 
a(i)=—f (2), g(2)= —] (I). (20.05) 


The first terms in Equations (20.02) and (20.03) describe cylin- 
drical waves diverging from edge 1 (y = a), and the second terms 
describe the cylindrical waves diverging from edge 2 (y = -a). The 
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nonuniform part of the current on each side of the strip also has the 
form of waves which diverge from edges 1 and 2, and are an “analytical 
extension" of the corresponding terms in Equations (20.02) and (20.03). 
The current wave encountering the opposite edge is reflected from it. 
Or else one may say that each of the cylindrical waves propagated from 


edge 1 or 2 undergoes diffraction by the opposite edge (secondary 
diffraction). 


If the strip's width 1s sufficiently large in comparison with 
the wavelength, then one may appr«ximately assume that the oncoming 
current wave near the strip's edgz will be the same as on a corre~ 
sponding half-plane excited by a linear source, the moment of which 
is selected in a definite way. It is also obvious that the current 
waves reflected from the edge will also coincide. Consequently, the 
problem of secondary diffraction by a strip may be reduced to the 
problem of the diffraction of a cylindrical wave by a half-plane. 

: The field created at the point P by a current filament parallel 
to the half-plane's edge and passing through the point Q (Figure 46) 
may be found by means of the reciprocity principle. In the case of 
E-polarization, it 1s determined by the relationship 


E,= ££, (Q), (20.06) 


and in the case of H-polarization 


ms 
He = a 4s (Q). (20.07) 
Here P, (m,) is the electric (magnetic) moment of the current fila- 
ment passing through the Q; Pos (m2) is the moment of the auxiliary 
current filament passing through the point P with the coordinates 


(o", R), and H, (Q) or E. (Q) is the field created by the auxiliary 
filament at the point Q. 


Now let us remove the auxiliary current filament to such a 
distance that the cylindrical wave arriving from it may be considered 
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to be a plane wave on the section 
from the edge of the half-plane 
to the point Q. In this case, in 
accordance with § 1 and § 2 the 
field created by it at the point 
Q will equal 


E,(Q)= E,1(0) [a (a, ¥ —¥") 4d. +7 )) 
Hz (Q) = Hye (0)[u(d. % — 9") bald, 779") 





(20.08) 


Figure 46. Diffraction of a 


cylindrical wave by a half-plane. The functions u introduced here 
Q is the source, Q* is the 
mirror image or the source, are determined (for the values. 


and P is the observation point. One’ <t) by the equations 


ae Vind eos t= 
u(d, & — 97) se Nt cmt) as x eltda-y 
s ; a cos woe’ 
+ eer HO Hee 
0 with +9 <7 <2, 
° ; 5 Vina cos EAE 
ald, 7 Hype ee rd x e'“dq 4+- 
weaey? 


+{ es cooty ee enO FE e” 
0 withr—— <P <2z, 


(20.09) 


and the quantities E, (0) and H, (0) are the values of the primary 
field created by the auxiliary filament at points corresponding to 

the half-plane's edge. In accordance with Equations (1.21) and (1.22), 
this field may be represented when kR >> 1 in the form 


i faa v) 


Ay, (O) = ik? m,: (20.10) 


ER 


FTD-HC-23~259-71 117 


Consequently, an electric current filament located above an 
ideally conducting half-plane excites, at the point F, the field 


e 
Roy 


Ox ga — 
E,=ik'p,lutd. ¢ --2")--ald eg tery V me ) (20.11) © 


and a magnetic current filament excites, at the point FP, the field 


7 (ae i) 


mee 2 
Hy =ikim,\u(d 9-9) La EP ee (20.12) 


-It 18 easy to see that the exponent e*"-4e*¥ =" 4n these expressions 
corresponds to the primary cylindrical wave arriving at the observa- 
tion point P, and the cxponent ¢#l&-4°" 49" +e"l corresponds to the reflec- 
ted cylindrical wave. 


The moments Mm and P. must be selected in such a way that in 
the dlrection 9" 2 m (Figure 4A) the filament would create a field 
equal to the field of the primary edze wave above an infinite, ideally 
conducting plane. We will conclude tnese calculations in the follow- 
ing sections, but for now let us make still one other ccmment on une 
formulation of the problem. 


In the previous chapters it was shown that the scattering object 
may be approximated by a series cf sources — "luminous" lines and 
points. Therefore, the problem of seccndary diffraction may be formu- 
lated as a rrotlem of searchin for functions which describe the 
continuous change of the field of each such scurce juring tl.e passage 
through the boundary of the lisht and shadcw torresponding to the 
source. 


§ 21. Seconuary Diffraction by a Strip (H-Felarization) 





A current filament with the moment m_ which *s positioned abcve 
an ideally conducting plane (h = 0, Figure 46) creates in space the 
field 


PUD=HC =? 3-25-71 ilf 


My - ek am 2H (RR). (21.01) 


Far from the flbament (when KR, >> LP), this ffleddo fa deserthed by the 


Aoympbot le oxpresston 


(ems <) 
Hos Azk'm, £.. ib ieaees Cob.) 
V 24k, 
Rut the peimary edye wave dn the direction @" = mn takes the value 


(are +) 


e 
HW, ie #40, (Q) 8 (Q) 'Y tneR, ) (2.0 3) 


where Hi. G2) Ia the Cfeld of the fnetdemt plane wave at the potut vs 
rma) ts the value of the angular Puneet bom of the primary evbindrefecdt 
wive In the direction towards the oeppostte edge of he steip.  Bqust- 
Io Expressdtoms (Ob.00) and (2beoo), we Cind the ¢(blament ts moment , 
the Pleld of whieh we use te approx timate the pieimary ede wave, tu 


the form 
t 
Ms anys er (QE (Q). (21.08) 
Aso a result, the tMledd erented by the Chlaments: beeated ohoave 


the balfeplane -@ayn and correspond iiys te ede 1 Ce ture yous 


be represented tor royton fpps r fro the form 
HAs AES). eee) 


The Cunet fon 


Wise (- 2) 


iar 
1 = ' ® ttaqina: ain 
Hiss 6 Mag (Qs f eC dy yage iv at air 


ae: 
deseribes the wave padiated by the seaureé om mod thre Mune t bor 


| ea 


PPD HC et sane I [he 


WV hacos (+-4) : 
4, ! = i i? ale fa (elm emote ot 
H7()= +Hug (1) X | eda i ae 


(#03) 


+Auall) Ee eltaitine-siner (21.07) 


describes the wave radiated by 
the source m2 The sum of these 





Figure 47. The problem of waves equals 
secondary diffraction by a 
strip. # ince 3-4) 


+ - 
m,. and m,_ are the sources 2 - : yy, em 
1z 1z > ALA)=— Agi el dg ~ el tite sin' 
the fields of which are used (rg Mat OX \ ” ¥ 27° + 
when approximating the primary os 
edge wave being propagated from ae lees 
edge 1 (y = a); tHe. 
: dm, are the sources ? 
m5, and m5, a F 


the fields from which are used 
when approximating the primary 
wave from edge 2 (y # <a). (21.08) 


itu grins -sin ED 


e 


The first term in this expression is the desired secondary wave from 
edge 2, and the second term is the field radiated by the filament 
which its loacted above the ideally conducting plane x ='0 and has the 
moment 


saa sina 


Saar = 
Me seat Hag (ie 


(21.09) 


where 


z= eth 


. ou 
¢a> - 


(21.10) 
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Summing the secondary wave which has been found with the unper- 


turbed primary wave from edge 1; we obtain 


H,(\—2)= 
Wiacos (7-4) 
- iav ; 
=H,, Fax f ed mae 


r 
i( te » = 


the (ein o— si8 ¢) 


e 
1b Aae Va (21.11) 


This expression reduces to zero 1f one assumes $ = = 1/23; consequently, 
the secondary diffraction eliminates the field discontinuities which 
occurred in the previous approximation when ¢ = - 7/2. However, in 
the direction » # 1/2 the field (21.11) is different from zero. Since 
H, is an odd function of the x coordinate, the relationship 

t!,4.2_,%9 means that the fringing field components H, and Eg will 
undergo a discontinuity with a transition through the direction 
¢ = 1/2. The reason for such a jump, as before, is that in our calcu- 
lations the plane x = 0 is a plane of currents. By finding the 
secondary wave from edge 2, we actually considered that the diffrac- 
tion takes place not on the edge of a finite width strip, but on the 


edge of an ideally conducting half-plane —aqcyqm. 


Again the resulting discontinuity has an order of magnitude of 
vEve . It is clear that one may completely eliminate the field 
discontinuities only with consideration of multiple diffraction. 
However, the calculation of fields arising with multiple diffraction 
requires specific consideration of the following terms in order of 
smallness in the expansion of the primary edge in inverse powers of 
ir (see, for example, [46]). All this greatly complicates the cal- 
qilations. Therefore, we, using the condition ka >> 1, will limit 
ourselves to an investigation of secondary diffraction, and in order 
to eliminate the discontinuities in the plane x = 0, we will proceed 
in the following way. 
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Let us consider the quantity g(1) in the Expression (21.11) to 
be a function of the angle ¢ [see Equation (20.05)], that is, let us 
replace @(2) by the function g(1). In this case the equation 


va ca(T-F) 5 
=e: ~ git) x \ e'"'dq wet 
_ Aa 
+ Hag (1) eineaine) : fab ads 


will give qualitatively correct results not only when PF» but 
also with all other values of ¢. Actually, the Fresnel integral is 
Glose to zero if Viacos (-*--+)>! » and in Equation (21.12) only the 
second term remains, as must be the case. Therefore, Equation (21.12) 
may be investigated as an interpolation equation, and it may be 
applied with any values of a(iPi<a,. It is easy to establish that 
now the fringing field does not undergo a discontinuity with the 


passage through plane x # 0, since Expression (21.12) becomes zero 
when g= . 


It is interesting to note that Equation (21.12) automatically 
follows from Equation (21.08) if in the latter equation one replaces 
g(1) by g(1). Essentially, this substitution is equivalent to the 
assumption that the moments of the filaments, the fields of which are 
used for approximating the primary edge waves, depend on the radiation 
direction (that is, on the azimuth 4 of the observation point) 


m= — mt Hag (eee (21.13) 


Such a determination of the moments of the auxiliary linear sources 
is used, for example, in the work of Millar [47]. 


Precisely in the same way that Equations (21.06) and (21.07) 
were obtained, we find (when x > 0) 


FTD-HC-23-259-71 122 


Wa eos (++3) 


Ht 2)-—= | H 2 : itdg ne wits isine—sin sy 
P= Muei2) | ela ee 


‘= (21.14) 
PRPS EST) 
Hy. (2) = + — Hu 2 ) old ee sem (ta lrtnen inet 
ieee (21.15) 





+ Hug (2)° Vis e~ ite (sin o—sin of 


These expressions give the field created by the filaments which are 
located above the ideally conducting half-plane —«<y<e and have 
the moments 


mo = -- mt =o Hag (2)e~mn*, (21.16) 


In accordance with Equation (21.04), here 


-0 


eae! (21.17) 


Furthermore, summing (21.14) and (21.15), we obtain 


dab a z 
ear 
Hy (2)= Hee g (2) \ e“dq yee 
i(ar pl 





= e mika iain e—aine) 
+ Hae heer ¢ (21.18) 


Here the first term is the desired secondary wave from edge 1, and 
the second term is the field radiated by the filament which is located 
above the ideally conducting Plane x = 0 and has the moment My 5» 
Summing the secondary wave which has been found with the unperturbed 
primary wave from edge 2, we have 
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Wihacos (4+4) 

H,(2—NHN=H, 2 R12) \ 
ea 
ifarte é 


4 
+ Hag (2)* wie abe amie (21.19) 


ier ; 
I © ..-ihe (stn o—sin¢ 
e'*dq rin + 





It is not difficult to see that the resulting expression becomes zero 
if one assumes 9= 5 in it. Consequently, the secondary diffraction 
eliminates the field discontinuity which we had earlier (§ 6) when 
= +; but at the same time it leads to a field discontinuity when 
e=— 5 Again the resulting field discontinuity may be eliminated 
by the above indicated method, replacing the quantity £(2) by @(2) 
— that is, by assuming the moments m5, and n>, depend on the obser- 
vation angle ¢. Actually as a result of such a substitution, we 
obtain from (21.19) the expression 


eltr 
H,(2-)=Ha 2 g(2) \ | eda meen 


e « J 
t (e+ 


+HusQ oa e 


2V' Fe cos (+4) 


—tha (sin ashe ¢) 
* 


(21.20) 


which vanishes when gay » This expression may be investigated as 
an interpolation equation which deseribes the field created in the 


region lel<> by the primary wave of edge 2 with consideration of 
its diffraction at edge 1. 


Now summing (21.12) and (21.20), we obtain the following expres- 
sion for the total field scattered by the strip: 


H,= He lG(i, pg (iecine—sin 


1 (a + 7) 
+G(2. pg (ayer eire—) ae : (21.21) 


Here 
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gt avons (3-4) 
Gi =zre | x f edg 
8 


(21.22) 
is the shading function of the primary wave travelling from edge l, 
and 


: Wtacee (743) 
G2y=soe © y { eda 
(21.23) 


is the shading function of the primary wave travelling from edge 2. 
These functions show that the primary wave from edge 1 undergoes the 
greatest perturbation when 92> » and the wave from edge 2 under- 
goes the greatest perturbation when p=>. 


An important property of Equation (21.21) is that it becomes 
zero when =2p — that is, the field discontinuity which we had 
earlier at the plane x # 0 is completely eliminated. 


In coneluding this section, let us return to Expressions (21.11) 
and (21.19) which lead to discontinuities of the fringing field in the 
plane of the strip (x = 0). One may show that the sum cf these 
expressions 
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agrees, when Vkacos(-+-$)>! , with the asymptotic solution obtained 

in the book [50] by means of integral equations. The solution found 

in [50] has the greatest precision when a=0,"9=0 , and it is completely 
useless if amcty or gutsy - 


§ 22. Secondary Diffraction by a Strip (E-Polarization) 

It is known that a current filament with an electric moment Po 
which 1s found at a distance h from an infinite, ideally conducting 
plane (see Figure 46) creates in space the field 


3 = ik"pye [HO"RR,) — HO (kR,)). (22.01) 


With small values of h (and R 
formed to the form 


1.2? kh”), this expression is trans- 
’ 





= i (amet , 
=~ iapwasing V Be! 4) (22.02) 
The primary edge wave is determined by the relationship 
(tet) 
E.=Eul9)! (55° (22.03) 


where E424) is the value of the incident plane wave field at the 

point q (Ry = 0). Consequently, the primary edge wave in the direction 
w = 0 may be investigated as the wave from a current filament located 
above an ideally conducting plane if one assumes the filament moment 

to be equal to 


i f@) 
Pe = Fagg Ee (9) “stag et (22.04) 





The fie'd, created at the point P by the current filament with 
a moment Po which 1s parallel to the half-plane's edge and passes 
through the point Q, is determined by Expression (20.11). Expanding 
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the right-hand member of this expression into a series in terms of the 
emall quantity h(h + 0) and limiting ourselves to the first term which 
is different from zero, we obtain 


E.= itp, 5 lald, ¥—¥)— 


—al(d, ret Vz (t-4). (22.05) 


By means of relatiorships (22.04) and (22.05), one may show that 
current filaments with moments Piz and Br. which are located on the 
ideally conducting half-plane —e@syqoo and correspond to edge 1 (see 
Pigure 47) create in the region lele> the field 


Wieew <—-+ 
fem a "0S 9 j e*dg-+ - 
oe 


+ pag tio (Tap jereetme [payee 
2 ‘ wes 
+ Euf (1) cos 94 ViE moma (22.06) 





The current filaments with the moments P5, and Po. which are located 
above the iceally conducting ualf-plane —eceyqa and correspond to 
edge 2 create in the same region the field 


Vive (54-3) 
Eyer But S cos¢ j edo 


+75 sin (+ + me ae Pt (2) eg tts (ola e—sin m4 


tle 7) 
e = ihe (sia oie ¢) 


+£al 2)cos9 yie ° 





(22.07) 
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The first terms is. Expressions (22.06) and (22.07) ere the desired 
secondary waves, and the last *erms in the expressicns are the fields 
from the current filaments located abcve the ideally conductinz piane 
Zz @ 0 and having the moments 


Pe hig uw NW, pie tg Ew Oe (22.08) 
where 
1 ty, 18 (22.09) 
119): cos tng as' or ar 


Summing the secondary waves which have deen found with the 
unperturbed primary waves, we obtain the total field scattered by the 
strip 


266 evn {z-4) 
Fy xs A, = Ew? Fag {F cuore x j edge 


\e Mae > cin "| inate rm simeh 


tiygt9(4-4) 
2 Vi6 cw (+ +4) 
+72) [cose eda + 
+s (4 +4) ona aaa 
Beal (1) ers MP9 4 7 (2) em stettn ome od |X 
ives) 
xR o (22.10) 


Now assuming, as in the case of the ti-polarization, that the 
moments Pie and Po. depend on the angle ¢, by replacing 


fay vy Sana fizy vy EE. (22.11) 


we obtain 


By — Hm EeclPh, 9h] (ee ey 


+ P(2, p)f(zye rine My ee 
(22.12) 
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where 
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(22.13) 
are the shading functions. They show that the primary wave from edge 
1 undergoes the greatest perturbation near Is » and the wave from 
edge 2 undergoes the greatest perturbation in the vicinity of =s° 


§ 23. The Scattering Characteristics of a 
ae Plane Wave by a Strip 


Expressions (21.21) and (22.12) which were obtained above for 
the field scattered by a strip approximately take into account the 
interaction of the edges and are valid when lel<> . However, they 
are not applicable with a glancing incidence of a plane wave on a 
strip (when a= ). 


In order to find equations which are applicable in this case, 
let us proceed in the following way. Let us write the expressions 
for the field radiated by the strip in the direction a with the inci- 
dence of a plane wave in the direction » (Figure 45) 

E,=—H,= Eu F (2, a) j (1) eftarinn— sted 
| i (net 7) 
a i a 
Ay B= Hus {G (2. a)g (iene ey 
| i {tre <) | 
} 


F e 
+G6(1. alg (2) eo imain a—sin¢) ] jn ‘ 





(23.01) 
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Here lel< 5 » but 9 cannot approximate 71/2. Now let us note that 
the expressions for the fringing field must satisfy the reciprocity 
principle — that is, they must not change with the simultaneous re- 
Placement of a by ¢ and ¢ by a. Comparing Equations (21.21), (20.12) 
and (23.01), it is not difficult to obtain the expressions 


E,=—'H) = El F (1, 9) F (2, a) (reese 
‘ +5) 


+ F(2, GF (1, apf (2ye reine Via : 


Hye Ey = Hel Q(1, 9)G(2, 3) g (1) rns 


(wee) 
+ G(2, 9)G(1, 2) g (2)¢nnime- siory Se : 








(23.02) 


which satisfy the reciprocity principle, have no discontinuities any- 
where, and are suitable for making calculations with any values of a 
and ? aes, lei<f)- From the second equation of (23.02), it 
follows that H, = Ey * Q when 9=t> — that is, the fringing field 
does not experience discontinuities in the plane x = 0. Moreover, 

H, * E, = 0 with any values of ¢ if a=+ — that is, a plane wave 


polarized perpendicularly to the strip does not undergo diffraction 
with a glancing incidence. 


The resulting Equations (23.02) may be investigated as interpo- 
lation equations. Actually, with Jal€ when Vika cos(7.= 3 )>1 the 
functions Fil, 2), F(2, 1) G(l,a) and G (2, a) are close to one, and 
Equations (23.02) change into the previous Expressions (21.21) and 
(22.12). But if ipl<€-$ and Viacos(-=-$)>1, then the functions 
F(l, 9) F(2, 9) , G (1, ¢) and G (2, #) are close to one, and Equation: 
(23.02) change into Equations (23.01). Let us recall that the fun:- 


tions F and G are determined by relationships (21.22), (21.23) and 
(22.13). 


In the direction of the principal maximum of the scattering 
diagram (¢ = a), Equations (23.02) take the following form: 
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= H,,; [ (2ite cosa NG 13) G02, app 

oG(2, a) - 4G(1, a) ‘ agg +) 

FOC, gy E8 9, wy Med yer 
F.= F|(2ik0 cosa... ccs PA, ayF(2. a)-- 


oF (2,3 : OF (1,4 «! : 
PAL 2) SOB (2, ay MMe |S eee (23.03) 


a 
Hence when a=—=>7F we have 


H,=0, 
at th me Bo joe 
ma Fa(2) e [-vR f e +e + 
So ge eit 
aaa tive f oda ie. os ae (23.04) 


It is interesting to observe that Expressions (23.02) to some 
extent take into account, in addition to Secondary diffraction, also 
tertiary diffraction. Actually, for the values I2!€3 and lees, 
we have 


G qa 9) G2, a) eitatsin a~—sin 9) = eldaisin em~sing) 


s 
efte(2+ sin atsin 9) ' 7 efka(2—sin a—sin : 7 ee 
_ e —-———_—________. e 
2 xka cos (4-4) 2 \/ xk cos (+4) 
i eftatt—sin atsing) 
beg (2-4) F +4) ‘ (23.05) 


The physical meaning of the four terms in the right-hand member of 
this equation is illustrated in Figure 48 (Figure 48a corresponds to 
the first term; Figure 48pb corresponds to the second term, etc.). 


Taking into account condition (6.15), one may write the eauations 


for the fringing field in the left half-space ($<le<< but fala ) 
in the same form as (23.02). Thus, the functions G(1, a), G(2, a) and 
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‘al by a? 


Figure 48, The schematic diagram of the 


waves corresponding to the various terms 
in Equation (22.05), 


F(l, a), F(2, 2) will, as before, be described by the relationships 


(21.22), (21.23) and (22.13). ae remaining functions in Equat‘ons 


(22.02) will be determined when - © <lel<e by the following equations: 


i a2¥8e sin (+ +#) 
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Calculations of the scattering characteristics were carried out 
based on the equations derived above. These scattering characteris- 


tics are the functions h(a, 9) and e(a, ¢) determining the fringing 
field by means of the relationships 


i 
é,= E,,kae (a, of a e . 7), 


_ i(wo%) 
Hea Habana. yf Bes". (23.13) 
The calculations were performed for the values kKa=Y28 and f£a=V80 - 
In Figures 49 ~ 62, the following designations were used: 1) the 
functions h and e correspond to the rigorous theory; 2) the functions 
hy and eg correspond to the field from the uniform part of the current 
(the physical optics approach); 3) the functions h, and e, correspond 
to the field from the uniform and nonuniform parts of the current, 
but without consideration of the interaction of the edges; 4) the 
functions hy and e5 correspond to the fringing field with considera- 
tion of secondary diffraction calculated on the basis of equations 
(23.13), (23.02) and (23.10). Thus, in accordance with § 6, 


con ge tlnilta (sto s — sia py > 
eo CON 2 Paisins — sin 9) ’ 


h,=cos gt! a [ka (sla a — ao 9) 


ka(sins — sin») (23.14) 





and 


& : slofka(sin 2—sin 2. joes {ha (sina — ein »)] 
"tka 


h, int cost (23.15) 
where ales. lel<z : 


The results obtained show that our approximation equations arree 
satisfactorily with the rigorous theory already when tu~-| 28 , 
although in the given case approximately one and one-half wavelengths 
are fitted into the width of the strip. In the direction toward the 
source (e=—e+2. Ccac 5 » and also with glancing irradiation of the 
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The 


Figure 49, 


Scattering diagram of a field by a strip as a function of 


The various 


the incident angle of a plane wave (a) and the parameter vka. 


curves correspond to various approximations. 
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he same as Figure 49 with «= — 
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Figure 51. The same as Figure 49 when 
@? = -7 +a, 


strip (a = =-1/2), when the functions €, and Nos and Qe, and hy lead to 
qualitatively incorrect results, the functions e5 and h, give, as in 
the remaining cases, fully satisfactory results. Actually, the curve 
|h,| coincides almost everywhere with the curve |h| (Figure 49-54) 
within the limits of graphical precision. But the calculated values 
of the function le] differ from the corresponding values of the 
function Je] only by hundredths of a percent (Figure 55 - 62). The 
better agreement with the rigorous theory associated with the E- 
polarization is explained by the weaker interaction of the edges in 
this case. A certain discrepancy of the curves [hs] and |h| in the 
vicinity of the principal scattering maximum is explained by the 
interpolation character of our equations. 


As a consequence of the interpolation character of Equations 
(23.02), the integral scattering diameter obtained from Expressions 
(23.03) when a = 0 does not coincide with the integral diameter found 
by Clemmow [46] in the form of the first terms of an asymptotic ex- 
pansion in inverse powers of “ka. However, our equations,as disttnet 
from the similar equations obtained by other authors, allow one to 
calculate the scattering characteristics with any Incident anrles of 
the plane wave. 
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The function h (a, 6) for a strip (a = 0, ka = 





0). 


Figure 52. 
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The same as Figure 52 when a = -7, 
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Figure 53. 
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Figure 54. The same as Figure 52 when ¢ ® -7 + a. 


Let us note that the functions e(u, ¢) and h(a, 9) for Figures 
49 =~ 52 were calculated on the basis of rigorous series which were 
obtained by the separation of variables in the elliptic coordinate 
system (compare (237) ')), 


4 


§ 24. Secondary Diffraction by a Disk 


Let us refine the approximate solution of the diffraction pro- 
blem for a disk which was found in Chapter II. 


Let an infinitely thin, ideally conducting disk of radius a be 
found in free space. Let us orientate the spherical coordinate system 
in such a way that the normal n to the incident wave front would lie 
in the half-plane » = 1/2, and form an angle 1 (0<1<4) with the z 
axis (Figure 63). Let us prescribe the incident plane wave field in 


EE, efMesinrtscotn) H=H,e*eset +econ qT) (24, 01) 


1) pootnote appears on page 162. 
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The same as Figure 55 when a = -n, 


Figure 57. 
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Figure 62. The same as Figure 59 when $ = -n + a. 


In accordance with § 12, the fringing field in the plane 
6 = tn/2 is described (when R >> ka*) by the equations 


== (02.9) — 1.0 -+} 
r eink 
FIM O10. M40} 


Hy Ey age (et. ) — (1, A+ 


a elt 
+ ilgl2, Y+g (lL. MAO ~~ - (24.02) 


These expressions are valid when <> and 1<F - The quantities 
included in them are determined by the relationships: 





cos oe 
F(t, y= ain@é—siad = ° 
cos tty sta = 
12. \=—— 


sind—sind ° 


(24.03) 
C=: ka(sin 8 — sin 8), (24.04) 
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|-1 with 9=—>- (24.05) 
Let us note that here 
FQ) with e=> 
f(1.3y= “ g? 
f(2) withg=— > 
(2) with =~ 
f(2.4)= «’ 
[F() with e=—> 
e(!) witne=> 
a(i, y= a2’ 
1 (2)with => 
fee with p= 
a (2, %?) = a2’ 
[eQ)withe= ~—= 


(24.06) 


and the functions f(1}, f(2), g(1) and g(2) are determined by the 
Equations (12.03) and (12.04). 


When ¢ >> 1, Expressions (24.02) take the form 





ink, (3) 
E,=—H,= Fiz |! Ae Whee 





—ae(lsje eT: (24.07) 


They show that the fringing field in this region may be investigated 
as the sum of spherical waves from two luminous points on the rim of 
the disk with the polar angle » = +7/2. The diffraction by a disk of 
each of these waves may be studied as was done in the case of a strip, 
but we shall proceed differently. 
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Figure 63. The cross section of pipure 64. 
a disk with the plane yoz; n 
is the normal to the incident 
wave front. 


Excitation of a half- 
plane by an elementary dipole 
which is located at the point Q. 


Let us compare the shading of spherical and cylindrical waves 
by a half-plane. Let an ideally conducting half-plane be found in 
free space, and let there be an elementary dipole at the point Q 
(Figure 64). Let us find the field in the plane perpendicular to 
the half-plane'’s edge and passing through the point Q. 


In accordance with the reciprocity principle, it 1s determined 
for the electric dipole by the relationship 


E.="* E:(Q) (24.08) 
and for the magnetic dipole by the relationship 
Hy = 5" HQ). (24,09) 
Here p,(m,) is the electric (magnetic) dipole moment found at the 


point Q; Poz and My, are the moments of the auxiliary dipoles which 


are placed at the point P; E,(@) and H,(Q) are the fields created by 
the auxiliary dipoles at the point Q. 
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Now let us remove the auxiliary dipoles to such a distance that 
the spherical wave arriving from them may be considered to be a plane 
wave on the section from the half-plane's edge to the point Q. In 
this case, in accordance with Equation (20.08), the field created by 
the wave at the point Q will equal . 


E,(Q) = E..(0) (a (4,97 — ¥")— u(d,7 +9"), 


; Hs (Q) = Hus (0) (u(d, & —9") + u(d, +9") (24.10) 
The expressions 
UF 
Fos(0) Des, Hye (0) = ktm (24.12) 


determine the fields created by the auxiliary dipoles in free space 
(with the absence of the half-plane) at the point 0. 


Consequently, the fields excited at the point P by the electric 
and magnetic dipoles which are found at the point Q above the half- 
plane equal respectively 


E.= lp, [u(d,9 —9")—uld¢ tee, | 
Hes hime [u (9 —9") ade tery. J (24.12) 


With the absence of the half-plane, these dipoles create at the point 
P the field 


eine 


E,= rey : Re en itd cos en 


Comparing Expressions (24.12) and (24.13) we find the shading 
functions 


F=[a(d,9'--3")—u (d, 9” 4 97)) ltd 809 ("21") 
64 [ 


u (d, 9 ~ 9") Lu (d, 9’ + 9")] elt eonter—ven, (24.14) 
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In the case wher the current filament passes through the point 
Q parallel to the half-plane's edge, the field at the point P is 


determined — in accordance with (20.11) and (20.12) — by the 
equations 


aa i{en— 
E,=ik'p,|u(d.¥ — 9") —u(¢ +e Ve | I 

—- t(ee—* 
Heit (ald — 91-40 +9 ee ( ‘| (24.15) 


With the absence of a half-plane, these sources create at the point 
P the field 


p=ih'p, Ve a(re— +) en itd cos (9-9) : 


H,=ik'm, Vz . ("-3) ere ee), (24.16) 


Comparing Equations (24.15) and (24.16), we obtain the same Expressions 
(24.14) for the shading functions. Consequently, a spherical wave in 
the direction perpendicular to an ideally conducting half-plane is 
shaded by it the same as a cylindrical wave. 


Let us note, however, that Expressions (24.14) arc not equivalent 
to Expressions (21.22), (21.23) and (22.13), since the first represent 
the shading function by a half-plane of a wave from a sirgle source, 
and the latter represent the shading function of an edge wave which 
we approximate by waves from two sources located on both sides of the 
corresponding half-plane. Since the shading functions of spherical 
and cylindrical waves are the same, the edge wave shading functicns 
of a strip and a disk also will coincide. 


Therefore, the approximation expressions for a field scattered 
by a disk which take account of secondary diffraction may be repre- 
sented in the region Pact. 0<d< > (with ¢ >> 1) in the following 
form: 
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Starting from Expressicns (24.18) and (24.19), it is not diffi-~ 
cult to write interpolation equations for the fringing field which 
are suitable for any values of y and ® in the interval (0; 172), 
but when ¢ = 2/2; 


ey Se Ll 2,8) F 1,2) 2,8) — 
— F (1.9) F (252) f(1,2)] J) PALF (2,0) F (1,2) / (2,804 
HF (1M F(2,2)/(1,3)4 0) (24.20) 


E,=H,="5 [ia2. 9) G(1,2)¢2,3) — 
— G1, 8) (23) a(1,2)] J, (6) +116 (2,0) G(1,2) g(2,8)-+ 
$G(1,9) 62,2) 201,214, 6} 9 - (24.21) 


Let us note that when y = 0 these expressions will be valid for 
any values of the azimuth ¢ , since then any point of space may be 
considered to be located in the incident plane. \ 


In the direction of the scattering diagram's principal maximum 
— that is, when %=1,9= — the fringing field (24.20) and (24.21) 
takes the form 


E,=—H, at By FY Fs nie ent| 


=H, = tte, -G(2,7)G(1, Dg C8 T- (24.22) 


However, these expressions have an interpolation character,and with 
small values of the angle y it is impossible to consider them to be 
more precise than the simple equations of § 9 and § 12. In particular, 
with y = 0, when the fringing field must not depend on the incident 
wave polarization, they give values which are different for the E- 
polarization and H-polarization by small quantities of the order of 
Ts . Therefore, in this case (when y = 0) it makes sense to use 
Expressions (24.20) and (24.21) only far from the z axis, switching 

to Equations (24.02) near the z axis. 
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Equations (24.20) and (24.21) have the following important pro- 
perties, They do not have discontinuities, they include the case of 
glancing incidence of a plane wave, and they satisfy the reciprocity 
principle. From them it follows that £,=H,=0 when $=+ — that 
is, the fringing field does not experience a discontinuity on the 
plane z = 0. Moreover, E,=H,=0 with any values of & if 1%=-5> 
— that is, a plane wave polarized perpendicular to the disk's plane 
does not experience diffraction with glancing irradiation of the 
disk. 


As in the case of diffraction by a strip, the new approximation 
expressions consider to some extent tertiary diffraction [see 
Equations (23.05) and Figure 48). 


{ 
Using Condition (9.04), it is not difficult to write equations 
*® 
for the fringing field in the left half-space (7 <o<5 e= +) 


Ey Hy PF (2, 20) F (168) (2.8) — 


— F(1,s— 9) F(2, 8) 7 (1,8) 4, (0) + 
bi LF (2, x —9) F(1,8)/ (2,8) 


HF (Le) FRA LMNAO (24.23) 
Ey=H,= “7? {le@.2—90(1,2)202,8)- 
— G(1, 2 — 8)G (2, 8) 2 (1, 2)J J, (C)+ 
+ 11G (2,=—9)G(1,2)¢(2,8) + 
4G (1,2—9)G(2,8) 0(1.21,.}, (24.23) 


where the functions f and g are determined by the equations 


g(1,8)=/ (2, eg r 





sind—smnet \ 
cos a —sin a=) i 
((1.)=g2,3—=——2 ___? | 
nae , (24.24) 
In the direction towards the source (a=2—7, e=—-F) , the 


fringing field equals 
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lak, ; 
E,=— Wy= 32 0.9/2.) — 


— FH (2,8) f (1, 20] A) EA LARC, 8) f (2,8) 
+ FP(2,8) F(1, 3,0} 

Ey= Hy 4? fic*(1, 2) (2,8) — 

— G* (2,8) g (1, 8yf 4, (6) + 6 (G7(1.*) g (2,8) 


+6" (2,2) 2(1, 1,0} 5 
2) a(L NA Ole. | (24.25) 





where 





t+siny, 
11,8) = 8 (2,8)= —- Et : 





ch _ '—siny 
112, 8) = g (1,8) = Yang’ | (24.26) 


As was already noted, it makes sense to use Equations (24.25) 
only far from the z axis, changing to Expression (12.15) of the pre- 
vious approximation in the vicinity of the z axis. A calculation of 
functions E(e—y) and S(e—7) (Pigures 65 and 66) which determine 
the effective scattering surface [see Expressions (12.17)] was per- 
formed on the basis of these equations when ka = 5. A comparison was 
carried out of this calculation with the results of measurements. 

The two experimental diagrams (the dashed lines) (2) depicted in 
Figure 65 characterize the experimental precision. As distinct from 
the previous approximations, which lead in this case to qualitatively 
incorrect results [see Equations (10.06), (10.07) and (12.15)], we 
observe a satisfactory agreement of theory with experiment. 


For verifying the results obtained, a calculation was also 
carried ou* of the functions vdc9 ) and v2) (4 ) [see Equations 
(9.07)] when ka = 5 (Figure 67 and 68) with normal irradiation of a 
disk by a plane wave. Curve 1 corresponds to the field calculated 
from the rigorous theory [34]; curve 2 corresponds to the field from 


(2) 


Footnote appears on page 162. 
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Figure 65. The diagram of a disk's effective 
scattering surface when the plane wave's 
Magnetic vector is perpendicular to the inci- 
dent plane. 
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Figure 66. The calculated diagram of a 
disk's effective scattering surface when 
the plane wave's electric vector is per- 
pendicular to the incident plane. 


the uniform part of the current (the physical optics approach). Curve 
3 corresponds to the field from the uniform and nonuniform parts of 
the current, but without the interaction of the edges. Curve 4 corre- 
sponds to the field with consideration of secondary diffraction. 43 
is seen from these graphs, consideration of the edge interaction re- 
fines the previous approximation and ensures better agreement with 

the rigorous theory results. 


The problem of secondary diffraction by a cylinder may be solved 
by a similar method. However, considering that the corrections which 
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(1) 
Figure 67. The function Vigy for a 
disk with normal incidence of a plane 
wave (curve 4). Curves 1, 2 and 3 from 
Figure 20 are drawn for comparison. 


depend on the secondary diffraction here are small (on the order of 
1 dB) when ka = 7, kl = 107, and the equations are substantially more 
complicated, we shall not cite them here. 


In the problems investigated above, the edge waves have the 
character of cylindrical or spherical waves — that is, they decrease 
rather rapidly with the distance from the edge. Therefore, in the 
case when the linear dimensions of the faces are approximately two 
wavelengths, it is sufficient to limit ourselves to a consideration 


FTD-HC=-23-259-71 153 


\ SERED 4e Se eee ee 
Pa ete NT 

See aes aee eee 
Bi (Sel aes ee ie 
2) SE Zee) Bee) ses 
Si BRP ee Ae! 


BEls//MRESEREREA\ Eee 
BEE) ZEEE 
tT 
Seni VERE 

vrWwititlit li wit 





CL IONS 


(2) 
Figure 68. The function Vy for a 
disk with the normal incidence of a 
plane wave (curve 4). Curves 1, 2 and 
3 from Figure 21 are drawn for 
comparison 


of only secondary waves. In Chapter VII we will investigate the 
problem of a dipole in which the edge waves decrease so slowly that 
it is necessary to consider multiple diffraction. 


§ 25. A Brief Review of the Literature 


In this and previous chapters, approximation expressions were 
obtained for the scattering characteristics of a plane wave by various 
bodies. These expressions were derived with the help of physical 
considerations which do not pretend to be mathematically rigorous, 
and they are adequate for sufficiently short waves. In the literature, 
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there are a number of works in which similar results were obtained. 

A majority of these works also are not characterized by mathematical 
rigor, and they are based on certain physical assumptions. Therefore, 
one may relate them to the physical theory of diffraction. Only in 

a few works (related to the simpler diffraction problems) did they 
succeed in obtaining specific results at a higher level of mathematical 
rigor ~~ more precisely, while developing asymptotic methods of mathe- 
matical diffraction theory. 


We will briefly list the most important results obtained in a 
number of papers and books, grouping the material in the following 
sequence: 


1. Diffraction by plane, infinitely thin plates (an infinite 
strip, a circular disk) and diffraction by auxiliary apertures in a 
flat screen (an infinite slit, a circular hole). 


2. Diffraction by three-dimensional bodies with edges (a finite 
cylinder, a finite cone, etc.). 


3. Other diffraction problems. 


When investigating the first group of diffraction problems, it 
is necessary to keep in mind the principle of duality [4] which enables 
one to easily change from a strip to a slit, from a disk to a circular 
hole, etc. In the literature as a rule, they preferred to investigate 
apertures in an infinite flat screen, whereas in our book, diffraction 
by a strip and a disk was studied. This approach facilitates the 
transition to three-dimensional bodies (see the remarks at the 
beginning of this chapter). 


Based on the time of appearance (if we do not consider the works 
of Schwarzschild [25] which we talked about in the Introduction), one 
should first of all mention the works of Braunbek [28 - 30] which 
were devoted to the diffraction of a scalar wave by a circular hole 
in a flat screen. Assuming that the plane wave is incident normal to 
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the screen, the author obtained an approximation solution in the form 
of a surface integral. The boundary values of the integrand were 
taken from the rigorous solution to the problem of diffraction by a 
half-plane which was found by Sommerfeld. The field was calculated 
in the far zone on the axis of the hole and far from it, and also on 
the axis near the screen. Using this approach, Braunbek recently 
solved the problem of scalar wave diffraction by an aperture in a 
concially shaped screen [31]. 


In the papers of Frahn [32, 33], this method was used for the 
diffraction of electromagnetic waves. Diffraction of a plane wave 
incident normal to an ideally conducting screen with a circular hole 
was investigated. The field was calculated in the hole and on the 
axis, and also the field in the far zone and the transmission coeffi- 
cient (the ratio of the energy passing through the hole to the energy 
falling on it) were calculated. 


In these works of Braunbek and Frahn, secondary diffraction was 
not considered. The expressions obtained by them for the fringing 
field intensity in the far zone agree with similar expressions 
following from our equations (§ 9). 


Karp and Russek [51] studied diffraction by a slit in the case 
when the incident wave's electric vector is parallel to the slit edge. 
They investigated each semi-infinite pert of the screen as a half- 
plane excited by the incident wave field and a "virtual" source 
localized on the edge of the opposite half-plane. The moments of 
these sources were determined from a system of two algebraic equations 
which were obtained by using the asymptotic expressions resulting 
from the rigorous solution for the half-plane. Secondary diffraction 
was considered, and partially the general interaction. Special 
attention was allotted to calculating the transmission coefficient, 
but equations for the scattering characteristics which would be 
suitable with all directions of incident wave propagation were absent. 
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Clemmow [46] and Millar (47 = 49] in their works calculated the 
transmission coefficients with normal irradiation of a slit and a 
hole, and also the field in the hole. The solution was sought by 
means of curvilinear integrals of the fictitious linear currents on 
the aperture edges. The interaction of the edges was considered. 

The case of inclined irradiation was not investigated, since it turned 
out to be too complicated for investigation by this method. 


The "geometric theory of diffraction" of Keller [42 - 44] which 
deals with diffraction rays is of special interest. The phase and 
amplitutude corresponding to each diffraction ray are determined at 
each ray point on the basis of geometric considerations and the law 
of the conservation of energy. The initial diffraction ray amplitude 
is assumed to be proportional to the incident ray amplitude at the 
point of its diffraction. The unknown proportionality constant be- 
tween the amplitudes and the initial phase difference -s determined 
from a comparison with the results of well-known solutions of diffrac- 
tion problems. In this way, the fields scattered with the normal 
incidence of a plane wave onaslit and hole in a flat screen are 
found. These fields are obtained with consideration of multiple dif- 
fractions, but they are not precise wave equation solutions, since 
their calculation was started from approximation relations. Moreover, 
geometric diffraction theory 1s not applicable near caustics, and also 
in the vicinity of the scattering diagram principal maximum. 


In a recently published paper of Buchal and Keller [52], a new 
method for the solution of diffraction problems for holes in a flat 
screen was proposed. The caustics and shadow boundaries here are in- 
vestigated as thin boundary layers, inside of which a rapid field 
change takes place. This method supplements geometric diffraction 
theory, and in particular enables one to find the field at caustics 
and on the shadow boundary. 


Recently, the method of integral equations has been applied to 


the solution of diffraction problems of holes in a flat secrecn. In 
particular, Greenberg [53, 54] reduced the solution of this problem 
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to an integral equation for a "shadow" current which 1s, in our termi- 
nology, half the nonuniform part of the current. The resulting 
integral equations may be solved (with any ratio between the dimensions 
of the hole and the wavelsngth) by the method of successive approxi- 
mations. Moreover, they allow one to obtain asymptotic expressions 
which are suitable for short waves. In Reference [55] Greenberg found 
an asymptotic expression for the current on a strip with ka >> 1 (2a 
is the strip's width). Greenberg and Pimenov [56] obtained a similar 
solution in the case of normal incidence of a plane wave on a circular 
hole. Using the same method, an asymptotic expression was found for 
the current on a flat ring [57], the width and inner diameter of which 
are a great deal larger than the wavelength. 


The above listed works [53 - 57] already relate to the mathemati- 
cal theory of diffraction: in them the first terms of the asymptotic 
expansions for the current were obtained with the desire evidently to 
also be able to calculate the following terms. Unfortunately, the 
asymptotic expressions which have been found up to now refer only to 
currents, and one is obliged to calculate the scattering characteris- 
tics by means of numerical quadratures [56]. As a consequence of the 
rapid oscillation of the integrands, such a method leads to rather 
unwieldy calculations and does not enable one to formulate a clear 
representation of the fringing field formation, and also does not 
allow one to study this field properly. 


Millar [58] investigated the problems of electromagnetic wave 
diffraction by slits in a flat screen. The system of integral equa- 
tions obtained by him for the current is solved by the method of 
successive approximations. The field in the hole is calculated from 
the currents which are found, and then on the basis of the field in 
the hole the field in the far zone and the transmission coefficient 
are calculated. All the indicated quantities are represented in the 
form of an asymptotic expansion in reciprocal powers of the parameter 
vka. A solution also is obtained in the case of glancing incidence 
of a plane wave. 
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Let us note that the asymptotic expressions obtained by the 
method of integral equations are distinguished by their considerable 
complexity, and frequently require tabuiation of the new special 
functions appearing in the expressions. 


In the recently issued volume of Handbuch der Physik [50], which 
is devoted to diffraction theory, the complex characteristic of plane 
wave scattering by a strip was studied directly, omitting the caleu- 
lation of the currents. For this characteristic, a singular integral 
equation was formulated, the solution of which was sought in the form 
of an asymptotic series in reciprocal powers of vka. The first term 
of the series corresponds to Equations (6.14) and (6.16). The follow- 
ing term takes into account the interaction of the edges, and becomes 
infinite with the glancing incidence of a plane wave and also for 
observation points lying in the strip's plane. Therefore, the simple 
expressions obtained in [50] do not allow one to construct the com- 
plete scattering characteristic. In [50] diffraction by a disk, a 
sphere, and an infinite circular cylinder was investigated, and also 
a review of the general methods of diffraction theory and a biblio- 
graphy encompassing a large number of works (mainly German and 
American) were given. 


The book of King and yu[59} presented (as a rule without deriva- 
tion) a series of asymptotic expressions relating to a slit and a 
circular hole and also to other diffraction objects. Here, however, 
equations from which one would be able to construct the scattering 
characteristics of a strip and a disk with any incidence of a plane 
wave also are missing. 


Works on diffraction by three-dimensional bodies having edres 
are comparatively scarce. In the paper of Siegel et al. [41], the 
effective scattering surface for a finite cone with the incidence of 
@ plane wave on it along the symmetry axis is calculated from clemen- 
tary arguments. The expressions obtained here do not fully character- 
‘fze the fringing field, and are suitable only for sharp cones to which 
we already referred in 5 17. In the papers of Keller [44], the 


\ 
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diffraction ray concept is used for calculating the scattering of 
scalar and electromagnetic plane waves by a finite circular cone with 
a flat base and also by a cone having a spherical rounding off instead 
of a flat base. The resulting expressions .are not applicable in the 
vicinity of certain irradiation and observation directions. In § 17 
we showed that the field scattered by a cone and by certain bodies of 
rotation is not expressed only in terms of the functions f and g, 
which refer to diffraction rays diverging from a wedge edge. This 
result evidently attests to the impossibility. of complete calculation 
of the scattering characteristic with the diffraction ray concept. 


Diffraction problems arising in antenna theory are usually dis- 
tinguished by their great complexity, since the corresponding metal 
bodies (mirror, horn, etc.) have a comspiicated shape. Since the dimen- 
sions of these bodies and the dimensions of the radiating apertures 
are considerably larger than the wavelength, the application of 
physical diffraction theory to antenna problems is very promising. 
only the first steps have been taken in this direction. Thus, Kinber 
[60, 61] performed a calculation of the decoupling and lateral radia- 
tion of mirror antennas. The feature specific to mirror antennas is 
that diffraction rays arising at the mirror's edge undergo multiple 
reflection on its concave surface. This multiple reflection was 
studied by Kinber in more detail as applied to the concave surface of 
a cylinder and sphere [62, 63]. 


Diffraction problems relating to an antenna dipole — a thin 
cylindrical conductor — are investigated in Chapter VII, and 
references to the literature are also given there. 


In conclusion, let us say a few words about diffraction of short 
waves by smooth bodies. The basic principles relating to such pro- 
blems were set forth in the fundamental works of Fok and Leontovich. 
These principles were established by the following methods of 
mathematical diffraction theory: 
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FOOTNOTES 


Footnote (1) on Page 138. These calculations were performed 


under the guidance of P. gs. Mikazan. 


Footnote (2) on page 150. See the footnote on page 86, 
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CHAPTER VI 


CERTAIN PHENOMENA CONNECTED WITH THE NONUNIFORM 
PART OF THE SURFACE CURRENT 


In the previous chapters, a theoretical investigation was 
conducted of the field radiated by the nonuniform part of the current. 
In this chapter we will discuss a method for measuring this field 
(§ 26) and we will investigate the phenomenon of the reflected 
signal's depolarization (§ 27), 


An experimental method for measuring the field from the nonuni- 
form part of the current was first proposed for bodies of rotation in 
the paper of Ye. N. Mayzel'ts and the author [12]. Later it was shown 
that this method has a universal character, and is suitable for 
measuring the field from the nonuniform part of the current excited by 
a plane wave on any metal body [13]. 


§ 26. Measurement of the Field Radiated by the 


ee 


Nonuniform Part of the Current 


Let an ideally conducting body of arbitrary shape be found in 


free space. A surface element of this body is shown in FBiraure 69, 


The coordinate system was selected in such a way that its origin would 
lie near the body, and the source Q would be located in the plane 

x #0. If the distance between the body and the source is a great 
deal larger than the body's dimensions, then the incident wave in the 
vicinity of the body may be investigated as a plane wav~. Let us 
reorssent it in the form 


Eg = E,,e8 Pen res coon E, 0. (26.01) 


Here y is the angle between the normal N to the wave front and the 
Z axis. 


Now let us place in front of the source, parallel to the radiated 
wave front, a polarizer P which transformed linear polarized radiation 
into a circularly polarized wave. Let the wave passing through the 
polarizer with an electric vector En lag in phase by 90° behind the 
wave with an electric vector E, (Figure 70). In this case, the polar- 
izer achieves a clockwise rotation » As a result, the incident 
wave field at the coordinate origin will equal 


i= 
‘ 
E.= 5 Eun H,= TE, . (26.07) 


The field scattered by the body may be represented in the wave 
zone in the following way: 


i-_ 

. a 1a 

B= — Hye VO ge 
pa A 








(26.03) 


where a is a certain length characterizing the body's size and (4,9) 
anc Z(A,9) are unknown angular functions. In the general case, the 


1 
es 
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Figure 69. 


The problem of elec- 


tromagentic wave diffraction by 


an arbitrary metal body. 


aS - is a surface element of the 


body. 


N 
wave front, 
Q - the source, 
P 
tion to a wave with 
circular polarization. 


ee 
45° . 
ep 
Figure 70. 


- the polarizer converting 
linearly polarized radia- 


- the normal to the incident 


field (26.03) is an elliptically 
polarized wave. In the direction 
this wave passes through the 
polarizer and creates behind it 
the field 


E.=—H ox “Fee res 
? Eus tae 
E,=H,= = Vr: 
(26.04) 
where 
v at /V. : 
=“. 2? (x=) (26.05) 


If the source radiates a 
wave of another polarization 


ey 


toward the source (O=*—1, 9=— > 


(H) | yoz), then the wave reflected 


by the body and passing through 
the polarizer is described at the 
point Q by similar relationships 





He=8, __ ite ith ae | 
ial, ." on { 
yam e | 


(26.06) 


Now let us investigate in the physical optics approach the 


diffraction of a plane linearly polarized wave by the same body. 


According to definition (3.01), 


ted on the body's surface by a plane wave with 


incident wave (Ey | yoz) equals 
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the 


165 


uniform part of the current excl- 


E-volarization of the 


f= — so Esa (n, sin y +2, cos pe’, | 


f=5 E,s-ng-sin ye, 


ome € . ad 
f= ge Eu n,cosye", | (26.07) 
and with H-polarization (Hy 1 yoz) 
ff =9, 
f= se Mee ree, 
= — Hf, n,-e 
hae Nene. | (26.08) 


Here Es and Hox are the electric and magnetic field amplitudes of 

the incident wave with E-polarization and H-polarization, respectively; 
y= kin'siny—-2'co7 1s the incident wave phase at the point (x", y’, 2’) 
on the body's surface; ny» Ny» n,» are the components of the normal 

to the surface at the same point. 


Furthermore, calculating the vector potential in the far zone 
on the basis of this current and substituting its values into the 


equations : 


E,=—H,= ik, ! 


E, =H, = ikA,, (26.09) 
we find the fringing field. With E-polarization, it equals 
A aR 
Ey Hy st Bye SE \tnssin 7 cos e+ 

Lin,sing + rn. cos Yisin 9] e*dS, (26.19) 

; eitk ; 7 . . 
Ey =H, = HE Exe Sg (Iota (sin cus 9 sing — cos sin ¥i — 
z A (26.19) 


— (a, sity.‘ Hy, cs p) eos Cos A] o'FAS, 


and with H-rolarization 


ee NC ee FIAT 7 


\~ 
aN 


( 
E=-H=> Hys 60s 9 ef n,e'*dS, 
By =H = Hye S[(ny sin ® + si : 
= HL= Hoe (nz sin ® + a, sing cos 8) eas. (26.11) 


Here R, 8 , @ are the spherical coordinates of the observation point, 
@=4y—krcosD , and integration is carried out over the 1lluminated 
elements of the body's surface. In the case of radar when the obser- 
vation and irradiation directions coincide (9==—1, e=—-y), 
Equations (26.10) and (26.11) yield 


i ak : 
£,=~-H,=—+ Eu (a, sin 7 + 
+1, 00s y) e* dS, (26.12) 
£,=H,=0 ; 


and 


ek 
Ey=H=H He Sin, sin + 


+n, cos y)e'* dS, 
E.=H,=0. : (26.13) 


Furthermore, assuming the incident wave amplitudes ar. specified by 
Equation (26.02), let us write Expressions (26.12) and (26.13) in the 
following way: 








am 
4 aR 
anaes, | 
En =H. a itEne ott gan ~ 
a yr Rk | (26.14) 
where 
Ba — B= pf imysiny + necos yeas. (26.15) 


Now let us represent the angular funetions of fringing field 
(26.03) in the form 
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a 


= YES, 


Sawer, (26.16) 


<0 


where the functions 7. r ri 


and rh X refer to the field radiated by 
the uniform and nonuniform part of the current, respectively. Substi- 
tuting these expressions into Equations (26.04) and (26.06) and taking 
into account relationship (26.15), let us find the fringing field 
passing through the polarizer P toward the source Q. In the case of 
E-polarization, 1t equals 


E,=—Ay= Mh (yr 4 8) ce 
é, Hs (alee (2s° 4 ~ SY en (26.17) 


and in the case of H-polarization 


_ ith IS 
Hy = Ey = ‘tes p ayy One? 
iah,, 


5 iaR " 
Hy=— E,= GUS EE BYE. | (26.18) 


The physical meaning of the result obtained is as follows. The 
field scattered by the body at the point Q is the sum of two waves 
polarized in mutually perpendicular directions. The reflected wave 
which is polarized the same as the primary radiation of the source 1s 
determined by the function t= (S24) » and is created only by the 
nonuniform part of the current. The reflected wave with the perpendi- 
cular polarization is described by the function %.==2%?-|-E'—%') and 
is the field radiated by both parts of the current. Let us note that 
in the general case the functions zi and zt do not coincide, and 
therefore they are not balanced out in the expressions for £_. In 
other words, the field radiated by the uniform part of the current in 
this case may not be separated from the fringing field. 


Thus, the investigated method allows one to separate from the 


total field scattered by any metal body of finite dimensions ‘hat part 
of the field which is caused by a distortion of the surface (the 


FTD-HC~2 3-259-71 18 


curvature, a sharp bend, a point, a bulge, a hole, etc.). One chould 
note that, in the case of electromagnetic wave scattering by a system 
of separate bodies, the separable part of the fleld is due not only 
to the surface's distortion, but also to the diffraction interaction 
of the bodies. 


It is necessary, however, to keep in mind that it is possible to 
realize the indicated fringing field distribution not in an arbitrary 
observation direction, but only in a direction for which the condition 
£9 = aro is fulfilled — for example, in the direction towards the 
source. 


Consideration of the nonuniform part of the current also enables 
one to explain the reflected wave depolarization which we will inves- 
tigate in the following section. 

(2) 


Figure 71 presents the results of measurements and calculations 


of the effective scattering surface 


at = natlS, [P= nat Et Bt, (26.19) 


which is dependent upon the nonuniform part of the current excited by 

a plane electromagnetic wave on a disk. The disk's diameter equals 

2a = =A (A is the wavelength). The calculations were performed with 
consideration of the secondary diffraction on the basis of the approxi- 
mation equations for the functions £ and © which were derived in § 24. 
Since it is difficult to prepare a thin disk with a sufficiently flat 
surface, the measurements were performed with an obtuse cone close to 
the shape of a disk and having a height approximately equal to one 


tenth of the dianeter. 


As is seen from Figure 71, the theoretical and experimental curves 
are fairly close together. A certain divergence between then, 
especially in the region of y values close to 90°, may evidently be 
explained both by the model's conical shape and also by the 


Dennen ee ee eee EEE IEEE tan ten ean 


(2) pootnote appears on page 174. 
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Figure 71. Diagram of the radiation from 


the nonuniform part of the current flowing 
on a disk. 


approximation character of the computational equations. The value 
y = 90° corresponds to the direction along the disk's surface, and 
the value y = 0° — to the direction normal to the disk. 


§ 27. Reflected Wave Depolarization 


Let us again return to the problem of scatterine of an electro=- 
magnetic wave by an arbitrary metal body. The relative position of 
the source Q, of a surface element of the irradiated body, and of the 
coordinate system 1s shown in Figure 69. Let us recall that the 
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source Q is in the plane y0z, and 
radiates a linearly polarized wave. 
Furthermore, we shall assume that 

° the polarizer P which is shown in 


* 
Figure 69 is now absent. 
& Let us designate by a the 
angle between the plane y0Oz and 
Figure 72. the incident wave electric vector 


Ey (Figure 72). The field of 
this wave will be represented in 
the form 


E,= H,, = E,,e8 sin 742 ¢08 3) 
H,=— E,,= 4,.0"" a+zros 7) (27.01) 


where 


Eur = Eqsin, Hya= — E, cos, p= tea, (27.02) 
The field scattered by the body is determined in the wave zone 
by the equations 


Ey = — Hy =F [Esta (1 8 9)+ 
ett 


+ Ash (1.9. 9)) Q- 





Ey=H, =F [Eeskalt. % 9) + 
ae 
+ HE, (1, 8. ge - 
(27.03) 


Here a is a certain length characterizing the body's dimensions, R, 
® , ¢ are the spherical coordinates of the observation point, 
£,.(7.9 9) and £,,(y,9 9) are unknown angular functions. 


It 1s obvious that the fringing field polarization — that Is, 


the orientation of its electric vector in space — depends in a con- 
plex way on the observation and irradiation directions. In the 


FTD-HC-23-259-71 171 


direction toward the source, it may not coincide with the polarization 
of the wave radiated by the source. Such a phenomenon is called 
reflected wave depolarization. 


It is easy to establish the reason for depolarization, if one 
investigates tno fringing field as the sum of the fields radiated by 
the uniform and nonuniform parts of the current. According to § 26, 
the uniform part of the current radiates the following field in the 
direction towards the source (92-1, =—+) 


£,=—H, =e 5 -E, 


‘\ 


iaE,, ** = 


icH,, #4? | 
E,=H =e a (27.04) 
The functions 7° and 29 satisfy the condition 29 = f°, and are de- 


scribed by Equation (26.15). From Equation (27.04), let us immediately 
obtain the equality 


fe 
Fy = 18% (27.05) 


which means that in the physical optics approach the reflected wave 
does not experience depolarization. Consequently, the reflected wave 


depolarization is caused only by the nonuniform part of the current 
or, in other words, by the surface distortion. 


Let us derive an equation for the magnitude of angle 5. This 1s 
the angle by which the electric field vector of the reflected wave is 
turned in respect to the electric vector of the wave radiated by the 


source. For this purpose, let us represent the functions T1(2) and 
E402) in the form 
Fo > 
B= Ey + Sar | 
=) =Sa + Bia’ | (27.06) 


where the terms Ty (3) and eee correspond to the fleld radiated by 


the uniform part of the current, and the terms ee and er 
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correspond to the field raaiated by the nonuniform part of the current. 
Comparing Expressions (27.04) and (27.03), we find that 


@ 
tI 


v >, z) =0, ) 
Bao, Lad. f (27.07) 


Therefore, the field scattered in the direction towards the 
source (@=«—1. y=— +). will equal 


E&,=— H,= Fiz (E°-+-5)4H,.5) Lat 
Ey= Hy = EES, — Hes (BP SS. | (27.98) 


This field's electric vector forms an angle & with the yoz plane. 
The angle 8 is determined by the equation 


_ Es _ S48 —Hetge 
I= BR poneaigs ©* (27.09) 


As a result, the desired angle 6 which characterizes the depolarization 
Magnitude will equal 


t=a— fp. (27.10) 


Thus, the field from the nonuniform part of the current, separatle 
"in a pure form" by means of a polarizer (§ 26), leads to depolarizatio: 
of the scattered radiation. 


Specific results from the depolarization calculation of wavec 
reflected from certain bodies may be found, for example, in the works 
of Chytil [75 - 77] and Beckmann [78]. In particular, in Reference [77 
it was shown that the depolarization effect on the effective scatter- 
ing surface of convex bodies in practice may be neglected only with 


the condition ka 2 4, 
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Footnote (1) on page 164. 


Footnote (2) on page 169. 
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FOOTNOTES 


A system cf metal riates parallel 
to the e, vector may serve as the 


simplest example of such a 
polarizer. 


See the footnote on page 386. 
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CHAPTER VII 


DIFFRACTION BY A THIN CYLINDRICAL CONDUCTOR 


In almost all the works devoted to the diffraction of plane 
electromagnetic waves by a thin cylindrical conductor, the current in- 
duced in the conductor was studied, and then, by integrating this 
current, the fringing field in the far zone was calculated. However, 
in view of the complexity of this problem, they succeeded in obtain- 
ing relatively simple equations only in the particular case when the 
observation direction and the direction toward the source coiricided, 
and was perpendicular to the conductor axis. In the general case 
when these directions did not coincide and were arbitrary, the expres- 
sions for the fringing field became very complicated and unsuitable 
for making calculations. Since they were obtained by integrating 
approximation expressions for the current, it turns out that they have 
still one other shortcoming — they do not satisfy the principle of 
duality. 


In this chapter, explicit expressions are obtained for the 
fringing field which are suitable for making calculations with any 
direction of irradiation and observation. We shall consider both the 
primary edge waves excited by the incident plane wave and also the 
secondary, tertiary, etc. edge waves. The total fringing field is 
found by summing all the diffraction waves. 
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§ 28. Current Waves tn an Ideally Conducting Vibrator 


The electrodynamic problem of determining the current in thin 
cylindrical conductors (vibrator) usually is reduced to an interro- 
differential equation. The latter is derived by means of boundary 
conditions on the conductor surface, and is substantially simplified 
in the case of thin conductors when the inequalities 


+ <1 and ka <1, (28.01) 


are fulfilled, where a is the radius and L is the length of the 


conductor and kat et : 
¢ 


Its solution may be found, for example, by the method of succes- 
sive approximations [79, 80] or by the perturbation method [85]. 
Recently, Vaynshteyn [81, 82] proposed a new solution for this equation. 
Since we will subsequentl; base our work on the results of References 
(81, 82], let us discuss them in more detail. 


Let us assume that the vibrator's symmetry axis coincides with 
the z axis, and its ends have the coordinates z = zy and z = Zo 
(L = Zo - z,). In the case of excitation of the dipole by a concen- 
trated external field 


Et = 3(z) (28.02) 


the current J(z) in the conductor may obviously be written in the form 
of the sum of the waves travelling along the conductor with a velocity 
e from the excitation point z = 0 and the ends z = zy and < = Z5° In 

Reference [81] it was shown that the complex amplitudes of these waves 
are slowly varying functions of the z coordinate. These functions 

may be approximately expressed in terms of the function wp(z), so that 


we obtain the following cxpression for the current J(z): 


S(ayHA ly iz yer! Ag (2 —2,eh Or + 
+ Ady (2,2) eM}, (28.03) 
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Here the quantity 





A= <- » Y= 1,781... (28.04)! 
determines the initial value of the current wave propagated from the 
excitation point ‘1), The function w(z) 1s the solution of the inte- 
gral equation, and in addition to the variable z it also depends on 
the parameters k and a. We will not list here all the properties of 
the function p(z), but let us note only that it satisfies the 
conditions 


4(0)-=1, p(co)=9, (28.05) 


and its absolute value monotonically decreases with an increase of z. 
This decrease, which is rather slcw and does not have an exponential 
character, is due to radiation. 


The constants Ay and A, determine the initial values of the 
current waves originating at the points z = Z4 and z = Zo» respectively, 
and travelling in the direction towards the opposite end of the 
conductor. These constants wre found from the conditions at the 
conductor ends 


vizi=ndiey = 0 (28.06) 
and equal 
A=- -bia(- 2— blegg bree. 
Ay Fy Hed b( zn e(Lye 8%, Neat 
where 
D==1--g'(Lye. (28.08) 


LA 


)pootnote appears on page 216. 
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Considering that the quantity 1/D is equal to the infinite geo- 
metric progression, 


pale tye +... (28.09) 


Expression (28.03) may be written in the exparded form 


F(Z) = Ale (lz pe” tl _ia(—2zye™ [y(z—2,) eter 
— 9 (bye y(z,— zy el gt Lye™* y (2 — 
— 2," J — oz Je (9 (2, — 2) ter 
—4o (Lye y (2 — ze V4 
+98 (Leg (2, — ze. I}. (28.10) 


The physical meaning of Expression (28.03) is seen from this. The 
first term in Equation (28.10) is the primary current wave which coin- 
cides with the wave excited by a concentrated emf in an infinitely 
long conductor, The second term (in all brackets) corresponds to the 
current resulting from the reflection of the primary current wave 
from the conductor end z 2 Za» and as a result of subsequent reflec- 
tions from the conductor ends which arise from this wave 

—Agl— ze "y(e— 2," | The third term (in all brackets) corre- 
Sponds to the current resulting from the reflection of the primary 
wave from the end z = Zo and as a result of the subsequent reflections 
Eeen the conductor ends arising from this wave ~-Jep(7,)¢"""p(z, ~ = zyeltteen2), 


It also follows from Equation (28.03) that external field (28.02) 
excites in the semi-infinite conductor (2:52) the current 


J(2y= JA be q2pe!?! —4(—2,)e oz —2,yeh (26.11) 
and in the semi-infinite conductor (--%<2<2z2,) the current 
Jz) =A (y(lz er! —ytz,e* yg (z,— 2p ey, (28.12) 


Comparing these expressions with the proper terms in (28.19), we see 
that the reflection of all the current waves at the end of a finite 
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length vibrator occurs in the same way as at the end cf a semi-infinite 


conductor. 


in the case of a passive vibrator (z,<27.). excited by the plane 


wave 


EL = Ege, a — ke cos bh, (28.17) 


the current also is represented in the form of the sum of waves (see 


[82]) 


J(z) =S [e“ = $i(z —2,) eit tient 
—, (2, — zeit 1 A a2 — 2,)e™ au 


A, $2, — 2) hy, (28.14) 


where the first term corresponds to the current excited by a plane 
wave in an infinitely lone conductor. Its complex amplitude S equals 


igus iwE ys ; 
aD Lee re ae (265.15) 


The second and third terms are primary edse waves arising as a consé- 
quence of the cut-off of the current set¥= | They are expressed in 

terms of the functions y + (z) and wv - (2) which depend, in additicn 
to the variable z and the parameters k and a, on the angle 8 . These 


functions satisfy the relationships 
7, )=1, 4, (20) =0, 


¥, (2)\ poe =. ¢. (2)| g29= (2) 


ke 
aN 
~~ 


9, (2)! peo 9- (2)1g = b- ee 


The initial values of the primary edge waves are Sucn that vieir sum 
2wz 


with the wave Se gives a current equal to zerc at the conductor 


ends. 


The last two terms in Equation (28.14) correspond to seccndarv, 


a 


tertiary, etc., edge waves, and have the same ferm as they do for 3 


transmitting vitrator (compare Equation (26.03)]. The unknown coeffi- 
efents x and k, are found from Canditions (28.06) and equal 


Ae + el" ded (L) -- 


-- 3. (L) a(Lye 


a-opl, she, 
.f 


Je 
Ase yey. tL) 


—%, (L) ai Lyell? oy et, (28.17) 


Using equality (28.09), Expression (28.14) may be written in the more 
graphic form 


J(2)= S{e'** ~—p.(2 ~2,) ele ta toned 4 
fy. (Lye eh" ty (z, — 2) efter 
~~ p(Lye**y (2 — 2,) my 
eyr(Lyett 4 (2, —z2) eM 
— be (2, — 2) eet 
te (Ly eet ty (z —2,) ef. 
—y (Lye y(z,—2) lM 
Pat (Lye Mt gz — ze — (28.18) 


Here besides the wave ge iz 


and the primary edge waves, which we 
talked about in connection with Equation (28.14), the secondary, ter- 
tiary, etc. waves diverging from the ends 2 = 24 and <= = Z5 are 
explicitly written out; they correspond to the first, second, etc. 


terms in the graphs. 


Passing to the limit in Equation (28,14) when Zo + @, we find 
the current in the semi-infinite conductor (215 «) 


Jz) 2: S-fetet —-y. (2 == zener (08.19) 


and, similarly, we find the current in the semi-infinite conductor 
(-2, Z5) 


\ J(2) = S-[eiv? — 4, (2p ale trey, (29.20) 
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It is not difficult to see that In the cace of @ passive vibrator 
the reflection of current waves at its ends occurs in the sume way as 


at the end of a semt-infinite conductor. 


Thus, the complex amplitudes of current waves ina thin, finite 
length conductor are proportional to the function: (2) and yw, (2) 
which monotonically decrzase wlth an increase of 2 as a consequence 
of radiation. Let us note several properties of current waves in a 
vibrator. Each advancing wave in sum with the reflected viave excited 
by it gives a zero current at the conductor's end. In the case 
‘Leva,--2, 2a-}-(n=1,2.3...) and D « 0, a current resonance begins in 


the vibrator. 


The precision of Expressions (28.03) and (78.14) obtained by the 
method of slowly varying functions is different in various sections 
of the conductor. It is comparatively low near the conductor's ends 
(and in the vicinity of the point z = C of a transmitting vibrator) 
where the current waves arise, and where their complex amplitude varies 
rather rapidly. As the distance from these vibrator elements increases, 
the precision of these equations increases without limit. 


It should be stated that with a more rigorous approach [7°, 20] 
the amplitudes of all the reflected waves will be determined by dir- 
ferent functions; however, the difference between them rapidly decreases 
with an increase of.the reflection number. The functions yw(s) anu 
v, (2) only approximately describe these current waves, but on the other 
hand they allow one to effectively sum them and to obtain closed 


equations. 


Using the variational method ror the functions p(s) and y,(2), 





we obtained the approximate, but on the other hand, simple cjyuntions 
(see [83)) 
=~! 
la 
v4 
= 
—_—— ~3:¢8 
la 1 E(2gxie-tie (08.01) 
(Equation continued on ‘ext vase) 
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{ ; 1, 

Y (z)=: Dies ae ee 

: ta E(ayx)e gs (28.21) 
where 

In(—N) = ix, Ini =ig. (28.22) 
° 

aatt, g=(tay', =F" (28.23) 
and 


E(y)=Ciytisiy=— (Sa. (28.24) 
, 


The integral cosine Ci y and the integral sine si y are determined by 
the relationships 





Cys— [Sih at, siya— (ear (28.25) 
rf eo 


and are thoroughly tabulated functions. 


The equations written above for the current in a finite conducter 
are distinguishable by their visualtsability, and they enable one 
to liken the conductor to a section of a transmission line in which, 
however, the attenuation of the current waves takes place, not accord- 
ing to an exponential law, but according to a more complicated law 
which is determined by Equations (28.21). In addition, the diffraction 
character of the problem is reriected in the equations. The conductor's 
specific features as a diffraction object are included in the very 
slow attenuation of the current waves. As a consequence of this, it 
is impossible to limit oneself to considering only secondary and 
tervlary diffractions, and it is necessary to sum all the reflected 
waves. As a result of such a summation, a “resonance denominator" 
D appears which takes into account the resonance properties of a 
thin, finite length conductor. 
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§ 29. Radiation of a Transmittin Vibrator 


The radiation characteristic of a transmitting vibrator way be 
calculated from a known equation by integrating the currents in it. 
However, such an approach is not advisable because, as was indicated 
above, the precision of Equations (28.03) is different in different 
parts of the conductor, and is low near its ends (z = Zy and z = Zo) 
and the point z = 0. The prineiple of duality gives more precise 
results. This principle leads to the following expression for the 
radiation field in the far zone [82]: 


E H ith (0) 
ia 2 sin 8 In - a é ; 


E,=H= i (29.01) 
The function 


1 (9) =!}— Ys (2,) elite! secon Oy pi(— 2,)e7teait feos Wo 
+ Byy, (Ly ett! “89 9) 2 Bay (Ly @ ite lt te08 (29.02) (29.02) 


is connected with the current (28.14) excited in a vibrator by plane 
wave (28.13) by the relationship 


J (0) == Sf (0), (29.03) 
The coefficients BL and By do not depend on the angle ®. 
Expression (29.02) enables one to trace the formation of the 


radiation. The first term (one) is the radiation field of an infinitely 
long conductor excited by a concentrated emf. Propagating in the 


direction 9 = 0, this fleld reaches the conductor's end z = =, and 
— being diffracted by it — generates a primary edge wave (the second 


term). Ina similar way, the primary edge wave diverging from the 
1 is excited (the third term). The last two terms 
in Equation (29.02) determine the waves arising as a result of subse- 


conductor's end z = z 


quent diffraction (secondary, tertiary, etc.). The amplitude: a and 


B, of these waves may be found from the conditions 


2 
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{(9)-= F(2)-29, (29.04) 


which means that the radiation of a finite conductor in the direction 
of its geometric extension must be equal to zero. These conditions, 
together with a consijeration of the relationships (28.05) and (28.16), 
lead to the system of equations. 


B, + B,y (Lye 7 2 4-20, 


By (Lye + B, = 4(2,) 0", (29.05) 
from which we find without difficulty 
! ¥ ‘ * hae _ 
B= tela) (hr ataye pens, 
By= alee 2th) (ee em, Cine 


Keeping in mind (28.09), let us represent the functions f(¥%) ina 
more graphic form 


f (8) == 1— yy (2, oP OO 4 (2) Cider ¢ 
* (). (Lye7tner | 4 (L) eft! a, (L) pitzeces t) as 
+> yt (L) ety (L) ew it cose é. J] ee 
mar 4. (—2,) @7 eat tens ¢) de " (= 2,) eittt-a x 
XE. (Lye Oo — (Lye og (Lye eet OL 
+ pt(Lyetg, (Ly ener, (29.07) 


where the secondary, tertiary, etc. waves corresponding to the first, 


second, and following terms in the brackets are explicitly written 
out. 


Thus, the field radiated by a transmitting vibrator arises as a 
result of multiple diffraction of edge waves at the vibrator's ends. 
Let us note in connection with this that the edge wave is diffracted 
by the opposite end of the vibrator in the same way as at the end of 
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a corresponding seml-infinite conductor. It fc not difficult to 
establish this by investijyating the radiation of a semi-infinite 


conductor excited by a concentrated emf. 


§ 30. Primary and Secondary Diffraction Ly a 


Passive Vibrator 


Let a plane electromagnetic wave fall at an angle % on a thin 
cylindrical conductor of length L = z, - Z4 and radius a (Figure 73). 
For purposes of generality, we will consider that the incident wave's 
electric field E, forms an angle a with the plane of the fifure. 
Then, its tanzential component on the conductor surface will.equal 


Et = E,,-e7* , (30.01) 
where 


Eu==Esin%, El=E,cosa, ty ="— k’cosd,. (30.02) 


The current induced in the vibrator by this field was investirated 
by us in § 28. As was already indicated above, Expression (28.14) 
which was obtained for it has a relatively low precision near tne con- 
ductor ends. Therefore, it is inadvisable to seek the frinceing field 
by integration of the current. Let us also note that the fringinre 
field found by such a method does not satisfy the principle of duality. 


We shall seek the scattering characteristic of a passive vibrator 
by starting from the follow. ng scattering picture which naturally 
follows from the previous results. An ineldent plane wave, beins dif- 
fracted at the conductor ends, excites primary edge waves which are 
radiated into the surrounding space. Being propagated alone the con- 
ductor, each of these waves experiences diffraction at the opposite 
end of the conductor and excites secondary edge waves. The latter, 
in turn, generate tertiary edge waves, etc. The total fringins field 
is comprised of the sum of all the edge waves being formed during 
sequential (multiple) diffraction at the conductor's ends. 
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In § 28 amt 29, we notea that current waves are reflected frori 
tie ends of a finite leneth eonduct r the same as from the end cf a 
semieinfinite conducter, and that the diffraction of these wavec at 
each end takes place in the came way ac at tre end of a semi-infinite 
ecniluctor. Vherefore, tne prinary edee waves may be found from the 
croblem of: scatterine of a plane wave by the semi-infinite conductor 
(a5 #) and the conductcr (-", 2,). Tne sum of such waves gives tne 
primary diffraction field - 





eae 
EN) a: H'" = «Ege FD), (30.03) 
where 
laut a ee 
H) a, 
; wll og ele sy <apaigsduce hye 
ON Sos Peasy Feost, - Eeos ity) 
ne 
i a7 89 _ gg itey (eur @ > cos dy 
—~"F (cost + cvs U,) O(a cos, & cos B,) : (30.04) ° 
<a The function (D(w, iw) may be 
ae calculated by means of the rigor- 
r % 2 7 ous solution to the problem of a 
i cicercan a semi-infinite vibrator (see [82] 


§ 3 and [83] § 4), and in this 
case, it satisfies the relation- 


Figure 73. The incidence of 2 ship 
plane wave on a thin cylindrical 
conductor; % is the incident 
angle. 


Dw, yyO(- egy tn Stn =° 
et ee oe (30.05) 


However, hencefcrth it will not be necessary to have the rigorous 
expression for the functicn 9. Let us note that Equations (29.93) 
and (30.04) are similar to Expressions (€.13) for a strip. These 
latter expressions do not take into account secondary diffraction. 


] 
cn 
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The secondary ede wave propasated from the end 2 = 4, Lo exe 
“ 


elted during diftractton at this ond of the primary current wave 
= Sy (2 —z,) eee ena, (30.06) 


where by y(z) we mean the functions obtained from the functions w, (z) 
by replacing ® by &. For the purpose of calculating the desired 
secondary waves, it is necessary for us, first of all, to find that 
external field which, when applied to an infinite conductor 

(--% 2%), Would excite the current (30.06) on its section (2,2: 9). 


For this purpose, let us study the current induced in an infinite 
conductor by the external field 


lwitn 2<2, 


Co c wer, >... _ = . . 
Eps Euelate cal malady ts (30.07) 


3 
Let us assume that Wo has a small negative imayinary part (Im We < 0). 
We may rerard the quantity Eye "ut as a concentrated emf which, in 
accordance with equation (28.03), creates in an infinite conductor 


CE ae Lyin lage, 
co . (30.08) 


Therefore, In accordance with the principle of superposttion, the 
total current created In the region ¢*,"z<e by tine external field 
(30.07) will equal 


, 


1 Ree eibiet 
Maye "(ape Gy elteree Mae, s 
‘Wita 
a “a 
iJ ~- Lise. git i + (8) ett ak | 
aka me (30.09) 


The resultiny intemrral may be expressed in terms of the funetdons 
w(z = 21) and vi (z - 2), and the corresponding relationships derived 
in § 2 of [82]. As a result, we find 
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: Feit: é _ 
J(2z)= a o(z —2,)e°O-” — 
Sh sta’ Zl 
— Set of (a2 et. 
2h sie? 0, la ————5— 
q4a cos > (30.10) 


Thus, it turns out that external field (30.07) excites, in 
addition to the wave y’, also the wave yw. In order to excite a "pure" 


y? wave, it is necessary obviously to apply an additional external 
field 


Ey = @,8(z —2z,), 





(30.11) 
such that 
te, . 
Sete leapt 
Sik sint “gta yko 
+ Ey y (2— 2,) e!* (se) =0. 
nea (30.12) 
Hence, 
= E,se!"e 
alc ean Tes (30.13) 
2ik sin os 


In order that the sum of external fields (30.07) and (30.11) would 


create the current (30.06), it is still necessary to fulfill the 
equality 


Ey, ioe, 
ea. oe ae ae 


rsa slat, 


aeons, eee aay 





27% sin? ), lo cLarianior whe 2k? sin? ), In 
qka cos ek 


(30.14) 


which determines the quantity 
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U 
ts -————__-- 


ry 
A yha cos 
Ey =e Pie . 


Consequently, for the excitation in an infinite conductor (with 


zZ> Zz) of current waves (30.06), it is necessary to apply the external 
field 





ry 

ha cos + 

EL = Eu 7] on 8(z—z,)— e'"a(2 2] ‘ 
In fa sind, 2k sin* > 

Iwith 2<z, 

Owith z>2z,. (30.16) 


2) =f 


In a completely similar way, one may show that the external 
field 








yka sin she i isses 
By Ew — 5? | tee) — ee oe)| 
: esti canes: nh ne 
In yea sin i, | Pit cost 5 
>e2 
rz,-2) 4] Pwith 2>2, 


() eel 2, 
eas (30.17) 


excites the following current in an infinite single-wire line (with 
z< Z5) 


= Soy" p> ~ Fe . (30.18) 
Now let us study the diffraction of current waves (30.06) by the 
semi-infinite conductor (-«, Zo). For this purpose, let us use the 
Lorentz lemma [4] 


PG BL-p ila =o. (30.19) 


Here j= —iwp,1k -R’) 18 the current of the auxiliary dipole with the 
moment Py which 1s located at point 1 with the coordinates (R, 8); Hy 
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is its field on the conductor surface, where the external currents 


35 are specified; E, is the field created by these currents at point 1 
(Figure 74). 


The external current 35 is 
determined by the well-known 
equation 





IP = — 4 In] (30.20) 
Pigure 74. 
in terms of the electric field E 
on the conductor's surface. In 
view of the boundary condition 


E,+ Ei =09 (30.21) 
we have 
e me 
iy — te Fa: (30.22) 


Furthermore, defining the dipole moment Py in terms of its field 
in free space (at the point x = y = z2 = 0) 
jar 
* me un fit © 
bes k*p, A sine (30.23) 
and changing from the magnetic intensity Ay to the total current 


J= OH, (30.24) 


induced by the dipole in the conductor, we obtain from the Lorentz 
lemma the following relationship: 


iar 


E, =H, = ay sin) £ \ Ef J(z) dz. 


“iw, eae (30.25) 
If the dipole p, is moved to a distance R >> z, ~ Zo» then the 
field radiated by it may be investigated on the section Za 7 2y of 
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the semi-infinite conductor (-2, Z5) as 2 plane wave. Then the 
current induced in this section of the conductor will be determined 


by the equation 


Jz) = {e'”* — ey, (2, res 2) e” (aa— ey ( 30 . 26) 
where 
Sx ’ w= -- k cos. (30927) 





24 
2k? sin® Alin jyha siad 


We will select the quantity Zo in such a way that, at a distance 

Z5 - Zp from the conductor end, the reflected current wave would be 
practically equal to zero (%.(z,—2)=90). Substituting the function 
(30.26) into the right-hand member of Equality (30.25) and taking for 
the quantity Ee the external field (30.16), we obtain 


boas 


2sinO lo at 


Ey =H. 
. yea sind 


” 


xs i Es fe™* — 9, (z,—z) Per dem 





prrerys aie r {6 sem — 9, (Lymer 


‘ mereoay * 
& 
nN f(r ore) % . * cae 
+E, sere — E, elt ters § eu ory ter-a)ee} (30.28) 


~—h 


An important feature of this relationship is that the interration 
is performed here not along the entire conductor (-«, Zo)» but only 
along part of it (-a, Z4)> where the function 4.(2,—2z) describes the 
current with good precision. The integrals here are calculated the 
same as in Equation (30.09). As a result, the field radiated by the 
semi-infinite conductor (-, Z5) will equal 


eftk 


R 








Ey = H,, = 


i 
210 Bia - yha sin 8 


(Equation continued on next. pare) 
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E,,e! (349o08, 
th (cos U-feos ty 1 


x 


dy le’ —¥ iL) efeetady tes 


sin? la ---— —-— 
E, eth eileen sees : 2 yha sin + ; 
ae CT ee ee a | (30.29) 


The terms in the braces having, the phase factor eit correspond to 
the desired secondary wave dlverging from the conductor's end z = Z5* 
Using Equations (30.13) and (30.15), this wave may be represented in 
the form 


(2, yo (ar 
EV (2) = HO (z,)= -— Lene = ae 


~ike, cos 
ri € 
2 sind to Tha slat 


(30.30) 


where 


egy > Ok 
Bg (z,) = —.- feos x 
ksin,(cosd + cos*)in 2 cog 


X | cos* e cos® + la —~, 9, (L)— 
° yea cos - , 


— sin! *e. sin® + la —L, #H) : (30.31) 


y-a sia 7 


In a similar way, let us find the secondary diffraction wave 
befine proparated from the end zg = 4° In order to do this, it is 
necessary to investirate the diffraction of primary wave (30.18) at. 
the end < = zy of the semi-infinite conductor (2,5 2<*e0). In this case, 
the principle of duality leads to the followin relationship 


Tein eth 
Big Hage OS [At Jtevae. 
2 : i@k ng 2 ( 30 - 32) 


which, after substituting the funetion (30.17) and the current 


iwE hs _ 
J(2= 4 Ni --—"—y ae {e“7—p_(z--2,) eit? bib tz “ 
2k? sin? 8 la 


, 


7ha sie : peyaas 


Pea Fee F | i 


in it gives us whe field radiated by tie cemi-infinite eonductor 
(2)> ), The wave radiated by the conductor's end is the desired 
secondary edge wave and may be represented in the form 








* 12) (2 ie, 
ED) (24) = HO (24) = — SE Se etenen®, (30.34) 
sta 0 tn Tha ov 
where 
fete Hal ’ 
GZ=a— pms. tal cere x 


h sin 0, (cos 8 +- cos 3,) lo saat, 


2 2% in? 8 eee fee 
x fe 3 sin* > In ron s y- (L)— 
» Be east”. tn ——!- 9? : 
aes is yea cos #14] (30.35) 


Otherwise, this expression may be written directly by replacing, in 
Equations (30.30) and (30.31), z, by 24> 8 by «—8@ and & by <«—4. 


§ 31. Multiple Diffraction of Edge Waves 


The secondary waves (30.30) and (30.34) which were found are the 
waves diverging from the ends of the semi-infinite conductors (-°, Z>) 
and (Z15 o), If one excites an infinite single-wire line bv the 
external field 


EM os: 5 (2 2,), (31.01) 
where 
. dine ONL 
B= GID (2 )] aif icc) eens eitl - ike, cos O 
2 ales k sin 8a —— oe : (31.02) 
fin 0 lf oir sin ds j1.0eE 


then a spherical wave arises which with @=g coincides with wave 
(30.30). With the excitation of an infinite line by the concentrated 


emf 
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E’= gP0(z —2,), (31.03) 


where 
2i tap $9, (L) 
by yea —lhycon%, 
SESE eae (31.04) 
: asia +10 ia stab, 
a wave arises which coincides with wave (30.34) when @=0. It is 


not difficult to see that these external fields actually excite in an 
infinite single-wire line current waves which are equivalent to the 
secondary current waves in a passive vibrator [that is, equivalent to 
those..waves which are expressed by the first terms in the brackets of 
Equation (28.18)]. Therefore, the tertiary waves may be investigated 
as edge waves radiated by the semi-infinite conductors (24> @) and 
(-2, 2) with their excitation by the external fields (31.01) and 
(31.03), respectively. From Equations (30.25) and (30.32), we find 
without difficulty the total field radiated with the indicated 
excitation by the conductor (245 o) 


aa 
E,= HP — pe gn] (icon) 
2sin 8 In aut . 
and the total field radiated by the conductor (-2, Zo) 
E.= A= of a 2k fe" 4, ( L) eftteblaly, 
"2slo 91a yea sin (31.06) 


As a result, we obtain for the tertiary waves diverging from the 
ends Zy and 25 che following expressions: 





. 80) (2) ott _ ins, cost 
E%(2,) =H (2) = ——— "a eee, 1.07) 
ou a asia Bin asia (3 
EV (2,) far 
E® (2,)=—= HW” 2, me ge = Sig witmecen® 
PO Nie oe » (31.08) 


where 
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Bey =— Pg ine, | 
BM (z,) = — 6-4, (Lye { (31.09) 


In the directions toward the cpposite end of the conductor, these 
waves are equivalent to the radiation of an infinite single-ewire line 
excited by the external fields 


Ev = GP Hz —2), B= G(2)] (31.10) 


E* = @).3(2, — 2), G; = G2), (31.11) 
Consequently, the quaternary waves again may be investigated as 

edee waves radiated by the semi-infinite conductors (-2, >) and 

(Zy5 ) with their excitation by the external fields (31.10) and 

(31.11). Using the reciprcecity principle, we easily cbtain 





OYE tea ‘ 
EW (2, = H®(2,)= ee) Speen 
2sin Btn Yasin 
aR 
EW? (z,) == HO (24) ——S ed _ EE etter, (32.12) 
2sin ON asin | 
where 
6 (z,) = — GP -9_(Lye, 
EM (2) = — 6-9, (De, (31.23) 
In a completely similar way, the nee order edge waves 
BO) (z= Ht (2,) = OE, _ SER g—tsscns 
° = i 4 R ° 
Ssia dla Yea sia 
Ain + ikR EA 
ae (z,) = Ht (z,) Sy —2ne ~ a ittaces a 
asin 8 lo pean o (31.14) 


t 
are found. Here } 
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6) (z,) = — BE-"-g (Lye, ) 


GI (2,) = — G9, (L)e iat | (31.15) 


and 
= '"(z,) is Bin = BE (2,) Ione” (31.16) 


Thus, the field arising with multiple diffractions (starting with 
the second) may be represented in the following form: 


SUES (ed + EY ea = 
aw? 
= ——_;- +3: B"'"2z, Jeter co Oy 
2sla 3 1o — cha a 
* a) ~— hs, cos 8 
| + Ee" as "|: (31.17) 
where 
‘ f 
¥ one) = 
and 
= 8 (2,)-+ 1-9 (Lyet— ay Eo et, (31-18) 
y EUV (2,) = 
a= - 
= GIN en-+ (Bg (Lye — BYE et, (31.19) 


and the functions @(z,;) and 6% are determined by Equations (30.31), 


(30.35), (31.02) and (31.04). We will not write out here the rather 
unwieldy final expression for this field, but we wil] proceed with 
a calculation of the total field scattered by a vibrator. 


§ 32. Total Fringing Fleld 


Before beginning the derivation of the expression for the 
scattering characteristic, let us make the following observation. The 
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functions 9 which enter into Equation (30.04) Satisfy relationships 
(30.05), and may be found by factoring. However, our investigation 

of the Successive waves arising with diffraction at the conductor's 
ends was approximate. Therefore, it makes no sense to use the precise 
Expression (30.04) for the primary field. We shall use the approxi- 
mation expressions for the function 96 


i 





®(—k cos ®, — kcos %)= In ——,-----, 
tha sin ~3- sin + 
O(kcos 9, kcosd,) =: In — yy 
yha cos ~ “9° cos r (32.01) 
' 


which were obtained by the variational method and have a precision 
which is sufficient for our purpose (see [83]). More precisely speak- 
ing, we will use approximation Equations (32.01) in conjunction with 
the rigorous Expression (30.05), and we will set 





lo 


®% 
\ , yka cos “x cos 3 
Oost, Skea) a 
1 8 Seasin® "ka sind, 
[nee sen 


] 
; qka sin | sia -> 


(a cos 9, Revs so 





=| 
| 








sata ™ searare (32.02) 
Then the primary field will equal 
a? iE . ‘git 
eal Teac Sai ni oa 
“(O89 cos tin a sin tg!" Sha sin ®, 
x ctg #e ty + In - weed —,'. aye oT (eo Bon %) 
yea cos “gr cus: r 
= tg mT ty _ In bes @ ikea (cos 8+4cos 0) P 
yea sin-> sta —y- (32.03) 


Now summing Expressions (31.17) and (32.03), we find the total field 
Scattered by a passive vibrator in the form 
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ibe 
E,=Hy=— ESF (0, 8). (32.04) 


where 





F(8, §,)= St —= x 


(cos ® + cos 9,) sin 9 sin % In - yin: 


ai 
yea sin ya ta ne, 


ens (cos 6+<os 8) __ 





; ® 
X Jos" —- cos? oe In 
yha cos “yr cos iy 


8 i ike, (cos @eos 8 
-~ sin? - sin? <3 a laisse = wore ika, (cos &+cos yee 


jha sin ya ry 


it 4 i 
‘Pe i ak -s sin? Be In uae 
qa sha 
® a, i — it a,con,+ 2, 
— cos?—- cos’ -}- Io ‘ i-% ew Nae oe! 
tka cos -3- 


+e" |sin’ = sin® oe tn d 


0 
a 
qhe sia 5 
8 8 _ 
: “asa e * cost In i , vile ik (2,003 Be, coe 8) 
eras 


cos @ ++ cos #,) la —— 
iach aii tie ad Oe 1 gith—itzy cos Oy 
Py.) 


me . en its.cos oy pens cos 





cos ® + cvs 8,) la —— 
icos BF cos Bo) In eg ent IkL—ihe, cos Oy _ 
20 ey ve 


Se i aa agement, 7 (32.05) 





in which all the functions v, and ve have the argument L. The result- 
ing expression, despite its complexity, has a clear physical meaning. 
Actually, the first term in the braces corresponds to the primary edre 
wave radiated by the conductor's end Zz = 243 the second term corre- 
sponds to the primary wave radiated by the conductor's end < = Zo- 

The terms included in the first pair of brackets refer to the secon- 
dary wave departing from the end 2 = Zi» and the terms in the second 
set of brackets refer to the secondary wave departing from the end 

Z = 25. The remaining terms describe the sum of all subsequent waves 
arising with multiple diffraction and have a resonance character. 

The resonance begins with L = Zo 7 2, % n-a/2(n = 1, 2, 3...) when N ~ 0. 
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Another important feature of the scattering characteristic is 
that it satisfie: the reciprocity principle — that is, it does not 
change its vaiue with a mutual interchange of % and . One may 
also show that the vibrator does not radiate in the directions along 
its axis, and it does not scatter electromagnetic waves with rlancing 
{rradiation, that is, 


F(0, 8) = F(x, 0) = F (8, 9): FC, 2) ==0. (32.06) 


Furthermore, using representation (28.25) for the functions yp 
and ¥,, we obtain the following expression for scattering character- 
istic (32.05) in the direction of the mirror-reflected ray (0=x—0) : 


k 
F(z—9), 8) = — —_*,—+ 
Ain yea slo a, 


% 
HVE (a sin? s) teu re (2kt< cos? +) 





eee =< 
é(1o yea rast) 


’ t ( 
+ a). eee a rT sin®, 2 rt 
(: m0, (n og so am) 
+ + y cost 42 —~ In ae gene mos ta 
tka sia 
+ + sin? “© — tn —+-,- go ne ftE Mtene Ob 
: pee Cos -B- "s 





i 
ta-— . 

TRE tht 40 , idk (1 —coe ® 07,0 - 
Dp & Ip,ve aie, YP 14 


"8 tea ust, ,0 jaZ (: cos 8 0 f 
-=>- [y" se (: cos ye . (32.07) 


With glancing irradiation of a vibrator, when 0,=—0 or %=+r, it 
follows from this that F(s, 0)=F(0, 2):=0 


Now assuming that = 5 in Equation (32.07), we obtain the 


relationship 








Equation continued on next pare 
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5 In ch a 
fethoe 
(32.08) 


which characterizes the reflected field magnitude with normal 
irradiation. 


Let us also write the expressions for the function F(®, 8) which 
corresponds to the radar case when the observation and irradiation 
directions coincide (8=%) 





F(8, 8)= 
= ————— cos? 5 i ee Tae cot 8 
(sina b= qha a r iss) tha cos? -3- 
—sin' > In fone atk = ew ike con oy 
yku sin? -,, 
2 ig 5 In ——-= = _(L)— 
. yea sin = 
— cost In --—L—— p, (Ly [eit treet on eat . 
yea cos + 
— 2280 In atl) flip tener rae 
. out Ins malt (Lye ™** cos 4. 


piwemnne (32.09) 


We may show that when Iaz Equation (32.99) leads to Expresston 
(32.03). 


The scattering characteristic (32.05) which was found above was 
obtained by summing all the waves formed with multiple diffraction. 
Such a method is very graphic, but somewhat lengthy. Cne may arrive 
at the same result more quickly if it 1s assumed vhat the edge wave 
diffraction process at the passive vibrator's ends takes place, start- 
ing with tertiary action, the same as at the end of a transmitting 
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vibavetioae, Mherefore, the positive vilbbatorts: ceantertae bsqeran ney 


be ddreet ly sought fir the Cons 


F (% 9.) — = a “4 f (9, 8,). 
S000 StH AU a cD V0 ea sind, (seed) 


where 


, 


119.) x 


co5 8} cond, 


a a “ov Oye 
‘w Jeus? py cust ey In t a ° eee Bf coe Ba 
ya con ‘y vue 2 : 
a eek re 
coe si ; sin’? in ‘ “ye shen Boe on 
(ta skn-., ste - 
t : | e 
-fe “| sin? y sito In : a Y_tL)-- 


° 
yea stn 2 


a a ; oo paw canal 
se cus? if cos?) In vs ’ a) Mentunre ta bie aa 
yhacos “y - 


pet faint pining any 
tha sla 2 
aes Secon te In ms v, ¥, (4) |e ea I 
yfa cos 2 


| Cy : (Le Ee | | Ch, (1) Pe . 


The bast we Germs tao Equation qae.d by are the sam of atl edb 
waves start ime wlth the tert lary waves whbete are preopemeitiod reer ttre 
eortduetorlss emda ua = a itt oF om yy Penpeet Pwedy . Me ormst antic wy 
bh Coy ctme de bernined: from dhe eomd bt fon 


AL Uat ee Fe 


1 in ery a) 
Whteb Tends te the oastem or ce quent bor 
oy who as { Psat cy ROY 
yp tdye™ (Cy Im rade : | 
oy at! t @ -4e Pa) ee { beers) 
C, Pe yeibbe y it ra (Lde ‘ 


Crom witteh 


Cy enna e -In- i [yi gern treoete 


“D yka 
ze oe eo ites }. 
t cee ‘ ° al~ihaeont, 
C,=— ae" ln re — 
— yi eTtreonte (32.14) 


§ 33. A Vibrator Which 1s Short in Comparfson With 
the Wavelength (a Passive Dipole) 


The theory of plane wave scattering by a thin cylindrical 
vibrator which is discussed In this chapter is based on a number of 
physical considerations. One good aspect of this theory is the fact 
that its precision inereases with the length of the vibrator, since 
the current waves whose diffraction we are investigating are expressed 
more clearly, the longer the vibrator. fowever, one may also show 
that for short vibrators, the length of which Is small in comparison 
with the wavelenrth, the equations derived by us have good precision. 


It is clear that a vibrator which is short in comparison with 
the wavelength acts as a dipole, creating a fringing field 


a 
E,=-Hy= — kD, ip sin 8, (33.01) 


where the dipole moment p.. may be calculated by solving the electro- 
static problem. This dipole Moment depends on the dimensions and 
shape of the vibrator. In accordance with [92], the dipole moment of 
a cylinder in a uniform electrostatic rmleld E., equals 


pe =Dit)( ip) Fe (34.02) 


where D(7) Js a dimenstonless funetton ~ = 1/32 whieh is shewn dn 
Fivure 75 by the cont tnueus curve, With 2 >> 1, one may ealeulate the 
function D(Z) by means of the asymptot ie expansion 


Dil) = t(a ty +. \ Q, =: 2( Indl - x): 


wy 
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If in this expansion one limits oneself to the first term, then 


D()\=——_—, . 
3 (tju—-5) (33.04) 


The results of numerical calculations based on Equations (33.04) are 
shown in Figure 75 by the dashed curve: we see that the latter 
equation gives a good precision already with 2 2 9. 
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Figure 75. Graph of the function D(Z)} which deter- 
mines the cylinder's dipole moment. 


Thus, the dipole moment of a vibrator which is short in conpari- 
son with the wavelength equals 


ees oe, (33.08) 


and its scattering characteristic must have the form 
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F(8, 0) = 40 a al 11-of( nz) “}. 


(33.06) 


In this section we find the first two terms of the expansion of 
the vibrator's scattering characteristic F in reciprocal powers of 
the large parameter L/a (with \ + ©), and we compare them with 
Expression (33.06). We shall limit ourselves to the case §=%= 5, 
when the function F is described by the simpler Equation (32.08). 


With small values of the argument z, the functions (z) and 
$(2)= 94), .2= los {see Equations (28.21)] may be represented in the 
form 


= (— #0) i 
t221— ae) 
b(2) = | ~€(2)— 20) on 
H2)= 1-H FO) +0(ar). | (33.07) 


The functions g and g included here are determined by the equations: 


£(2)— 2 (0)= 10 pe Ea, gO) = Ine (33.08) 


@ 

(2) 71) sa tn We is fe ‘e = yy 

8 (2)—g(O)-= in BF pe | Fda, B= In Se. (33.09) 
& 

Let us note that Expressions (33.07) completely agree with the 
corresponding terms of the aysmptotic expansion for the functions ¥ 
and ~, which may be obtained from the initial integral equations which 
determine these functions (see, for example, [81], § 4). 


Limiting ourselves in the expansion for the functions (z) and 
$(z) to terms of the order of (xz)3, we have 
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2yhe 
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[ine (in ane —1)+ 


pe (in WF) =i ee (in vez = *)}: 


In addition, terms of the order (in as) 
(33.10) and (33.11). 


Equation (32.08), then one should omit terms of the order (in aw) 3 


Now if we substitute these expressions into 


(33.10) 


(33.11) 


are omltted in Expressions 


it. Therefore, the function aes +) may be represented in the form 


F(y +)> wey y2 


2i t t 
edad 
tea 


O(a tn 


— E(kL) |+ 


~ 25 (tio 24 ot 


al 
tipper we}: 


(33.12) 


Furthermore, taking into account Equations (33.10) and (33.11), 


we find 
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(33.13) 


(33-14) 


(33.15) 


and finally 
i 
2ia-—— 
‘ qka i 
re eapete y (Lye = 
od (aie aed tad ie a 
— Fo TE EE tin) en BE AS 142) + 
- ae Le nope l4 , Sind 
+ig RPL In —— — hPL i? im Ps (33.1€ (33.16) 


Using these relationships, it is not difficult to show that 


“(> +) {c+ z (zy toll’ ty}. (33.17) 


The equation which has been found may be rewritten in the form 


5. 8) ol fin ZV. 
PCy 7) afm — >) pray ‘) h (33.18) 


It completely agrees with Equation (33.06) which follows from [92]. 
This result confirms the correctness of scattering characteristic 
(32.05) calculated by us, and shows that it 1s applicable for vibrators 
of any length. 
§ 34. The Results of Numerical Calculations 
The function F(0, 0) enables one to calculate the integral scatter- 


ing thickness S and the effective scattering area o of a passive 


vibrator. The integral scattering thickness is determined by the 
relationship 


° P 
3=F° (34.01) 
where 
pee 
ne 9 (34.02) 
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is the energy flux density in the incident wave averaged over an 
oscillation period, and 


a ‘ 
P= ERe fnieHejas— +. Esin 9, Re| J(ay elt ee gz (34.03) 
&% 
is the value of the energy scattered by the vibrator into the surround- 


ing space averaged over a period. Since one may represent the fring- 
ing field in the far zone in the direction ®=x—® by the equations 


‘ , % 
E,=H,=— * sina, J (2) elt# °° yp (34.04) 
. a 
and 
a < 
E,=H, =~ ESF (20, 8), (34.05) 


then, having determined from this the integral 


a ‘ 
Jaye rds = oe F(x—9,, §,) (34.06) 


BS) 


we obtain 
S= * cost alin F(e 9. 0,). (34.97) 


Calculations of the quantity S/Le (with a=0, %&=:-) ), performed by us 
for vibrators with the parameter tee taking the values xy = ~0.05 
and x = -0.1, are found to be in agreement with the results of 
Leontovich and Levin [85]. With x = -0.1, our curve (the dotted line 
in Figure 76) is only slightly displaced in the direction of longer 
wavelengths and gives slightly higher resonance peaks. 


The effective scattering area o according to the definition 
equals 


(8, 9) Seek (34.08) 


where p is the known quantity (34.02) and 





Figure 76. The integral scattering thick- 

ness of a vibrator as a functior. of its 

length (with normal incidence of a plane 
wave). Curves 1 were calculated by Leontovich 
and Levin [85]. Curves 2 were calculated 

on the basis of Equation (34.07). 


represents the average value of the energy flux density scattered by 


the vibrator in the direction fh . Consequently 
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Eytan te (34.09) 


+ 2(0, 0) = 7 cus? a-| F(, 8) I". (34.10) 
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If the receiving antenna operates with the same polarization as 
the transmitting antenna, then rhe corresponding value of the effec- 
tive area es equal 


4 . 


a 9,00, y= Acosta 10, OI" (34.12) 
pk Meee OBAG i 

In i ithe ease my, sedan when the transmitting and receiving antennas are 
combined arid’ the polarization is arbitrary, the vibrator's scattering 


properties ‘are characterized by the average vanue 


ce tof eo. 9,)da= P| F (9, Ot. (ana 
| - = ; 


In Figures 77 and 78, the results of calculations performed on 
the. basis.of Equations (34.12) and (32.08) for the case of normal 
irradiation (*=+) are represented by the ‘dotted lines. Figure 77 
illustrates the sepencence of the function 6 on the quantity kL with 
a given: value of a, =2107% = 15, — that is, when the ratio of the 
vibrator’ s length to its diameter equals L/2a = 452.. In Figure 78 
the graph of the function ¢ is constructed as a function of the fre- 
quency f = c/A°10” (in megahertz) for the prescribed parameters 
L = 5 cm and a, = 15. Here the curves plotted by Lindroth (79] are 
drawn with a continuous line, and the curve in Figure 77 calculated 
by Van Vleck et al. [86] is traced by the dash-dot-curve. 

-. 
The curves of Lindroth and Van Vleck were calculated by inte- 
grating the current which is found as a result of the approximate 
solution of the integral equation. However, this procedure was per- 
formed in [79] and [86] in a different way. Lindroth obtained an 
expression for the fringing field in the form of an expansion in 
reciprocal powers of the parameter 2_. The expression includes terms 
of the order of a In [86] a different kind of approximation was 
used which led, as can be seen from Figure 77, to rather rough result 
especially in the resonance region. Our curve (the dotted area) 
agrees almost everywhere within the limits of graphical precisicn wit 
the curve of Lindroth. A noticeable divergence is observed only in 


the magnitude of the first resonance peak. 
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Figure 78. The effective scatter- 
ing area of a vibrator as a 
function of the frequency 
f = c/A-1079 (in megahertz) with 
normal incidence of a plane wave. 
The designations are the same 
as those in Figure 77. 


INS SS 


N. 


4 In Figure 79 and 80 radar 
alagrams are constructed for vi- 
brators of a length L = 0.5. and 


Figure 77. The effective scatter~ Ll = 2\ with the specified value 


ing area of a vibrator as a L/a = 900. Curves 1 were calcu- 
function of its length with 


normal incidence of a plane wave, lated by Tal using the variational 


Curve 1 was calculated by method [87]. Curves 3 were ob- 
Lindroth [79]; curve 2 was 


calculated by Van Vleck [86] by tained by the method of induced 


means of the method of integral emf [86]; curves 4 were obtained 
equations; curve 3 was calcu- 
lated on the basis of Equation by the above-indicated method of 


(34.12), Van Vleck. The results of calcu~ 
lations based on our Equations 
(34.12), (32.08) and (32.09) -age 


ae shown by curves 2. 
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Figure 79. A comparison of the diagrams 
for the effective scattering area of a 
half-wave vibrator calculated by various 
methods. 


Curve 1 was calculated by Tai [87] by the 
variational method; 


Curve 2 was calculated on the basis of 
Equation (34.12); 


Curve 3 was calculated by the method of 
rhedie” emf (in the work of Van Vleck 


Curve 4 was calculated by Van Vleck [86] 
by the method of integral equations. 


In the cited references, the fringing field was calculated by 
the direct integration of the current. In order to determine the 
current, various approximation methods were used. In the variational 
method [87] a functional was constructed for this purpose which was 
stationary in respect to small current variations. Then the current 
was sought in the form of some function containing undetermined con- 
stants. These constants were found from the condition of the func- 
tional's stationarity. This method enables one to rather easily 
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Figure 80. Diagrams for the effective 
scattering area of a vibrator calculated 
by various methods. The designations 
are the same as those in Figure 79: 


L 2 2A 


obtain the first approximation; however, its results, especially for 
long conductors, may depend in a substantial way on the form of the 
trial function. In the induced emf method [86], the current is 

sought in the form of a combination of trigonometric functions with 
unknown coefficients. These coefficients are determined by using the 
law of conservation of energy. This is the simplest method, but it 

has a number of serious defects. Thus, as a consequence of incorrectly 
accounting for the current component having the incident field phase, 
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it leads to inaccurate results in the case of odd resonances (cspe- 
cially for long conductors), and it does not give the displacement of 
‘the resonance peaks from the values \ = 2L/n (n = 1, 3, 5....) in the 
direction of longer wavelengths. 


The results obtained by us are also approximate. However, our 
Equation (32.05) satisfies the reciprocity principle, and is applicable 
for any length vibrator. For very short vibrators L << A, it changes 
into the asymptotic expression for the scattering characteristic of 
a dipole (see § 33). For vibrators with a length of several wavelengths 
(L ~ndA,n #1, 2, 3, 4), Equation (32.05) gives satisfactory results. 
Calculations performed on its basis for radar reflection with normal 
irradiation agree with the results of Lindroth. Good agreement is 
also observed with the results of Leontovich and Levin for the inte- 
gral scattering characteristic. With an increase of the vibrator's 
length, the precision of this equation increases, and in this way it 
is favorably distinguished from the equations proposed for the 
scattering characteristics by other authors. 


Moreover, the divergence between the various approximation 
methods indicates the necessity of performing rather detailed calcu- 
lations based on precise methods, for the purpose of evaluating the 
actual error of the approximation methods. Such calculations may be 
performed, for example, by means of the method discussed in 
References [88, 89] or [91]. 
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FOOTNOTE 


Footnote ‘1) on page 177. Let us note that one may refine 
Equation (28.04) by multiplying its 
righthand member by the factor 5 


(usually or x 1) caleulated in 
Reference [84], 
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CONCTUS LON 


In this book, the gotutfon of aonumber of di ttract fon problems 
wis obtained based on the approximate: consftderat bon of the theta pera 
furbat ton fn the vietotty of ao eharp bend of the surtiaee or a aheurp 
ede. Equations were derived for the seatter ine eharactertout bea, ar 
In certaln eases (Chapter IV), for the radar pefleet bon thbeknesses 
with a speetfiled Leradtation direction. The espresstons whieh have 
been found have acelear physteal meaniog. They satdefy the ree tproe bby 


prineipte, and Chey are eonventent Cor mak dae ealbeubat tors. 


The results obtatned enable us to form a onere complete concept. 
or the applfeabliity Pimtta of the physteal optles approgeh. ot ts 
Usterl ly customary tooassume that this approaeh elves: rediabbe pesults 
eoly PP the bedyts diinens fons are Taree ti eomparhoen whth the waves 
Temeth.  Sueh oan ooptiden is based on the Collowlme: arecument. The 
Physteat opties approgeb eomsfderss onty the radiat fom trom the 
untform part of the current, and does not Tnelude in the ealeulat lone 
the nonunfform part of the current whieh ds eoneentrated tn the 
Vietnity of the bends and the sharp edges. ‘therefore, when the body fs: 
dimens fons are eons bderab Ly Daaeer) tdian fhe owavelometh, Che nenuab torn 
part oor the current: oceeuphes aovedattvedy saath part of the body fs: 
surface. Therefore, one would think that bts laf bacitee womkd be 


smell. 


Mut foo aetuadbiéy Pt turns out that the Pediatr t Pty of phys tect 
opt bes. results: depends substan! fatdy, met oondy oom tte body fs dimer. 
Sfons, but alse on the bodyts shape (nd fhe Prerad bert bert and: abeseryert pony 
dfreet bons. For example, with (te ycbanme fies fie bderee ef a weawe cor 
the Plat Paee of ao body, tlie edge sote ocenp bed Dy fhe qemu tb terre 
port oof the current is constderably broadened aad othe oetfeet oor this 
current: becomes substantial. Therefore, plysdeal opt bes retvers quadli~ 
tatively Encorreet results: for the @leld sent dered Py Plat ptertes 
with plane tis Irradiation Pndependent of the vatdio between their 


dimen fons and the waveleneth. The efteed ef the nomad form part oot 


Sd oe eS 


the current becomes noticeable also in those directions where, accord- 
ing to physical optics, the fringing field must be equal to zero or 
have a small value. 


The problem of diffraction of a plane wave with its incidence on 
a cone along its axis (§ 17) serves as a clear example of how impor- 
tant the above-indicated factors are. Although in this case the non- 
uniform part of the current, concentrated near the cone's vertex, has 
practically no influence on the scattering, nevertheless, the physical 
optical approach gives values for the radar thickness which are tens 
of decibels smaller than the experimental values, even with large 
dimensions of the cone. The deciding factor here is the nonuniform 
part of the current flowing in the vicinity of the sharp circular base 
rim of the conical surface; the nonuniform part of the current has an 
especially large value for sharply pointed cones. 


Another interesting example of a similar nature Is the scattering 
of a plane wave by a finite paraboloid of rotation (§ 18) where the 
physical optics approach leads to qualitatively incorrect results. 

The effective scattering area calculated in this approach turns out 
to be a periodic function of the paraboloid length, and with certain 
lengths 1t becomes zero which most certainly does not correspond to 
reality. 


The investigation of the diffraction of edie waves shows 
(Chapter V) that for flat plates one may limit oneself to considera- 
tion of secondary diffraction, if their linear dimensions are one- 
and-a-half to two times larger than the wavelen;yth. 


Let us note that we attempted te obtain equations for the 
scattering characteristics which would possess physical visualiza- 
bility and which would be convenient for makin cateculatlons. In 
keeping with this, we were obliyed to Introduce various kinds of 
interpolation equations and simplifled equations whieh satisfy the 
formulated requirements, but in return are not tn the gxeneral case 
the dominant terms of the rtyorous asymptotic expanston in powers of 
the small parameter \/a. 
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Our purpose was not to calculate the current on the body's sur- 
face, the field In the near zone, or the Iinterral seatterine thickness. 
These questions are investigated in a number of other works hased on 
the physical theory of diffraction which were already enumerated in 
§ 25. In them, in particular, the first terms of asymptotic expan- 
sions in powers of A/a were obtained for the integral thickness which 
characterizes the total power scattered by a body. !lowever, in these 
works, as a rule, equations are missing for the scattering character- 
istics wiiich are valid with any direction of irradiation and observa- 
tion. Therefore, the results of this book and the indicated works 
mutually supplement one another. 


At present, only a limited number of diffraction problems have 
ylelded to theoretical studies, as a result of which experimental 
studies of diffraction by various bodies have taken on great importance. 
In Chapter VI an experimental method was discussed which enabled one 
to Isolate in a "pure form", and to measure, the field from the non- 
uniform part of the current excited by a plane wave on a metal body 
of any shape. — In the same chapter, it was shown that the well-known 
phenomenon of depolarization of the wave reflected from a body which 
is found in free space is produced by the nonuniform part of the 
current, or, in other words, by the surface distortion. 


The investigation carried out in Chapter VIT for the problem of 
diffraction by a thin, finite length cylindrical conductor represents 
a natural development and completion of the method of consideriny: the 
multiple diffraction of edge waves which was applied in Chapter V. 

In Chapter VII equations were derived for the scattering diarram whieh 
are suitable for vibrators of an arbitrary lenrth with anv irradiitien 


and observation directions. 


The results obtained in this book show the fruitfulness of 
physical diffraction theory, and enable one to arrive at the solution 
of other more complicated problems. Such problens may be divided 
into two classes. froblems which may now already be solved on the 
basis of the known results of diffracticn theory are related to the 
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first class. As an example of such a problem, one may point to the 
osroblem of diffraction of a plane wave by a frustum of a cone or by 

an infinitely long cylinder with a polygonal transverse cross section. 
Those problems whose solution requires obtaining (using the methods 

of mathematical diffraction theory) a whole series of new results 

nust be referred to the second class. In particular, in order to give 
a complete solution to the diffraction problem of a finite cone, it 

is necessary to have more precise knowledge on the diffraction laws 

of a semi-infinite cone. 


Summing up, one may say that physical diffraction theory aids 
one in analyzing and sorting out the diffraction phenomena for complex 
bodies, poses problems for treatment by mathematical diffraction 
theory, and enables one to effectively apply the rigorous results of 
mathematical diffraction theory for the solution of new problems. 


In conclusion, I express my deep thanks to L. A. Vaynshteyn for 
his valuable advice and regular discussion of the questions to which 
this book is devoted, and also for his attentive reading of the 
manuscript and for a number of useful remarks. I also take this 
opportunity to express sincere thanks to M. L. Levin for his interest 
in this work and his helpful remarks. 
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Russian Typed Meaning 
u e evlindrical 
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